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1. - Introduction 

The recognition of the key role entanglement plays in the understanding of the radical depar- 
ture of quantum from classical physics came historically remarkably late. In the early years of 
quantum mechanics starting from the mid twenties of the last century, often referred to as the 
'golden years', this aspect was not quite in the center of activities: Researchers were occupied 
with successfully applying the new theory to a wide range of physical phenomena, continuously 
adding to the impressive list of theoretical studies matching experimental findings. It was not 
until the year 1935, when Einstein, Podolsky and Rosen expressed their dissatisfaction with the 
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state of the theory, constructing a Gedanken experiment involving a measurement setup in a 'sep- 
arated laboratories paradigm' that should identify the description provided by quantum mechan- 
ics as incomplete yj . This Gedanken experiment involved local measurements on constituents 
of a composite quantum system prepared in an entangled state in a distributed setup. In the same 
year, Schrodinger, also one of the major contributors to the theory, formulated some of the math- 
ematical implications of entanglement on the statistics of measurement outcomes; and actually 
coined the term 'entanglement' both German and in English j^l- The program envisioned by 
Einstein and colleagues - to demonstrate the incompleteness of a quantum mechanical descrip- 
tion - may be fairly said to have essentially failed. They nevertheless could pinpoint the aspect 
of quantum theory in which it would crucially depart from a local classical statistical theory. 

This was fully realized in the 1960ies, when Bell reconsidered the situation discussed by 
Einstein, Podolsky and Rosen, restated in a setting involving spin-1/2 degrees of freedom due 
to Bohm 01 • He demonstrated the vahdity of bounds to correlation functions of measurement 
outcomes of dichotomic measurements, provided that they would be resulting from a 'local re- 
alistic model', meaning from a local classical statistical theory 0|. These bounds are indeed 
violated by the predictions of quantum mechanics. After the advent of reliable sources of en- 
tangled states, many experiments were performed, all consistent with the predictions of quantum 
theory, and none with the bounds put forth in form of Bell's inequalities (see, e.g., ref. jslla]). It 
can be said that it is in the role of entanglement where the departure of quantum from classical 
physics is most manifest, indicating that the intrinsic randomness in quantum theory can not be 
thought of as resulting from mere classical ignorance in an underlying classical statistical theory. 

In the meantime, it has become clear that entanglement plays a central role also from a dif- 
ferent perspective: it can serve as an essential ingredient in applications of quantum information 
processing 0]. For example, entanglement is required for an established key to be uncondi- 
tionally secure in quantum key distribution Ist l9i ll0llllll . Entanglement is also believed to be 
responsible for the remarkable speedup of a quantum computer compared to classical computers 
[Tsiflill . the underlying logic of which being based on the laws of classical physics 

Formally, entanglement is defined by what it is not; a quantum state is called entangled, if 
it is not classically correlated. Such a classically correlated state is one that - in the distant 
laboratories paradigm - can be prepared using physical devices locally, where all correlations 
are merely due to shared classical randomness. The kind of correlations in such preparations 
are hence of the same origin as ones that one can realize in classical systems by means of com- 
municating over telephone lines. This is in sharp contrast to the situation in entangled states, 
which cannot be produced using local physical apparata alone. This facet of entanglement hence 
concentrates on the preparation procedure. The concept of distillable entanglement fill in turn 
grasps directly entanglement as a resource, and asks whether maximally entangled states can be 
extracted, distilled, within a distant laboratories paradigm. 

Part of the theoretical challenge of understanding entanglement lies in the fact that is also 
(at least partly) responsible for the quantum computational speedup: state space is big. The 
dimension of state space, the set of all quantum states corresponding to legitimate preparation 
procedures, grows very rapidly with the number of constituents in a composite quantum system. 
In fact, it grows exponentially. Consequently, in the whole development of quantum informa- 
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tion theory, it has been a very useful line of thought to investigate situations where the involved 
quantum states could be described with a smaller number of parameters, while still retaining 
the essential features of the problem at hand. So certain 'theoretical laboratories' facilitated the 
detailed investigation of phenomena, properties, and protocols that arise in quantum information 
theory, while keeping track of a set of states that is at most polynomially gro wing in dimension. 
So-called stabilizer states Werner states l23ll . matrix-product states j24r . or quasi-free 

states are instances of such 'laboratories'. In the center of this review article are the graph 
states, the structure of which can be described in a concise and fruitful way by mathemati- 
cal graphs. They have been key instrumental tools in the development of models for quantum 
computing, of quantum error correction, and of grasping the structure of bi- and multi-partite 
entanglement. 

Graph states are quantum states of a system embodying several constituents, associated with 
a graph(^). This graph may be conceived as an interaction pattern: whenever two particles, orig- 
inally spin- 1/2 systems, have interacted via a certain (Ising) interaction, the graph connecting 
the two associated vertices has an edge. Hence, the adjacency matrix of a simple graph, a sym- 
metric N X N matrix for a system consi sting of N qubits with entries taken from {0, 1}, fully 
characterizes any graph state at hand 11341 l36l l37l l38l 13911 . In this sense the graph can be under- 
stood as a summary of the interaction history of the particles. At the same time, the adjacency 
matrix encodes the stabilizer of the states, that is, a complete set of eigenvalue equations that are 
satisfied by the states (^). Thus graph states are actually stabilizer states lEiil . This class of graph 
states play a central role in quantum information theory indeed. 

To start wit h, grap h states form a universal resource for quantum computing based on mea- 
surements 1I32II33LI34II35II . In such one-way computing, one starts off with a cluster state, which 
is a specific instance of a graph state, and performs von-Neumann measurements at single sites 
associated with vertices. In fact, it was in the form of such cluster states Js^, when graph states 
have first been considered, in a narrower sense, with respect to a graph reflecting a cubic lattice. 
On the subspace that is retained any unitary can be implemented, thereby realizing universal 
computation without the need of making use of any controlled two-system quantum gates. The 
cluster state hence forms a universal resource for quantum computation. The performed mea- 
surements introduce a probabilistic aspect to the scheme; yet, the overall set-up of the one-way 
computer is deterministic, as the process can be de-randomized by means of appropriate local 
unitaries taken from the Pauli group applied before the readout step js^-ESl- 

Such one-way computing is interesting at least from two perspectives: on the one hand, it 
provides a computational model I33I different from the original gate-based model, which resem- 



(^) Note that in the literature one finds several inequivalent concepts of quantum states that are in one 
way or another associated with graphs I25l27ll . For example, entanglement sharing questions have been 
considered in a multi-partite quantum system based on quantum states defined through mathematical graphs, 
seerefs. IHIilllilll. 

(^) In some sense, this graphical representation plays a similar pedagogical role as Feynman diagrams in 
quantum electrodynamics: The latter provide an intuitive description of interaction processes in spacetime, 
but, at the same time, they have a concise mathematical meaning in terms of the corresponding propagator 
in an expansion of the scattering operator. 
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bles more closely the gate model of classical computation. In fact, questions of complexity and 
simulatability often become more transparent in the picture of the one-way computer than in 
the gate model. For instance, operations from a certain set of operations, the so-called Clifford 
operations 101, can be kept track of classically in an efficient manner. When measuring Pauli 
operators in the course of the computation, the states of the whole system will be a sequence of 
graph states that can be described by a sequence of adjacency matrices . However, if in the 
first step of the gate-based model, a non-Clifford operation is applied, unlike in the picture of the 
one-way computer, it is no longer obvious in the gate-model that the dynamics of the rest of the 
network can be described efficiently. 

On the other hand, there are good reasons to believe that the undoubtedly tremendous chal- 
lenges associated with actually realizing a quantum computer can be lessened when relying on 
an architecture based on graph states. In many quantum systems in a lattice, nearest-neighbor 
interactions are natural, rendering the preparation of cluster states through one dynamical process 
a relatively feasible step. Furthermore, and maybe more importantly, one realizes a certain dis- 
tinction between the process of creating entanglement and the process of consuming it. Hence, 
even if one exploits a lossy or even probabiUstic process in order to prepare the graph states, in 
many set-ups one can, in principle, end-up with graph states free of errors, then to be used as a 
resource of the actual computation. Even if the entangling operations are themselves faulty, e.g. 
subject to photon loss in a linear optical implementation, fault tolerant schemes can be devised, 
in principle up to a loss rate of several tens of percents. 

In quantum-gate-based quantum computation, graph states also play a prominent role as code- 
words in quantum error correction, allowing for reliable storage and processing of quantum in- 
formation in the presence of errors. This is done by appropriately encoding quantum information 
in quantum states of a larger number of quantum systems. This is the second branch how the con- 
cept of graph states originally came into play, namely in form of graph codes l4o[l4lll . These 
instances of quantum codes are determined by the underlying graph. 

Finally, the idea of the 'theoretical laboratory' itself allows for a wide range of applications. 
Aspects of bi-partite and, in particular, multi-partite entanglement are typically extraordinarily 
hard to grasp 14211 . This is again partially due to the fact that state space is so rapidly increasing 
in dimension with a larger number of constituents. The decision problem whether a mixed state 
is entangled or classically correlated is already provably a computationally hard problem in the 
dimension of the constituents. Apart from the classification of multi-particle entanglement, even 
for pure states, a complete meaningful quantification of multi-particle entanglement is yet to 
be discovered j3- Already at the origin is the question in what units to phrase any result 
on quantification matters: there is no multi-particle analog of the maximally entangled pair of 
spin-l/2-particles, to which any pure bi-partite entanglement is equivalent: any pure state can be 
transformed into such maximally enta ngle d pairs in an asymptotically lossless manner, one of the 
key results of entanglement theory 14414511 . This means that the achievable rate of distillation and 
formation is the same. The multi -partite analogue of such a reversible entanglement generating 
set (as it is formed by the maximally entangled qubit pair in the bi-partite setting) has not yet 
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been identified, at least none that has a finite cardinaHty('^). 

Within the setting of graph states, some of the intriguing questions related to aspects of multi- 
particle entanglement can be addressed, expressing properties of entanglement in terms of the 
adjacency matrix. Central tasks of interconverting different forms of multi-particle entangle- 
ment, in particular of 'purifying ' it, can and have been studied: here, the aim is to extract pure 
graph states from a supply of noisy ones, having become mixed through a decoherence pro- 
cess Isill . Based on such purification protocols, even protocols in quantum cryptography have 
been devised, making use of the specific structure of this class of quantum states 1471 ksll . In 
turn, graph states form an ideal platform to study the robustness of multi-particle entangled 
states under such in realistic settings anyway unavoidable decoherence processes. In a nutshell, 
multi-particle entangled states may exhibit a surprising robustness with respect to decoherence 
processes, independent from the system size I49I1 . Finally, questions of inconsistency of quantum 
mechanics with local classical statistical theories become very transparent for graph states Isoll . 
In all of the above considerations, the class of graph states is sufficiently restricted to render a 
theoretical analysis feasible, yet often complex enough to appropriately grasp central aspects of 
the phenomenon at hand. 

Departing sUghtly from the original formulation, weighted graph states have been considered 
1I51I1 . where the interaction is no longer a fixed Ising interaction with a constant weight. Such 
states can be thought of as resulting from a semi-classical Boltzmann gas: classical particles 
carrying a quantum spin- 1/2 degree of freedom would in an idealized description give rise to 
such a graph state through collision of particles. Such weighted graph states find numerous 
applications in the description of random quantum systems. They can also be taken as a basis set 
of states to approximate ground states of many-body systems in a variational principle. Owing 
the structural similarity to the fact that mathematically, discrete and continuous Weyl systems 
have so much in common, graph states in harmonic infinite-dimensional quantum systems have 
finally been studied that resemble very much the situation of graph states for finite-dimensional 
quantum systems [52, 53J. 

This review article aims at providing a self-contained introduction to the theory of graph 
states in most of these aspects, both conceived as a tool and as a central resource in quantum 
information theory in their own right. We will introduce the basic notions of graph states, discuss 
possible descriptions of their entanglement content, their interconversion and purification, and 
various applications in the context of quantum information science. 

I'l. Outline. - We start with a detailed introduction of graph states, whose entanglement 
properties are analyzed in the following chapters. After setting basic notations in sec. Othat 
are frequently used throughout this article, we give essentially two alternative definitions for 
graph states in sec.|5] namely, in terms of the underlying interaction pattern and in terms of the 
stabilizer We illustrate how elementary properties of a graph state and basic operations, such 
as (Pauli) measurements, on these states, can be phrased concisely in terms of the underlying 
graph. It is shown that the action of Clifford operations on graph states (sec. |3J and the reduced 



{^) For a brief review on multi-particle entanglement, see, e.g., ref. i43l . 
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states for graph states (sec.|6j can be determined efficiently from the underlying graph. These 
relations will allow for a classification of graph states in sec.Qand for an efficient computation 
of entanglement properties in sec. |8| We discuss some examples and applications to quantum 
error correction, multi-party quantum communication, and quantum computation in sec.|3| We 
briefly discuss some generalizations of graph states in the language of discrete Weyl systems in 
sec.|5] Sec.|5]contains also a short review about possible reahzations of graph states in physical 
systems. 

In sec.Qwe will then discuss the classification of graph states in terms of equivalence classes 
under different types of local operations and provide a complete classification for graph states 
with up to seven vertices. The results indicate that graph states form a feasible and sufficiently 
rich class of multi-party entangled states, which can serve as good starting point for studies of 
multi-party entanglement. 

In sec.|S|we discuss various aspects of entanglement in graph states. We briefly review some 
results about the 'non-classicality' of graph states and how their entanglement can be detected 
with Bell inequalities. The genuine multi-particle entanglement of graph states is characterized 
and quantified in terms of the Schmidt measure, to which we provide upper and lower bounds in 
graph theoretical terms. 

Finally, we introduce two possible extensions of graph states. On the one hand, in sec. |5| 
the class of weighted graph states is introduced, which comprises particles that interact for dif- 
ferent times and provides an interesting model for the study of the entanglement dynamics in 
many-particle systems. Here, we consider N initially disentangled spins, embedded in a ring 
or d-dimensional lattice of arbitrary geometry, which interact via some long-range Ising-type 
interaction. We investigate relations between entanglement properties of the resulting states and 
the distance dependence of the interaction in the limit — > oo and extend this concept to the 
case of spin gases. 

On the other hand, in secQ^I graph diagonal states serve as standard forms for mixed states 
and occur naturally whenever pure graph states are exposed to decoherence. We show that the 
lifetime of (distillable) entanglement for GHZ-type superposition states decreases with the size 
of the system, while for a class of other graph states the lifetime is independent of the system 
size. These results are largely independent of the specific decoherence model. Finally, the con- 
cept of entanglement purification is applied to graph states and possible applications to quantum 
communication are described. 

1'2. Notations. - At the basis this review article lies the concept 
of a graph I26I EtIi . A graph is a collection of vertices and a de- 
scription which of the vertices are connected by an edge. Each graph 
can be represented by a diagram in a plane, where a vertex is repre- 
sented by a point and the edges by arcs joining two not necessarily 
distinct vertices. In this pictorial representation many concepts re- 
lated to graphs can be visualized in a transparent manner In the con- 
text of the present article, vertices play the role of physical systems, 
whereas edges represent an interaction. Formally, an (undirected. 




Fig. 1. - This ring depicts a 
graph with 7 vertices 
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finite) graph is a pair 

(1) G^iV,E) 

of a finite set V = {L . . . , A^} and a set £^ C [V"]^, the elements of which are subsets of V with 
two elements each 12611 . The elements of V are called vertices, the elements of E edges. In the 
following, we will mainly consider simple graphs. A simple graph contains neither loops (edges 
connecting vertices with itself) nor multiple edges. We also regard a generalization of these 
simple graphs, where each edge {a, b} is associated with a weight (fab representing the strength 
of the respective interaction. Although the concept of a weighted graph is more general than that 
of a simple graph, we will use the notion of a graph in the more narrow sense of a simple graph, 
unless we explicitly mention that a certain section is devoted to weighted graphs. 
Since there are in general 



(2) 



fN\ _ N{N-1) 



different possibilities for choosing set of edges E' in a graph of = N vertices, the number 

of distinct graphs is 2 V 2 ; . Graph theory is mostly interested in problems that are invariant 
under permutations of the vertices, when these permutations respect the neighborhood relation, 
i.e., map neighbored vertices onto neighbored vertices. Such permutations are called graph 
isomorphisms. Two graphs Gi ~ {Vi,Ei) and G2 = (V2,i?2) are called isomorphic if there 
exists a bijection j :V\^Vi such that 

(3) {aM^Ei ^ {/(a),/(5)}ei;2 . 

Note that the number of non-isomorphic graphs still grows exponentially with the number N of 
vertices js^l- 

Vertices a,b G V that are the endpoints of an edge are referred to as being adjacent. The 
adjacency relation gives rise to an adjacency matrix Tq = T associated with a graph. F is a 
symmetric x iV-matrix, with elements 



(4) 



1, if {a,b}eE, 
otherwise. 



In the case of weighted graphs, the adjacency matrix also specifies the weights of the edges, i.e., 
Tah = (flab- We will make repeated use of the neighborhood 

(5) Na:={beV\{a,b}eE} 



of a given vertex a €z V. The neighborhood is the set of vertices adjacent to a given vertex. 
The number jA^al of neighbors is called the degree of the vertex a. A vertex a G V of degree 
\Na\ = will be called an isolated vertex. An {a, 6}-path is an ordered list of vertices a = 
ai, a2, . . . , flji-i, a„ = b, such that and aj_|-i are adjacent for all i. A connected graph is a 
graph that has an {a, 5}-path for any two a,b G V. Otherwise it is referred to as disconnected . 
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When a vertex a is deleted in a graph G, together with all edges incident with a, one obtains a 
new graph, denoted by G \ a. For a subset of vertices U C y of a graph G = ( V, ) let us denote 
with G\U the graph that is obtained from G by deleting the set U of vertices and all edges which 
are incident with an element of U. In a mild abuse of notation, we will also write G \ F for the 
graph that results from a deletion of all edges e E F, where F C E C [y]^ is a set of edges. For 
a set of edges F C [V]'^ we will write G U F := {V, E U F) and G + F {V, E + F), where 

(6) E + F = {EUF)\{EnF) 

is the symmetric difference of E and F. Note that the symmetric difference corresponds to the 
addition modulo 2 or the component-wise XOR if the sets are considered as binary vectors over 
the integer field F2 modulo two . Similarly, an induced subgraph G[A] of a graph G = (V, E), 
where A C V,is obtained by deleting all vertices (and the incident edges) that are not contained 
in A. 

Graphs may be colorable. A proper two-coloring of a graph is a labeling V — > {1,2}, 
such that all adjacent vertices are associated with a different element from {1, 2}, which can be 
identified with two colors. In graph theory these graphs are also called 'bi-partite graphs', since 
the set of vertices can be partitioned into two disjoint sets, often called sinks or sources, such that 
no two vertices within the same set are adjacent. It is a well known fact in graph theory that a 
graph is two-colorable if and only if (iff) it does not contain any cycles of odd length. 

In the remainder of this article each vertex stands for a two-level quantum system H" ~ 
or qubit. The state vector of the single-qubit system can be written as = ajO)" + /3| 1)° 
with |ap + |/9p = 1. The vectors |0) and |1) are the eigenvectors of the Pauli matrix with 
eigenvalue +1 and —1. The matrices (Tq = la, cr" = cr", erf = cr," and = cr° are the Pauli 
matrices of this two-level system, where the upper index specifies the Hilbert space on which 
the operator acts. Note that these operators form an orthogonal basis of Hermitian operators 
with respect to the scalar product {A, B) := \i{A^ B). Up to the phase factors ±1 and ±i they 
also generate the Pauli group V := ({±1, ±j} x {o-q, '^x^ o'y, o'z})- We will frequently use the 
projectors onto the eigenvectors of the Pauli operators. For example, 

denotes the projector onto the eigenvector \z, ±) of cr? with eigenvalue ±1 (similarly for cr° and 

To simplify notations, we use subsets ?7 C y as an upper index for states, operators and sets. 
They denote the respective tensor product of a given state or sets, e.g. 

(8) = (8) 1+)" or p^ = (g)7'^ 

aev aGV 

where |+) = -j= (|0) + The subsets are also used to label those vertices where the operator 
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acts non-trivially, for example 



(9) 




Moreover, we identify sets U and their corresponding binary vectors U = {Ub)b£V — [Ui , • ■ • , Un) 
over the binary field with the same symbol. Finally, a refers to both the vertex and the cor- 
responding one-element set {a} ensuring that cr^ = frl"^. These notations allow us to use set 
and binary operations in the same formula. For example, for A,B G V{V) = J-2 we will write 
A\J B, B and A\B {A V\A) for the union, intersection and difference (complement) 
as well as A + B and (A, B) for the addition and the scalar product modulo 2. 

In the multi-partite case one can group the vertices into different partitions and, for example, 
study the entanglement with respect to these partitions. Here, any tuple (^i, Am) of disjoint 
subsets Ai C V with Ai ~V will be called a partition of V . We will write 



if {Ai, ...Am) is a. finer partition than (Si, Bujr). which means that every Ai is contained in 
some Bj. The latter is then a coarser partition. 

2. - Definitions for grapli states 

With the notations introduced in the previous section we can provide some definitions for 
graph states. Throughout this article, we mainly consider two alternative descriptions. Most nat- 
urally, graph states can be regarded as the result of an interaction of particles initially prepared 
in some product state. Certainly not all imaginable interaction patterns can be represented rea- 
sonably by a simple graph. In sec. l2'll we introduce the description of graph states in terms of 
the interaction pattern and show that such a definition is also meaningful if all particles interact 
with the same Ising-type interaction but possibly for different interaction times. This description 
generalizes to so called weighted graph states, which are introduced in sec. |5| The alternative 
definition proposed in sec. I2'2I on the other hand, is restricted to the class of states that corre- 
spond to a simple graph. Such states can be described efficiently in terms of their stabilizer, 
which is a subgroup of the Pauli group. We briefly address the question of local unitary equiv- 
alence, discuss the relation to stabilizer states and illustrate an alternative representation of the 
stabilizer formalism in terms of its binary representation(*). We sketch a possible extension of 
the stabilizer formalism to d-level systems and finally summarize further alternative approaches 
to graph states in sec. 12 31 

2'1. Interaction pattern. - In this subsection we give a careful motivation for the concept of 
graph states in terms of interaction patterns, concluding with a precise definition given at the 
end of this subsection. Let us consider a set of 2-level systems (qubits) that are labeled by the 



(*) Although the remainder of the article will not be based on the binary representation. 



(10) 



{Ai,...Am) < {Bi,...,Bm'), 
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vertices of a graph G = (V, E). The qubits are prepared in some initial state vector \'^) and 
then are coupled according to some interaction pattern represented by the graph G. For each edge 
{a, 5} e E the qubits of the two adjacent vertices a and h interact according to some (non-local) 
unitary Uab = e""''"''^''''. Here, Hab denotes the interaction Hamiltonian and ipab represents the 
coupling strength or (with appropriate physical units) the interaction time. The most general of 
such setups, in which the qubits can interact according to different 2-body interactions Hab, has 
to be described by graphs, whose edges carry a labeling that specifies both the different unitaries 
Uab as well as the ordering in which interactions occur. 

Under which conditions can the outcome of this interaction pattern be completely specified 
by a simple graph G? If the graphs shall give a sufficient characterization for a large class of 
interaction patterns, we can pose the following constraints: 

(1) Since the graph G does not provide any ordering of the edges, all two-particle unitaries 
Uab involved must commute: 

(11) [Uab,Ubc\^0 ya,b,cev. 

(2) Because we deal with undirected graphs{^), the unitaries must be symmetric: 

(12) Uab = Uba ya,b&V . 

(3) If the edges are not further specified by weights, all qubits should interact through the same 
two-particle unitary U: 

(13) Uab^^^"''^ ya,beV. 

In the case of qubits condition (1) is already sufficient to restrict the analysis of particles to the 
case where the qubits interact according to the same Ising interaction, e.g. iJ^j, = <^z'^z- This 
statement is reflected in the following proposition. 

Proposition 1 (Standard form for commuting interactions). With an appropriate choice of the 
local basis in each individual qubit system, any set of commuting two-particle unitaries, i.e., the 
unitaries fulfill (1), does only contain interactions of the form 

(14) ipabHab = Vab(j''^ol + aa(T'i + abal . 

In other words, any interaction pattern in which the qubits interact according to some two- 
particle unitaries chosen from a set of commuting interactions, is up to local z-rotations{^) an 



(^) In a directed graph the set of edges E is given by ordered pairs (a, b). The order imphes that vertices 
a and b are connected by a directed edge from a to b. 

(®) I.e., y = 6*''°'^= to be performed before or after the interaction pattern. 
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Ising interaction pattern 

(15) Ui,{^ab):^e-^^^^'^>" . 

Proof: It suffices to consider condition (1) for two different unitaries U ~ e^*^ and U = 
g-i/f jjjg settings of three vertices a, 6 and c: 

(i) Uab = t/^"'^^ and C/^ = JJ^''^'^^: [ijl'''''}, i7{''^'=}] = 0, 

(ii) Uab = C/^"'''^ and Ubc = [ffl".*-}, = 0. 

Note that here H and denote the complete Hermitian generator that includes the interaction 
time or coupling strength ip. We have also used the fact that [f{A), I{B)] = iff [A, B] = 0. Ev- 
ery such Hermitian operator H allows for a real decomposition with respect to the basis of Pauli 
operators {ao, a^, a.y, a^}, i.e., H^°-'^^ = erf ct^. Moreover, a local unitary transfor- 

mation at a single qubit system translates to an orthogonal transformation of the corresponding 
operator basis cr; i— > a[, i.e.. A' = OAO^ for some orthogonal matrix O. By local unitaries 
we can thus diagonalize one of the Hamiltonians, say H^'^'^^ = Aiafa^ and represent the 
other Hermitian matrix H with respect to this basis, i.e., H^^'^^ = '^jk ^jk '^j'^k- With these 
decompositions (i) reads 

(16) ^ A.Bjfe (T,'' ® [a„ ® (7^ = , 

ijk 

from which 

(17) \fi,j = 1,2, Swithi^j : A, = OVSjfc = Vfc = 0,1,2,3 

follows. If H corresponds to a non-trivial two-body interaction up to a (local) change of basis 
we can assume that A3 ^ 0. Rewriting (ii) with these decompositions one essentially arrives at 
two different cases: If another component, say A2 7^ 0, then all components in B except Bqq 
have to vanish, which would imply that H is a trivial interaction. If instead A\ = A^ = 0, 
then at least all components in B apart from Bqq, Bq^, B^q and £^33 have to vanish. Since the 
component Bqq correspond to some negligible global phase factor, we thus have shown that any 
two commuting interaction Hamiltonians have to be of the form eq. il4\ . Any terms due to Bq^ 
or ^30 correspond to local 2-rotations and all these rotations commute with the Ising interaction 
terms i/^^^ = '^z'^z- Thus the interaction pattern can alternatively be described by an interaction 
pattern with pure Ising interaction according to the same graph and some local z-rotations to be 
applied before or after the coupUng operation. |--| 

The remainder of this article is largely devoted to the entanglement properties of states that 
result from an interaction pattern described by a simple or weighted graph. We can omit the 
2; -rotations, since they do not change these entanglement properties. In the following we thus 
consider an interaction pattern of qubits that are coupled only by pure Ising interactions. Note 
that the Ising interaction i?^^ already symmetric and hence (2) does not yield an 

additional constraint. Without condition (3) the state, which results from the application of the 
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interaction pattern, is determined by (a) the initial state vector j^*) and (b) by a weighted graph. 
This weighted graph identifies the pairs {a, b} of qubits which interact together with the in- 
teraction strength ipab (interaction time) of the respective interactions. The resuhing states are 
weighted graph states as they are introduced in sec.|5] However, in the remainder of this section 
we will restrict to states that can be described by simple graphs. Now (3) implies that we have 
to fix the interaction strength ip in eq. il5\ . For graph states according to simple graphs we will 
from now on choose if = j. Together with the choice of 

(18) I*) = 1+)^ 



for the initial state this ensures that this interaction creates maximal entanglement between to 
qubits in the state with state vector |+), i.e., t^i(,|+)|+) is maximally entangledij). In sec. l2'2l 
we will see that this choice also allows for an efficient description of the resulting states in terms 
of their stabilizer 

To further simplify notations we will not use the Ising interaction in eq. M5\ but rather the 
f controlled) phase gate 

(19) Uabi^ab) ■■= e-^'^^'"'' with Hab.= \lT{l\®\l)\\\ 

as the elementary two-qubit interaction. Note that the corresponding interaction strength now is 

'4>ab — TT, because from 

(20) H^b = = \ {lab < + Hi,) 
we find 

(21) UabK'^ab) = e 4 e -1 ^ e ■» - U^,{ — ) . 

In other words, the phase gate corresponds to the Ising interaction up to some additional ^- 
rotations around the z-axes at each qubit. For simple graphs, i.e., (^ab = tt, we find that 

(22) Uab ■■= UabM = (8 1^ + Pl_ ® (T,^ . 



This gate corresponds to a controlled cr^ on qubits a and b, i.e. 



/ 1 

10 

10 

\ -1 / 



(J) A state vector |v]>)"'' is maximally entangled iff the reduced state at one qubit is maximally mixed, i.e.. 
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The choice ipab = t' ensures not only that the state vector 



(23) 



Uab 



(10)1+)" 



is maximally entangled {Bell state) but also that U^j^ = lab or Uab = Uab ■ Consequently, the 
phase gate Uab creates as well as deletes the edge {a, 6} in a graph G depending on whether the 
edge is already contained in G or not. We conclude the above findings into our first definition for 
graph states: 

Graph states (I) 

Let G = (V, E) be a graph. The graph state \G) that 
corresponds to the graph G is the pure state with state 
vector 



(24) 



{a,b}eE 



We will also refer to the state vector \ G) of the pure state 
as a graph state. The preparation procedure reads: 

1. Prepare the qubits at each vertex in the pure state 
with state vector |+) as eigenvector of Ox with 
eigenvalue +1. 

2. Apply the phase gate Uab to all vertices a, b that 
are adjacent in G. 




Fig. 2. - The preparation procedure to ob- 
tain a graph state that corresponds to a ring 
graph with 5 vertices. 



Since Uab is the unitary two-qubit operation on the vertices a, b, which adds or removes the 
edge {a,b}, the initial state with state vector |+)^ of the preparation procedure can also be 
regarded as the graph state that corresponds to the empty graph. 



22. Stabilizer formalism. - It is often more con- 
venient to work with the stabilizer of a quantum 
state (or subspace) than with the state (or subspace) 
itself. Quantum information theory uses the sta- 
bilizer formalism in a wide range of applications. 
Among those, quantum error correcting codes (sta- 
bilizer codes) are a very prominent example |22]. 
Here, the stabilizeri^) 5 is a commutative sub- 
group of the Pauli group that does not con- 
tain — ly (and thus not ±ily). Apart from the in- 
teraction pattern, graph states can also be defined 
uniquely in terms of their stabilizer: 




Fig. 3. - The correlation operators for a graph 
state that corresponds to a ring with 5 vertices 



(*) We refer the reader to ref. 0| for an introduction to the stabilizer formalism. 
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Proposition 2. Graph states (II) Let G = (V, E) be a graph. A graph state vector \G) is the 
unique, common eigenvector in (C^)^ to the set of independent commuting observables: 

where the eigenvalues to the correlation operators Ka are +1 for all a £ V. The Abelian sub- 
group S of the local Pauli-group generated by the set {Ka \ a € V{ is called the stabilizer 
of the graph state. 

Proof: The fact that \G) is actually uniquely defined by its correlation operators, follows 
from the subsequent Proposition |3] since the set of eigenstates to all possible eigenvalues for 
Ka is a basis for (C^)^. Nevertheless, Proposition ^provides also an alternative definition for 
graph states. Hence, we have to proof that this definition is equivalent to the definition in the 
previous section. Note that the graph state for the empty graph actually is stabilized by the set 
of independent Pauli matrices {cr° | a £ V}. Hence, by induction over the set of edges E it 
suffices to show the following: Given a graph state vector \G) stabilized by Ka, the application 
of the phase gate Uab in eq. i22\ leads to a graph state vector |G") with a new stabilizer generated 
by K'a, which corresponds to the graph G' that is obtained after the edge {a,b} is added (or 
removed). This certainly holds for all vertices c G V \ {a, b}, since Kc commutes with Uab- For 
the remaining vertex b, we find 

(25) Uab Ka Uab^ - Uab (P"- + Ka ^ Ka , 

because axPz,± = Pz.^cTx- Due to Uab ~ Uba, we similarly obtain for a 

(26) UabKbUab^ ^<Kb, 

so that the transformed stabilizer corresponds indeed to a graph G', where the edge {a,b} is 
added modulo 2. q 

The generators Ka of the stabiUzer S have a straightforward interpretation in terms of corre- 
lation operators: Consider a graph state vector |G) that is measured according to the measure- 
ment pattern given by Ka — a^a^", i.e., the qubit at vertex a is measured in x-direction and 
the vertices b in Nb in ^-direction. Then Ka provides constraints to the correlations between the 
measurement outcomes m° = ±1 and ?7i^ — ±1, namely 

(27) n < = 1 • 

Since all elements a E S stabilize \G) they give rise to different constraints to the correlations 
of the spin statistics for graph states. 

That the set of correlation operators has a unique common eigenstate, is most easily seen by 
considering the graph state basis: 
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Proposition 3 (Graph state basis). Given a graph state vector \G), the set of states 

(28) \W)^aY\G) 

is basis for (C^)^. The states \ W) are the eigenstates for the correlation operators Ka according 
to different eigenvalues Wa where W = {Wi, . . . , Wn), i.e., 

(29) Ka\W)G = i-l)^- \W) . 



The projector onto the graph state has a direct representation in terms of the corresponding sta- 
bilizer S: 

(30) \G){G\^^Y.^ 

0-G5 



Proof: For the verification of eq. ( I29> it suffices to consider a single cr" operator at some 
vertex a. cr^ commutes with all correlation operators Ka for a ^ b and anti-commutes with Ki,. 
For all a £ V,we obtain 

(31) Ka\W) = Ka n ^l\G) - (-1)'"^"' n = (-1)""'° 1^) ' 

bew bew 

where 6aeu = 1 if a G U and zero otherwise. Thus, any two distinct sets W, W C V correspond 
to eigenvectors \ W) , \ W') for the set of generators {Ka \ a G V} but with eigenvalues that differ 
in at least one position a G V. Hence = Sww, where Sww = lifW = W and zero 

otherwise. Since there are 2^ possible sets, the eigenvectors {|M^)}vi/cy form a basis of (C^)^. 
A similar calculation verifies eq. OOt : 

(32) {W\ J2ct\W') = 2^SwoSw'9 

creS 

for any basis vectors \W) and \ W'). The normalization constant is due to tr(|G)(G|) = 1 
and because the number \S\ of stabilizer elements is 2^. q 

In the following we will briefly address local equivalence for the class of graph states and 
relate this class to the more general concept of stabilizer states and codes. We also present an 
alternative description of the stabilizer of a graph state in terms of its binary representation and 
review a possible generalization of the stabilizer formalism to d-level systems. 

2'2.1. Stabilizer states and codes. There exists a natural generalization of the description of 
graph states within the stabilizer formalism. Each stabilizer S, i.e., any commutative subgroup 
of the Pauli group that does not contain —ly, uniquely corresponds to some stabilized 
subspace H5 C (C'^)^, which is the largest subspace satisfying SUs = H5. The minimal 
number 



(33) 



7-5 rank(5) = min{n | ({si, ...,s„}) =5, Sj G 5} < iV 
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of generators for a stabilizer is a well-defined quantity and is called the rank of the stabilizer. 
Thus, a necessary requirement for some stabilizer S to represent a graph state is that it has rank 
A^, or, equivalently, that S is generated by N independent stabilizer elements. More generally, 
any full rank stabilizer S stabilizes exactly one (up to an overall phase factor) pure state, which is 
called a stabilizer state and which is in short denoted by \S). This stabilizer state is the pure state 
with the unique common eigenvector with eigenvalue 1 of all elements of S, which is denoted 
by S\S) = \S). Thus, every graph state is a stabilizer state; however, the class of stabilizer states 
is strictly larger than the class of graph states. 

It is clear that an A^— qubit stabilizer state vector \S) is completely determined by a set of 
N independent generators {sa}aev of S. Note that, for computational reasons, it is often much 
more efficient to deal with such a set of independent generators rather than with the complete 
stabilizer itself. However, this leads to an ambiguity in the description of a stabilizer state, 
since there are many independent generating sets for every stabilizer. Therefore, the question 
frequently arises whether two given sets of generators {sajaev and {s'^}aev generate the same 
stabilizer S. In section U'2.4l we will see an efficient approach to answer this question. 

If a stabilizer S does not have full rank r < iV, it only stabilizes an N — r dimensional 
subspaceC of (C^)^ ItI E^I . In principle, such a subspace corresponds to an [r, N ~ r]-stabilizer 
code encoding N — r into N qubits. For a decent stabilizer code the N~r degrees of freedom are 
used to detect possible errors. The main idea is to arrange the code in such a way that, under the 
influence of errors, the complete N — r-dimensional space containing the encoded information 
is mapped onto an orthogonal eigenspace of S. The coherent information encoded in this r- 
dimensional space can then be maintained by some error correction procedure. More precisely, 
suppose that some error operator a E occurs, i.e., the underlying noise process has a Kraus 
representation with a as one of its Kraus operators (^). Then if a G S the stabilized subspace 
C, and thus also any encoded quantum information, is not affected at all. On the other hand, if 
a £ \ N(5), where N(5) = {a E \ aSa^ C S} denotes the normalizer of the subgroup 
S, then (7 anti-commutes with at least one element of the stabilizer S and thus transforms the 
complete subspace C into an orthogonal subspace. By measuring a generating set of stabilizer 
elements Si the corresponding error thus can be detected. Only if the error a e is an element 
of the normalizer N(iS) but not of the stabilizer itself, i.e. a G N(iS))\iS, then this transformation 
remains within the codespace C and thus cannot be detected. More generally, for a correction of 
a set of possible errors {Ei] C , i.e., a noise process with Kraus operators Ej, the effect of 
different errors has to be distinguishable by the error syndrome obtained through measuring 
the stabilizer generators s;, i.e. SiEjC = PiEjC. One finds ItI EzI that the set of errors {Ei} is 
correctable if EjEk ^ N(iS)) \ iS for all j and k. If there is a unique error Ej associated with 
a given error syndrome Pi, the error can be corrected by applying E^^. If, however, two errors 
Ej and Ek correspond to the same error syndrome Ei, which implies EjC = E^C, both errors 
can be corrected by applying either of the operators ijj and e]^, since if Ej occurred but e\ is 



(®) Note that a restriction in the error considerations to Pauli errors is legitimate, since error correction 
capabilities of a code C can be determined w.rt. any basis of operation elements Ei (see e.g. Theorem 10.1 
and 10.2 in ref. H). 
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applied for error correction one nevertheless finds Ej^EjC = C by assumption. 

Supplementing the N — r generators of the stabilizer by r additional elements Zi , . . . , Z,. e 
to form a full rank stabilizer corresponds to the choice of a basis of codeword vectors(^°) 
in the codespace C. This basis is frequently called the 'logical computational basis'. As we 
will discuss in sec. I4'3I graph states, or more generally stabilizer states, also appear as codeword 
vectors in stabilizer codes. Together with similarly defined logical X-operators Xi, the logical 
Z-operations Zi allow for concise manipulations(^^) of the underlying code space such as error 
detection and correction, and the concatenation of codes in order to improve error correction 
capabilities. 

2'2.2. Local Clifford group and LC equivalence. Each graph State vector \G) corresponds 
uniquely to a graph G. In other words, two different graphs G ~ {V,E) and G' = {V,E') 
cannot describe the same graph state: the interaction picture tells us that \G) ~ \G') would yield 
a contradiction 

(34) 1+)^= n Uab\G')= n Uab\G) 

{a.b}&E' {a,fc}6-E' 

= n n Uab\+)^ = n uab\+v. 

{a,b}eE' {a,b}eE {a,b}eE+E' 

Here, E + E' denotes the symmetric difference eq. ^ of the edge sets, which is by assumption 
not empty and thus yields a non-vanishing interaction. 

However, graph states of two different graphs might be equal up to some local unitary (LU) 
operation. We will call two graphs G ~ (V, E) and G' — {V, E') LU -equivalent, if there exists 
a local unitary U e U(2)^ such that 

(35) \G')^U\G). 

Locality here refers to the systems associated with vertices of G and G' . Denoting 

(36) E' := USU^ = {UsU^ \ s e S}, 

where S is the stabilizer of the state vector \G), one finds that s'\G') ~ \G') for every s' E S'. 
In this sense the group S' is a 'stabilizing subgroup' of the state vector |G"), being a group of 
(local) unitary operators that have |G") as a fixed point; however, in general S' is not equal to 
stabilizer of |G"), since in general S' is not a subgroup of the Pauli group(^^). In view of this 
observation, it is interesting to consider the subclass of those local unitary operators U satisfying 



(^°) Choose \W) = \{Wi,...,Wr)) G C such that = {-1)^'\W). 

(^^ ) For details we refer the reader to ref. ItI E^I . 

(^^) This issue is closely related to the problem of local unitary versus local Clifford equivalence of graph 
states, which is discussed in sec.|7| 
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= U U\ meaning that U maps the whole Pauli group onto itself under conjugation. 
The set 

(37) CX := {U e U(2)^ | UV'^U^ - V^} 

of such unitaries is a group, the so-called local Clifford group (on qubits). If \G) and |G") 
are graph states such that U\G) = \G') for some U G CY , then the group S' in ( I36t is equal 
to the stabilizer of \G'). Therefore, the action of local Clifford operations on graph states can 
entirely be described within the stabilizer formalism - and this is one of the main reasons why the 
local Clifford group is of central importance in the context of graph states. In the following, we 
will call two graph states |G) and |G') LC- equivalent iff they are related by some local Clifford 
unitary U G , i.e., \G') = U\G). 

The local Clifford group is the iV— fold tensor product of the one-qubit Clifford group Ci 
with itself, where Ci is defined by 

(38) Ci:={U G\J{2)\UVU^ =V}. 

One can show that, up to a global phase factor, any one-qubit Clifford operation [/ e Ci is a 
product of operators chosen from the set {iJ, S*}, where 

(39) = ^1 (Hadamard gate) ~ Q ^) (single-qubit phase gate). 

The action of the Clifford group Ci under conjugation permutes the Pauli matrices cti, (T2 and 
CT3 up to some sign ±1. This can be shown as follows: First, the matrices ±cto and ±i(To are 
left unchanged under conjugation. Secondly, the set {±(Ti, ±(T2, ±173} has to be mapped onto 
itself, since U<JiU^ is Hermitian iff ai is Hermitian. Because the conjugation is invertible, the 
conjugation permutes the matrices cti, a2 and 0-3 up to some sign ±1. Also, note that it suffices 
to fix the action of U for two traceless Pauli matrices, say cti and (T2, since the action for the 
other matrices follows from linearity of the conjugation and the relation 0-3 = —i(7i(j2- 

If one disregards the overall phases of its elements, the one-qubit Clifford group has finite 
cardinality. In Tab. U we have itemized all 24 single-qubit Clifford unitaries, disregarding such 
global phases. For each unitary we have also included a possible decomposition in terms of Pauli 
operators and the ^-rotations 



(40) y±I^-e±^T-. J = 1,2,3, 

that we frequently use throughout this article. These rotations correspond to the elementary 
permutations {1,2,3} ^ {1,3,2}, {1,2,3} ^ {3, 2, 1} and {1, 2, 3} ^ {2, 1, 3} that only 
permute two indices. Instead of H and S any two of these elementary permutations can be used 
to generate the Clifford group Ci. 

An important result in the theory of graph states and stabilizer states is that any stabilizer 
state is LC-equivalent to some graph state. This statement was first proven in ref. llssi and llsill 
for the more general setup of stabilizer codes over d-level systems. 
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Table I. — All 24 single-qubit Clifford unitaries and their decomposition into elementary permutations. 



Proposition 4 (Stabilizer states). Any stabilizer state vector \S) is LC-equivalent to some graph 
state vector \G) , i.e., \S) — U\G) for some LC-unitary U € . This unitary can be calculated 
efficiently. 

Proof: A proof for the qubit case in terms of the binary framework (see sec. l2'2.4l can be 

□ 



found in ref. |3 



A similar statement holds more generally for all stabilizer codes: Any stabilizer code is LC- 
equivalent to some graph code. Thus, graph states can be regarded as standard forms(^^) for 
stabilizer states, since many properties, such as entanglement, are invariant under LC operations. 
Note that this standard form is however not unique. A stabilizer state vector \S) can be LC- 
equivalent to several graph states |Gi) = Ui\S) and IG2) = U2\S), whenever these graph states 
are LC-equivalent |Gi) = UiUl\G2). Thus, the study of local equivalence of stabilizer states 
reduces to that of local equivalence of graph states. 

Note that in general there are 24^ different Clifford unitaries (up to global phases) to relate 
two graphs states with vertices. Therefore, the difficulty to decide whether two graph states are 
LC-equivalent seems to increase exponentially with the number of parties. However in sec. 
we will briefly mention a method due to ref. 15811 that scales only polynomially with the number 
of vertices. 

Interestingly, the action of local Clifford operations on graph states can be described in terms 
a simple graph transformation rule, called local complementation Is^ : letting G = {V,E) 
be a graph and a G V, the local complement of G at a, denoted by Ta{G), is obtained by 
complementing the subgraph of G induced by the neighborhood Na of a and leaving the rest of 



(^^) In ref. IstIi some normal forms for stabilizer states are suggested that do not rely on graph states, but 
which also allow for an efficient calculation of various (entanglement) properties. 
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the graph unchanged: 

(41) Ta : G ^ Ta{G) := G + Na . 

With this notation the following result can be stated 

Proposition 5 fLC-rule). By local complementation of a graph G at some vertex a ^ V one 
obtains an LC-equivalent graph state \Ta{G)): 



(42) 

where 

(43) 



T,{G)) = U:{G)\G) 



KiG) 



is a local Clifford unitary. Furthermore, two graph states \G) and \ G') are LC-equivalent iff the 
corresponding graphs are related by a sequence of local complementations, i.e. G' = Ta-^ o . . . o 



Tq^ (G) for some ai 



, a„ g V. 



Fig-Sdepicts an example for such a successive 
application of the LC-rule. Starting with the first 
graph the complete orbit can be obtained by ap- 
plying the LC-rule to the vertices in the preceding 
graph that appear above the arrow of the following 
diagram: 





No. 1 
No. 4 
No. 7 
No. 10 



No. 2 
No. 5 
No. 8 
No. 11 



No. 3 
No. 6 
No. 9 




• Apply LC-Rule 



Fig. 4. - An example for a successive applica- 
tion of the LC-rule, which exhibits the whole 
Proof of Proposition^ equivalence class associated with graph No. 1. 

Let G be a graph with correlation operators The rule is successively applied to the vertex, 
and G' = r„(G) the corresponding graph under which is colored red in the figure, 
local complementation at vertex a with correlation operators K'^. For c e F \ Na we find 



(44) 

For b e Na, we compute 



u:k,{u:)^ = = K 



(y„ a, " 
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Thus the stabilizer S [UD^ is generated by {K'c}c£V\Nc, U {^a^'^blfcGWa- ^^Y muhiplication 
of the generators K'^K^ with K'^ (since a G V\Na) it follows that C/J S {U^Y is also generated 
by {K'^}aev and therefore stabilizes the graph state \G'). This proves that a sequence LC- 
rule applications yields an LC-equivalent graph state. That the action of any unitary within the 
Clifford group on graph states can be decornposed into a sequence of LC-rule applications is 
however more involved and we refer to refs. I39ll60ll for a proof. q 



2'2.3. Clifford group. Whereas the LC group CY is defined to consist of all local unitary 
operators mapping the PauU group to itself under conjugation, one is also often interested in the 
group of all unitary operators with this property, i.e., the group 



(45) 



Cn ■■= {u e u(2^) I ur^u^ = r^}, 



which is called the Clifford group (on qubits). By definition, Clifford operations map stabilizer 
states to stabilizer states. Up to a global phase factor any Clifford operation U can be decomposed 
into a sequence of 0{N'^) one- and two-qubit gates L7> .22l in the set {H, S, CNOT}, where H 
and S are defined as before and 



(46) 



CNOT: 



/I 0\ 

10 

1 

Vo 1 0/ 



(controUed-NOT gate) 



For example, the controlled phase gate Uab can be decomposed into Uab = H^CNOTatH^ 

2'2A. Binary representation. We now briefly 
review an alternative representation of the stabi- 
lizer formalism in terms of its binary representation 
lEjUioll . This description is frequently used in the 
literature, since it allows one to treat the properties 
of the stabilizer in terms of a symplectic subspace 
of the vector space Fj^. 

Any element u of the Pauli group can be 
represented uniquely, up to some phase factor, by 
a vector U = (U^, U^) e F|^, where U^, G 
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Fig. 5. - The generator matrix for a graph state 
that corresponds to a ring with 5 vertices. 



(47) 



n - n 



where U° , U° e F2 for every a gV . At each qubit we have used the following encoding of the 
Pauli matrices as pairs of bits: 



(48) 



(0|0), 



(1|0), 



(1|1), a, ^ (0|1) . 
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For example, the correlation operator A'l of the ring in fig.|3lhas the binary representation 

(49) A'l = cr^^crf , — , (10000|01001) e F^" . 

It is important to note that the binary representation captures the features of a Pauli operator 
only up to a phase factor 

The binary representation has the following two important properties: letting u,v,x £ 
with corresponding binary vectors U, W, X e ¥2^, one finds that 

(i) uw^ X < — > U + W = X (mod 2) , 

(50) (ii) [u, w] = < — > P W = (mod 2) , 

where ~ denotes equality up to a global phase factor and where the 2 x 2N matrix 

(51) P 














defines a symplectic inner product on the binary space F^^^. Property (i) shows that the encoding 
u G U G F^^ is a homomorphism of groups. Note that the multiplicative structure of 

the group is mapped to the additive structure of F2^, where addition has to be performed 
modulo 2. Property (ii) shows that two Pauli operators commute if and only if the corresponding 
binary vectors are orthogonal with respect to the symplectic inner product. 

It follows from (i) and (ii) that, within the binary representation, the stabilizer S of any 
stabilizer state \S) on N qubits is an N-dimensional, self-dual linear subspace S of F2^. By 
self-duaUty it is meant that 

• U^PV ^ for every U, V G S, and 



if X G F^^ and X^ PU = for every U G S, then X G S. 



The subspace S is usually presented in terms of a generator matrix (X|Z) (where X and Z 
are N x N matrices), which is a full rank 2N x N matrix, the rows of which form a basis of 
S; a generator matrix is obtained by assembling the binary representations {S]^}aev of a set 
of independent stabilizer generators {sa}aev as the rows of a 2N x matrix. Note that any 
generator matrix (X|Z) of a self-dual subspace S satisfies 

(52) (X|Z)P(X|Z)^ = 

from the self-duality of S. The generator matrix for a graph state \G) has the standard form 

(X|Z) = (i|r) , 



where T is the adjacency matrix of the graph G. Fig.|5]displays the generator matrix for the ring 
on five qubits (see also fig.|3}. 
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Choosing a different set of generators in S corresponds to a transformation of the generator 
matrix of the form 

(53) (X|Z)k^(X'|Z') = R(X|Z), 

where R is some F2-invertible N x A^-matrix. From the definition of self-duality, two generator 
matrices (X|Z) and (X'|Z') correspond to the same stabilizer iff 

(54) (X'|Z')P(X|Z)^ = . 

In section U' 2. II we encountered the problem of recognizing whether two sets of commuting 
Pauli operators are generating sets of the same stabilizer This issue has a simple solution within 
the binary representation as follows. First we must test eq. ( I54> for the corresponding generator 
matrices. In addition we have to compute the transformation matrix R = {Rat) in eq. i53\ . 
which can be achieved by Gaussian elimination over F2. Finally, we ensure that, at the level 
of the actual stabilizer rather than its binary representation, {sa}aev is indeed transformed into 
{s'ajaGV, 1-6., whether 

(55) < = n (""bf^" yaeV . 

bev 

The action of (local) Clifford operations on stabilizer states also has an elegant translation 
in terms of the binary stabilizer framework. Let U G Cjv be an arbitrary (possibly non-local) 
CHfford operation and let fu : F^^ — > F^^ be the unique function such that /(X) is the binary 
representation of UxW when x E with binary representation X G F^^. First, it follows 
from the property 

(56) Uxx'U'' = {UxU^){UxU^) 
for every x' € , that 

(57) /c/(X + X') = AKX) + /c;(X') 

for every X, X' e Fj^. In other words, fu is a linear transformation of and we write 
fuO^) ~ QX, for some (nonsingular) 27V x 27V matrix Q over F2. Secondly, the property 

[UxU\ Ux'W] = [x, x'] imphes that 

(58) X^Q^PQX' = X^PX' 

for every X, X' G F2^, showing that Q is a symplectic transformation, i.e., Q-^PQ = P. One 
can also prove the reverse statement, i.e., that every symplectic transformation can be realized as 
a Clifford operation. 
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It follows that conjugation of the stabilizer S' = U SU^ corresponds to the linear transforma- 
tion 

(59) (X|Z)k^ (X|Z)Q^ . 

Letting S, S' be full rank stabilizers on N qubits with generator matrices (X|Z), (X'|Z'), 
respectively, it is straightforward to prove the following chain of equivalent statements: 

S' = USU^ for some U Cn with corresponding symplectic matrix Q 

(60) 4=» (X'|Z') = R (X|Z) for some invertible R 

(61) ^ (X'|Z')PQ(X|Z)^ = . 

In the special case where U is a local Clifford operation, i.e., U G CY , the corresponding 
symplectic matrix Q has the following particular structure: 



(62) Q = 



A B 
C D 



where A, B, C, D are diagonal N x N matrices. The property Q-^PQ = P is then equivalent 
to AD + BC = I. This is in turn equivalent to stating that the TV 2 x 2 matrices 



(63) Q„ 



A R 



are nonsingular (over F2) for every a G V. Note that the matrices Qa correspond to the one- 
qubit tensor factors of U and, up to a simultaneous permutation of rows and columns, the matrix 
Q is equal to Qi © ... © Q^v- 

Now, suppose that G and G' are two graphs with adjacency matrices T and V, respectively. 
Then, from eq. ( 16 1> . the graph states |G) and \G') are LC-equivalent iff there exist N x N 
diagonal matrices A, B, C, D satisfying AD + BC = I, such that 

(64) r'Br + Dr + r'A + c = o. 

Thus, in order to check whether \G) and |G") are LC-equivalent, one has to decide whether the 
linear system of equations ( I64l i. together with the additional quadratic constraints AD+BC = I, 
has a solution. This approach leads to an efficient algorithm to recognize LC-equivalence of 
graph states Isiil as follows. First, note that the set V of solutions (A, B, C, D) to the 
linear equations j64> . with disregard of the constraints, is a linear vector space. A basis B = 
{bi, . . . , bd} of V can be calculated efficiently in 0{N'^) time by standard Gauss elimination over 
F2. Then we can search the space V for a vector which satisfies the constraints. As (|63 is for 
large a highly overdetermined system of equations, the space V is typically low-dimensional. 
Therefore, in the majority of cases this method gives a quick response. Nevertheless, in general 
one cannot exclude that the dimension of V is of order 0{N) and therefore the overall complexity 
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of this approach is non-polynomial. However, it was shown in ref. 15 9r that it is sufficient to 
enumerate the subset 

(65) V := {b + b' \b,b' e B} CV 

in order to find a solution which satisfies the constraints, if such a solution exists, where one 
observes that |V"| = 0{N'^). This leads to a polynomial time algorithm to detect LC-equivalence 
of graph states. The overall complexity of the algorithm is 0{N'^). We note, however, that it 
is to date not known whether it is possible to compute the LC orbit of an arbitrary graph state 
efficiently. 

2'2.5. Generalizations to d-level systems. The Stabilizer formalism can be generalized to d-level 
systems, where d is a prime power, see refs. In such a generalizations 

to systems where the Hilbert spaces of constituents are (C*, a lot of the intuition developed 
for binary systems carries over We will here sketch the situation only in which the individual 
constituents take a prime dimension. Ironically, in this more general framework, the case of 
binary stabilizer states even constitutes a special case, which has to be treated slightly differently 
than other prime dimensions. Actually, the language reminds in many respects of the setting of 
'continuous-variable systems' with canonical coordinates Iz^ . The familiar real phase space in 
the latter setup is then replaced by a discrete phase space. Also the Weyl operators, so familiar 
in the quantum optical context, find their equivalent in the discrete setting. 

At the foundation of this construction is the unique finite field Wd of prime order d. All 
arithmetic operations are defined modulo d. We may label a basis of (D^ as usual as 1 0) , . . . , | d— 1) . 
The shift operators and the clock (or phase or multiplier) operators are then defined as 

(66) U^\m) := \m + x), 

(67) Vp\m) := e'^P"'\m). 

for x,p G Wd- The number u; := e*^ is a primitive d-th root of unity. Let us assume that d is 
prime but exclude the case c? = 2. Using the above shift and clock operators, one can associate 
with each point {x,p) £ in phase space a Weyl operator according to 

(68) wix,p) ■.= VpU^. 

These Weyl operators correspond to translations in phase space. In analogy to the previous 
considerations, we may define the elementary operators 

(69) X^Ui, Z = Vi, 
satisfying X"^ ~ 1 and Z"^ = 1. The operators 



(70) 



'P'^wix^p) 
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for (.T, p, v) £ form a representation of the Heisenberg group with its associated group com- 
position law. The Weyl operators in this sense can be conceived as generalized Pauli operators 
familiar from the binary setting. They satisfy the Weyl commutation relations 

(71) w{x,p)w{x' ,p') — Lu^P '^w{x + x\p + p'), 

as can be readily verified using the above definitions, so the product of two Weyl operators is 
up to a number again a Weyl operator. This is the discrete analog of the familiar canonical 
commutation relations for position and momentum for Weyl operators in continuous phase space, 
which takes essentially the same form. It follows that two Weyl operators w{x,p) and w{x' ,p') 
defined in this way commute if and only if 

(72) =0, 

so if and only if the standard symplectic scalar product vanishes, which is defined as 

(73) [(x,p), (x',p')] := i^.P) f ? 7>M (^''P')^ = P^' - V 



.1 . 

for {x, p) G F^. This is a antisymmetric bi-linear form in that 

(74) [{x,p),ix',p')] = ~[{x',p'),ix,p)]. 

Hence, the discrete phase space is a symplectic space over a finite field. In turn, the linear 
combinations of all Weyl operators form an algebra, the full observable algebra of the system. 

The composition of N constituents of a composite systems, each of dimension d, can be in- 
corporated in a natural fashion. We now encounter coordinates in phase space {xi,pi, xn,Pn) 
with X G F^ and p G F^. The above symplectic scalar product is then replaced by the one 
defined as 

(75) [(x,p),(x',p')] = (x,p).a.(x',p')^, 
with 



N 

(76) <^ := 



-1 

1 



Similarly, the Weyl operators become 

(77) w(x,p) = w{xi,pi) (g) ... (g) w{xn,Pn), 
and let us set 

(78) Win, p) ;= cj-T^w(x, p). 
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We now turn to the actual definition of stabilizer codes and stabilizer states 

SEilEi- At 

the foundation here is the definition of an isotropic subspace. An isotropic space S C is a 
subspace on which the symplectic scalar product vanishes for all pairs of its vectors, so where 

(79) [(x,p),(x',p')] =0 

for all (x, p), (x', p') e S. Now, let x ^ character as a map from S into the circle group (for 
example, the map mapping all elements of S onto 1). Let us denote the dimension of S with 
k. Then, the projector onto the stabilizer code associated with the isotropic subspace S and the 
character x can be written as 

(80) p = ]k x*(x,pmx,p). 

In particular, the state vectors \tp) from this stabilizer code are exactly those that satisfy 

(81) x*(x,p)M^(x,p)|^) = ItA), 

for all (x, p) £ S. In other words, this state vector is an eigenvector of all of the operators 
X* (x, p)VF(x, p) with the same eigenvalue +1, as we encountered it in the binary setting. Again, 
it is said that is stabilized by these operators. The above Weyl operators are, notably, no 
longer Hermitian. Hence, they do per se allow for an interpretation in terms of natural constraints 
to the correlations present in the state. This setting can be naturally been generalized to prime 
power dimension d = p"^ with p being prime and r being an integer If, however, the underlying 
integer ring is no longer a field, one loses the vector space structure of F^, which demands some 
caution with respect to the concept of a basis (^^) 

Similarly, a stabilizer code can be conceived in this picture as the image of an isotropic 
subspace S under the Weyl representation. If a stabilizer code is one-dimensional, it is a stabilizer 
state. 

Again, any stabilizer state can be represented as a graph state, up to local Clifford operations. 
This has been shown in refs. jsslE^I- The notion of a Clifford operation still makes sense, as a 
unitary that maps Weyl operators onto Weyl operators under conjugation, 

(82) [/iy(x,p)[/taW^(Q(x,p)), 

where Q is an element of the symplectic group, so preserves the above symplectic form. 

The respective graph state corresponds to a weighted graph with weights Fab G F^, with a, b 
again being associated with the vertices of the underlying graph. A graph state is now a state 



(^'') If d contains multiple prime factors the stabilizer, consisting of d'^ different elements, is in general no 
longer generated by a set of only A'^ generators. For tlie minimal generating set more elements N < m < 
2N of the stabilizer might be needed l63ll . 
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stabilized by the operators 

(83) Ka = Ut n {V^f-^ =^"11 (^')'^"- 

The symmetric adjacency matrix F contains elements Tab G and, thus, it has no longer binary 
entries as in the case of a simple graph as in qubit systems. The interaction is instead specified by 
a strength r ~ Tab given by the weight of the edge {a, h] in the weighted graph. Note, however, 
that this concept of a graph state based on a weighted graph is different from the one used in 
the remaining part of this review article (see sec. When conceiving the preparation of the 
graph state via the successive application of phase gates, the associated unitary J/^j, acting on the 
Hilbert spaces of the systems labeled a and h is given by 

(84) Ul,\mr\n)<' = oj-^^^\mr\n)\ 

This picture of graph states in discrete Weyl systems, embodying the case of d-dimensional 
systems, as well as their processing in the context of the one-way computer, has been considered 
in detail in refs. (illl . Also, quantum error correcting codes have in this setting been described 
in ref. j^. This language of discrete Weyl systems provides a clear-cut picture to describe 
finite-dimensional systems in phase space. 

2'2.6. Remarks on harmonic systems. Finally, it is worth noting that the close analogy between 
discrete and continuous Weyl systems suggests the existence of similar structures as graph states 
in the setting of quantum systems with canonical coordinates, so systems in a real phase space 
with position and momentum coordinates. Variants of such an idea have been considered in a 
number of publications Is^ IsSL l65i. i661 : to describe them in detail, yet, would be beyond the 
scope of this review article. Here, we rather note the structural similarities to the previous finite- 
dimensional setting. 

The phase space of a system with N canonical degrees of freedom - N harmonic oscillators 
- is R^^, equipped with an antisymmetric bi-linear form defined by 

This form originates from the canonical commutation relations between the canonical coordi- 
nates of position and momentum, which can be collected in a row vector as (ri, ■■■,r2N) = 
{xi,pi, xtvjPat). These canonical coordinates satisfy the canonical commutation relations 
between position and momentum, although these variables can, needless to say, correspond to 
quadratures of field modes. As before (in this now real phase space) we may introduce Weyl 
operators embodying translations in real phase space, defined as 



(86) 
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for (x, p) G R^^. This Weyl operator is, in a number of different conventions, a frequently used 
tool in quantum optics under the name of displacement operator. These Weyl operators inherit 
the canonical commutation relations: it is easy to see that they satisfy the Weyl relations 

(87) W{x, p)W{x', p') = e-'('''P)'^(^''P') VF(x + x', p + p'). 

The structural similarities are obvious. The characteristic function is here, just as in the discrete 
case, defined as the expectation value of the Weyl operator, 

(88) /(x,p) =tr(W^(x,p)p) 

for (x, p) G R^^. This is a generally complex-valued function in phase space, uniquely defining 
the quantum state. Hence, the description in terms of Weyl systems serves also as a language 
appropriate for the description of both the discrete and as well as the infinite-dimensional setting. 

A certain class of states for which the assessment of entanglement is relatively accessible is 
the important class of Gaussian states. They are those quantum states for which the character- 
istic function is a Gaussian. Then, the first moments, dj = tr{rjp) and the second moments 
fully characterize the quantum state. The second moments, in turn, can be embodied in the real 
symmetric 2 x 2Af -matrix 7, the entries of which are given by 

(89) 7j- fc = 2Re tr ( {r, - dj ) (r^ - 4 )p) , 

J, fc = 1, iV. This matrix is typically referred to as the covariance matrix of the state. Simi- 
larly, higher moments can be defined. 

Analogues or 'close relatives' of graph states in the Gaussian setting now arise in several 
context: (i) They can be thought of as orig inating from an interaction pattern, similar to the 
interaction pattern for Ising interactions These interactions may arise from squeezing and 
Kerr-like interactions, (ii) They can also be resulting as ground states from Hamiltonians which 
are specified by a simple graph, in turn reflecting the interaction terms in the Hamiltonian 

(90) F = pp^/2 + xVx^, 

where againp = (pi, ...,pAr), x = (xi, xjv), and the real symmetric matrix A^ x A^-matrix V 
incorporates the interaction pattern as the adjacency matrix of a weighted graph Eilillillil. 
Then, the resulting covariance matrix is nothing but 7 = V^^/^ © V^/^, when ordering the 
entries in the convention of {xi, ...,xn,Pi, ■■■,Pn)- (iii) Also, the direct analog of stabilizer state 
vectors (conceived as state vectors 'stabilized by a stabilizer') in the setting of continuous Weyl 
systems still makes sense, yet one has to allow for singular states l63l which can no longer be 
associated with elements of the Hilbert space of square integrable functions, but can conveniently 
be described in an algebraic language (or within a Gelfand triple approach). 



Entanglement in Graph States and its Applications 



31 



2'3. Alternative approaches. - Due to the description in terms of their stabilizer, graph states 
can be represented, under suitable interpretations, by various mathematical structures, which 
connect these objects also to other areas of applications in classical cryptography and discrete 
mathematics. For example, a graph code can be described by a self-dual additive code over the 
field ¥4 = GF{4:) or by a (quantum) set of lines of the projective space over ¥2 = GF{2) 
i6o[l69ll . Graph states are also equivalent to quadratic boolean functions Foll . which are used in 
classical cryptography. 

In the remainder of this section we focus on another 
description of graph states in terms of Valence Bond 
Solids (VBS), which does not rely on the stabilizer for- 
malism and, hence, can be extended to weighted graph 
states (see sec.|5J. This representation was recently in- 
troduced by Verstraete and Cirac iFlll and has already 
found interesting applications in density-matrix renor- 
malization group (DMRG)(^^) methods. The VBS pic- 
ture has its roots in the Affleck-Kennedy-Lieb-Tasaki 
(AKLT) model iFsll . which allows to find exact expres- 
sions for the ground states or exited states of some par- 
ticular Hamiltonians. In DMRG variational methods are tj- • u ^ . <■ -^u c 

Fig. 6. - A graph state tor a ring with 5 

applied to a generalization of these AKLT-states, the so- vertices as valence bond solid, 
called matrix-product states I24I1 . in order to perform numerical studies of various physical sys- 
tems, especially within the field of condensed matter physics. Lately, VBS states have attracted 
some attention, since they allow for a clear reformulation of DMRG algorithms leading to im- 
provements of the DMRG methods for the simulation of many-body systems in two or higher 
dimensions or with periodic boundary conditions Fsll . 

In the context of this article graph states can also be regarded as particular VBS states: Here, 
the graph state \G) arises from a set of Bell pairs (bonds) 




(91) \Br^'' =Ua.bA+r^+f' 

between some virtual qubits after some suitable projections P onto the real qubits (see fig.|^: 

(92) \G)=\{Pa n 1^)"'"' 

{a.,b)eE 

More precisely, the graph state can be obtained following the procedure: 

1. Replace the real qubits at each vertex a by da virtual qubits a^, . . . , a''", where da = \Na\ 
denotes the degree of the vertex a. 

2. For each edge {a, h] in G create a Bell pair \B)'''^' between some virtual qubit a* at vertex 
a and some virtual qubit V at vertex h by using the Ising interaction U . 



(^^) For a review of these methods we refer the reader to ref. It^I . 
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3. Project all virtual qubits at each vertex a into the real qubit (sub-)system by 

(93) Pa :="io)(or\..(or'" +''|i)(ir\..(ir'° . 



That this procedure provides an equivalent description for graph states, can be shown inductively 
using the fact that the phase gate Uaibj on the level of the virtual qubits 'commutes' with the 
projection onto the real physical qubits, i.e. 

(94) [Ua^t^,Pc]=0 Va^6^c= (c\...,c'^=) . 

For the twisted four-qubit ring this is depicted by the commutative diagram in fig.El 



3. - Clifford operations and classical simulation 



The stabilizer formalism is not only 
suited to describe states (or codes), but 
also to calculate the action of Clifford 
operations on these states. Clifford op- 
erations are (possibly non-local) Clifford 
unitaries U £ Cn (see eq. ( I37t ) and pro- 
jective measurements of a Pauli opera- 
tor s G , which we will call Pauli 
measurements. The restriction to projec- 
tive measurements in the Pauli basis en- 



U 




P 



P 



u 




Fig. 7. - The phase gate on the level of the virtual 
qubits 'commutes' with the projection onto the real physi- 
sures that such measurements performed cal qubits. 

on stabilizer states [codes] yield again stabilizer states [codes] as measurement results iTi Eiil . 

It is not necessary to consider measurements of arbitrary Pauli operators s G . Since it is 
possible to efficiently decompose an arbitrary Clifford unitary U G Cjv in terms of the one- and 
two-qubit gates H, S and CNOT (see sec. l2'2.2V any Clifford operation can be simulated by a 
sequence of at most 0{N'^) of these gates together with one Pauli measurement (say s = a^) 
at a single vertex. Thus a sequence of Clifford operations acting on some stabilizer state can 
be replaced by a sequence of one- and two-qubit gates H, S and CNOT and single-qubit Pauli 
measurements with only a polynomial overhead 0{N'^) in the number of gates. In the circuit 
model for quantum computation an equivalent scheme is often considered. The class of quantum 
computations that involve only 



state preparations in the computational basis. 



• the one- and two qubit gates H, S and CNOT and 



• measurements of observables in the Pauli group V, including the classical control of gates 
conditioned on the outcome of such measurements. 
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is called the class of stabilizer circuits. All the states of the 'quantum register' in each step of 
such a stabilizer circuit are stabilizer states. These states can be characterized by their set of sta- 
bilizer generators. A formal representation of this set of generators in the memory of a classical 
computer(^^) allows one to efficiently keep track of all changes by pure classical computation. 
The effect of the one- and two-qubit gates as well as the one-qubit Pauli measurements to the 
generating set can be calculated using 0(A^^) (^^) steps on a classical computer. In this way, any 
Clifford operation can be efficiently simulated on a classical computer, which is the content of 
the Gottesman-Knill theorem ll7l l76ll . 



Proposition 6 f Gottesman-Knill theorem). Any stabilizer circuit on a quantum register of N 
qubits, which consists of M steps, can be simulated on a classical computer using at most 



0{N^M) elementary classical operations. 



□ 




Although this result was already known for 
a few years , only very recentl y, su ch clas- 
sical simulator was implemented ITTIl that ac- 
tually requires only 0{N'^) elementary opera- 
tions on classical computer In the remainder 
of this section we will see how graph states 
can provide an alternative algorithm. This al- 
gorithm is based on elementary graph manip- 
ulations and was proposed and implemented 
in ref . jjil . Its complexity for the elementary 
gate operation requires 0{(P) basic steps on a 
classical computer, where d denotes the maxi- 
mal degree of the graph representing the quan- 
tum register This proposals is hence advan- 
tageous if this maximal degree remains com- 
parably small to for the different register 
states throughout the computation. 

According to sec. l2'2.l1 anv stabilizer state \S) can be represented as a graph state \G) up to 
some LC-unitaries U G CY . Thus, in order to keep track(^^) of the different steps i = 1, . . . , Af 
in the stabilizer circuit computation, one has to store a local Clifford unitary Wi as well as the 



Fig. 8. - Example for a Ui, -measurement at vertex 1 in 
graph No. 1, which is followed by a cr^ -measurement 
at vertex 2: In graph No. 1 a cr^ -measurement is per- 
formed at the vertex 1. For the application of the x- 
measurement rule, vertex 2 was chosen as the spe- 
cial neighbor ho, yielding the graph No. 2 up to a lo- 



cal unitary C/^_4 = (±*'^h ) ■ As stated in Table 
nil the subsequent -measurement on the new graph 
state is therefore essentially another CTj; -measurement, 
now at vertex 2 with a single neighbor bo = 5. The 
final graph is then graph No. 3. 



( ) In the binary representation a computer has to store the generator matrix and additional phases at each 
qubit. The matrix requires 0{N^ ) memory size, whereas for the phases a register of size 0{N) is sufficient. 
With this information the complete stabilizer can be recovered (see sec. l2'2.lV 

(^'^) Note that the update of the stabilizer can be determined in only 0{N'^), but the determination of the 
exact measurement result in the case of measuring a Pauli-matrix ±(Ti G S seems to require some Gaussian 
elimination, which needs 0{N'^) time in practice fTTil . 

(^*) Usually a stabilizer circuit is required to start in some kind of standard input state |0)^ of the compu- 
tational basis, which has a trivial representation in terms of graph states, i.e., |0)^ — H^\-\-)^ . But the 
following argumentation will also hold for an arbitrary stabilizer state as the input, if one allows for a poly- 
nomial overhead at the beginning of classical simulation in order to determine the corresponding graphical 
representation. 
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Px,±crz 


— (^zPx,^^ 




~ (^zPy,^, 


Pz,±az 


= '^zPz,±, 


Px,±{-iazY'^ 


= {-ia^f/^Py,^, 


Px,±{iayf''' 


= {iOyY^^Pz,±, 


Px.±{-i<TyY'^ 




Px.±{iOzf/^ 


= {ia,f/^Py,±, 



Table II. - The relevant commutation relations for 
of Pauli measurements is applied to a graph state. 



Py,±{-ia^ 


1/2 ^ 


{-iTzY'^Px,±, 


Py,±{icry 


1/2 ^ 


{tayy/''Py,±, 


Py,±{-iOy 


1/2 ^ 




Py,±(i(7z 


1/2 ^ 




Pz,±(—icrz 


1/2 ^ 




Pz,±{iO-y 


1/2 _ 


{iay)^^^Px,±, 


Pz,±{-i(Ty 


1/2 = 


{-'iayY''^Px,±, 


Pz,±{i(Tz 


1/2 = 


{ia.Y'^P.,±, 



i projections and Clifford operators if a sequence 



graph Gi of the graph state \Gi), which is LC-equivalent to the actual stabiHzer state \Si) in 
the quantum register of step i, i.e., \Si) = Wi\Gi). Note that the storage of a graph Gi on 
N vertices requires only bits for the entries of the corresponding adjacency matrix. 

Moreover as discussed in sec. l2'2.l1 at each vertex a the list of single-qubit LC-unitaries Wi can 
be characterized by one of 24 'permutations' of the Pauli matrices depicted in TableU 

In order to be an efficient representation for the classical simulation of the stabilizer circuit, 
the graphical description remains to be provided with a set of graphical rules that account for the 
changes of the stabilizer when a one- or two-qubit gate or some single-qubit Pauli measurement 
is applied to it. The one-qubit Clifford unitary occurring at a vertex can easily be dealt with by 
updating the corresponding unitary according to some fixed 'multiplication table'. For the two- 
qubit unitaries we will, instead of the CNOT gate, consider the phase gate Uab in eq.|22j since on 
'pure' graph states it simply acts by adding or deleting the corresponding edge {a, 5}. However, 
the case where Uab does not act directly on \G) but on M^|G) for some LC-unitary W that is 
non-trivial at the vertices a and h, i.e., ^ la or ^ Ih, requires a more careful treatment. 
Remember that each of the possible single-qubit unitaries has a decomposition in terms of ele- 
mentary ^-rotations given in Table |I] Since all unitaries Wa = la, z'^jX'^ jU"" , y/Haz can be 
'commuted through' the phase gate yielding at most some additional ct^ on the vertex a or b, e.g. 
Uab<^x = cTx^^Uab, wc are left with the analysis of the 14 additional cases, for which at least one 
unitary Wa or Wh is of the type ^/Sa^ or ^J±i(7y. The graphical rules for these cases can be 
obtained using the LC-rule derived in sec. l2'2.2l in order to remove these unitaries. For example, 
one finds that for an arbitrary graph G 

(95) Uab Vii^lG) = x/t^^"^" UabHG)) , 

since = \f~%a~J^°' {UJ^' (similarly for \/ —ioj^'), where U^ denotes the LC-unitary in 

eq. ( I43> for the LC-rule. In fact in ref. ItsIi it is shown that, in this way, any of the remaining LC- 
unitaries Wa at some vertex a can be removed by means of at most five local complementations 
r applied at this vertex or at one of its neighbors. 

Let us finally examine the effect of single-qubit Pauli measurements in more detail, since 
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the graphical rules will be used in the subsequent sees. We will at first consider the case of a 
projective measurement of some Pauli operator ctx, ay or at a singe vertex a in a graph state 
without additional LC-unitaries at this vertex and will later mention how to cope with the general 
case. For a Pauli measurement of the graph state |G) at a vertex a we find that the graph \G') on 
the remaining unmeasured vertices can be obtained from the initial graph G by means of vertex 
deletion and local complementation: 

az ■ deleting the vertex a from G; 

(Ty : inverting G[iVa] and deleting a; 

ax ■ choosing any bo G Na, inverting G[Nbg], applying the rule for ay and finally inverting 
G[Nba] again. 



This is the content of the following proposition 137116 111 . 

Proposition 7 fLocal Pauli measurements). A projective measurement of ax, ay, or a^ on the 
qubit associated with a vertex a in a graph G yields up to local unitaries Uf^ a new graph state 
\G') on the remaining vertices. The resulting graph G' is 

(96) P;_±|G) = ^ \z, ±T ® K±\G - a)> 

(97) P;,±|G) = \v, ±r ® U'^MG) - a), 

(98) Pl^\G) = |x, ±r ® C/:,±|t,„ (t. o (G) - a)) , 

for any choice of some bo G Na, whenever the ax-measurement is not performed at an isolated 
vertex. If a is an isolated vertex, then the outcome of the a^-measurement is +1, and the state is 
left unchanged. The local unitaries t/f are 



(99) 


^ z,+ — 


1, 


VI- 




(100) 









1 : 

= V+iCT^ 


(101) 


Ux.+ = 


-. b„ NaMNigUbo] 

-lay az 


\ K- 


= ^J-iay az ° 



For a measurement of ax the local unitary Ux.± depends on the choice of bo. But the resulting 
graph states arising from different choices of bo and b'o will be equivalent via the LC-unitary 

"o bo J- 

For a sequence of local PauU measurements, the local unitaries have to be taken into account, 
if the measured qubit is affected by the unitary. We have summarized the necessary commutation 
relations in Table HII which denote the transformation of the measurement basis, if a subsequent 
measurement is applied to a unitarily transformed graph state. Fig. |8] shows two subsequent 
applications of the rather complicated (T^; -measurement. An exhaustive table can also be provided 
for the general case that the Pauli measurement occurs at some vertex a that has an arbitrary non- 
trivial LC-unitary Wa attached to it. 
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Proof {of Proposition0: The -measurement rule follows directly from the definition of 
graph states in terms of the underlying interaction pattern (see eq. MAV ): 

(102) pi±\G)^pi± n n '^-'^1+)"' 

{a.b}eE i<^.d)eE 

1 J |z,+)°(8)|G \ a)^\° if measurement result is to" = +1 

\ |z, —Y ® cr^° |G \ a)^\° if measurement result is to" = -1 

In other words, with probability i a cr^ -measurement at a vertex of some graph state gives either 
|G \ a) as the graph state on the remaining vertices if the measurement outcome is m.^ = +1 or 
(T^" \ G\a) as the graph state on the remaining vertices if the measurement outcome is to^ = — 1 . 

With the LC-rule at hand (see Proposition|5} one can now derive the measurement rules for 
a Gx- or cTy-measurement from this cr^ -measurement rule. For this, one can use commutation 
relations, which are similar to those in TableHl] in order to show that 

(103) Pl^ = Ul (G)Pl^ {VI (G))t 

(104) P;,± = [/:(G)P,"^([/:(G))t 

where 6o is a neighbor(^^) of a in G. With eq. ( I104l i and using e^*^'^'' = \/±iak = 
can now compute 

Pl^\G)^Ul{G)Pl^\Ta{G)) 

= jr(r\^ / ® ^z^\^a{G) \ a)^" If = +1 

"^^V2 I \z,+r®\Ta{G)\a)''\- ifm^ = -l 

±_ j ® Vi^a^^kaiG) \ a)^\" if = +1 



(105) 



V2 \ V^k, +)" ® Vi^haiG) \ a)^\" if ml 

J_ i |2/,+)"® v/-»fTf"|T,(G)\a)^\" ifm;; = +l 
V2 [ |y,-)"® v/i^|r„(G)\a)^\" if m" = -1 



-1 



This is the CTy -measurement rule, from which the CTx -measurement rule now can be derived along 
the same lines using eq. (I103> . q 

4. - Examples and applications 

In this section we give some prominent, not necessarily distinct classes of examples for graph 
states. We also sketch important applications of these states in multi-party quantum commu- 
nication, quantum computation and quantum error correction. These examples illustrate that 



(^^) Note that if a is an isolated vertex the graph state is \G) — |+)" IG')^^" for some graph on the 
remaining vertices. In this case a era; -measurement yields +1 with probability 1. 
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graph states do not only provide an interesting model to study multi-party entanglement, as it 
will be carried out in the following sections, but can also be an important resource for quantum 
information processing. 



4' 1 . GHZ-states. - We start by considering the 
A/^-qubit Greenberger-Horne-Zeilinger states 



(106) 



1 



\GHZ) = ^(|0)«^ 



\®N\ 




which were introduced in ref. I79II for the case of 

three qubits and since serve as a 'text book' ex- Fig- 9- - The GHZ-state is LU-equivalent to the 



graph state corresponding to a star graph or the 
complete graph. 



ample for multi-party entangled 'Schrodinger cat' 
states. These states are special examples of states 
that maximally violate multi-partite Bell inequalities llsoll and can, for instance, be used to im- 
prove frequency standards 118 ill . GHZ states have also become an interesting resource for multi- 
party quantum communication, e.g. in the context of secret sharing and secure function evalu- 
ation I82I1 . The multi-party GHZ-state corresponds to the star graph and the complete graph. 
This is easily seen by applying Hadamard unitaries H^^"" to all but one qubit a in the GHZ- 
state, which yields the star graph state with a as the central qubit. A further application of the 
LC-unitary for the LC-rule Ta in Proposition [S] then transforms the star graph state into the 
complete graph state. Thus the star graphs for different central vertices as well as the complete 
graph are LC-equivalent representations of the GHZ-state. 

4'2. Cluster states and the one-way quantum computer. - The initial resource f or q uantum 
cornputation in the one-way quantum computer (QCg), as it was introduced in refs. j32ll33ll34i 
I35L 13611 . is the cluster state in two dimensions, which corresponds to a rectangular lattice. More 
generally, a d-dimensional cluster state is represented by a graph of the form of a d-dimensional 
lattice. But only cluster states of dimension two and higher can serve as a universal resource for 
quantum computation Issll . 

In the following we consider the example of the quantum Fourier transform (QFT) on three 
qubits as depicted in fig. ^| to sketch the basic concept of the QCq model. In the standard 
framework of quantum computation (^") this QFT is implemented by performing a sequence of 
elementary (2-qubit) gates operating on the 71 = 3 input qubits. In the QCq model we start 
instead with the preparation of a highly entangled resource state via the same entangling pro- 
cedure as for the cluster state, except that a subset of n of these qubits representing the input 
qubits initially are prepared in the desired input state. Then the actual computation consists in 
a sequence of local measurements which turns the state on the remaining, non-measured qubits 
into the desired (pure) output state. 

In general, the measurement outcomes of the local measurements have to be stored in a classi- 
cal register ('information flow vector') of 2n bits, since the direction of subsequent measurements 



{^^) For details we refer e.g. to ref. 0|. 
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Fig. 10. - -Quantum Fourier Transformation on 3 qubits- The upper left graph is a 13 x 5-cluster state 
and represents the initial resource for implementing a QFT on three qubits in the framework of the one-way 
quantum computer. After removing some superfluous vertices by means of -measurements the right up- 
per graph is obtained, whereas the lower graph is achieved after performing all Pauli measurements within 
the local measurement protocol associated with the QFT, except the measurements on the input and output 
vertices and corresponds to the 'non-classical' resource for the QFT. The form and color of the vertices 
indicates the different measurement directions involved or whether the respective vertex corresponds to an 
input or output qubit. 



have to be adapted according to these results Moreover, each measurement 

round requires an update of the information flow vector that contains the 'algorithmic infor- 
mation', from which, at the end, the outcome of the quantum computation(^^) can be read off 
directly. The different dependencies of the measurement directions upon previous measurement 
results induce a temporal ordering on the set of vertices that can be measured simultaneously, 
whose total length represent the temporal complexity of the quantum computation in the QCq 
model. In a first measurement round all Pauli-measurements can be performed at the same time, 
since these do not depend on previous measurement results. This fraction of the quantum com- 
putation corresponds to the part of the network model that solely consists of Clifford operations. 
According to the Gottesmann-Knill theorem(^^) this Clifford part can be efficiently simulated on 
a classical computer Note that stabilizer circuits are not even universal for classical computation 

El. 

In the framework of the QCq the stabilizer circuits can be dealt with by using only a single 
measurement round. If the input state is some stabilizer state, e.g. some state of the compu- 
tational basis, then such a Pauli-measurement round still ends up in a stabilizer state, whereas 
the states obtained after some subsequent measurements can in general no longer be described 



(^^) E.g. as it would have been obtained by a corresponding network of quantum gates. 
(^^) See Proposition|6|in sec.|3| 
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within the stabilizer formalism. Note that the first measurement round might as well contain the 
input and output(^'^) vertices, which in general are measured in cr^- and tr^-direction if the quan- 
tum computation is carried out with respect to the computational basis. In other words, input and 
output qubits can also be measured long before the whole quantum computation is completed, 
since then the actual result of the quantum computation is given by the information flow vector 
after the last measurement round. 

In fig. ^1 we have depicted the graphs that are obtained after performing all measurements 
and after all Pauli measurements except those at the input and output vertices. The corresponding 
graph states can still be used to implement the QFT. For this the input state is teleported(^^) into 
the input vertices and the usual measurement protocol is applied. Alternatively one can regard 
the graph as a preparation procedure for an initial resource that already incorporates the input 
state. More precisely one prepares the input qubits in the input state and the remaining qubits 
in the |+) state, entangles the qubits according to the interaction pattern given by the graph and 
finally carries out the measurement protocol as described above. 

In I3^l33[l34ll3^ it is shown that any quantum algorithm within the network picture requiring 
only a polynomial(^^) amount of temporal, spatial and operational resources can be efficiently 
simulated by a QCq that requires (i) a polynomially overhead of elementary operations in the 
classical preprocessing to derive the equivalent setup and measurement protocol for the QCq, (ii) 
a polynomially bounded amount of classical and quantum resources, and finally (iii) a polyno- 
mial increasing time cost for the classical and quantum processing during the computation (^^). 
In this way the QCq can serve as a universal model for quantum computation that has quite 
promising scaling behavior for those practical implementations, in which the resource cluster 
state or more generally the initial graph state can be prepared by a homogeneous Ising interac- 
tion, i.e., independently of the system size. Furthermore, the QCq is equivalent l7lll84ll to other 
measurement-based schemes for quantum computations. 

First results towards fault-tolerant quantum computation with the QCq were also obtained 
in refs. For a reasonable noise model, including noisy cluster state preparation 

as well as erroneous measurements, a quantum computation subjected to decoherence below a 
certain threshold can be implemented on a QCq in a fault-tolerant way. This means that for 
sufficiently small noise of a specific type the quantum computation can be implemented on a 
QCq with 'any' desired accuracy if one allows for a reasonable overhead in the computational 



(^^) If the objective of the quantum computation is not only to compute a (coherent) quantum state on the 
output qubits but also to perform a read out in order to obtain some classical result. 

(^'') Note that this teleportation can be described within the QCq picture. For this the qubit holding the 
input state at qubit ai, some auxiliary qubit 02 in the |+) state and the qubit at the input vertex b are 
entangled to form a chain attached to the graph at b . Then a\ and 02 are measured in era; -direction. The 
resulting state on the graph vertices now coincides with the state after some teleportation of the input state 
into the vertex b (up to LC). 
(■^^) In n input qubits. 

^26 -J rpj^g complexity of the quantum processing is given by the number of local measurements. For the 
classical processing the number of elementary classical steps for the update of the information flow vector 
and for the determination of subsequent measurement directions is logarithmic in n and proportional to the 
number of measurement rounds. 
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resources. But note that the derived error threshold seems rather unrealistic for practical pur- 
poses. Current research focuses on an improvement of this threshold and a generalization of 
the underlying noise model combining standard concepts of quantum error correction with some 
intrinsic error detection possibilities of the QCq and purification methods for the underlying 
graph state as the computational resource (see sec. l4'4t . Further results on fault tolerant quantum 
computation can also be found in ref . Issll . 

43. Quantum error correcting codes. - For quantum error correcting codes based on stabi- 
lizer codes li^ the codewords as well as the encoding procedures can be represented as graphs 
jiOiEBE^- The latter can be understood along the lines of the previous subsection, because the 
graph state is the computational resource for implementing the encoding process in terms of the 
QCq model. The graphs depicted in fig.^Jfor example correspond to the encoding procedures 
for the five-qubit Steane code and the concatenated [7, 1, 3]-CSS-code that encode a state on one 
qubit into some state on five and 49 qubits respectively. 

The encoding consists in preparing the qubit 
at the (red) input vertex in the input state and the 
remaining qubits in the |+) state, entangling the 
qubits according to the graph and finally measuring 
the (blue) auxiliary vertices (^^) and the input ver- 
tex in the (Ta; -direction. In this way the input state 
is encoded into the state on the remaining (black) 
qubits. Alternatively one might as well prepare the 
graph states depicted in the figures as such, tele- 
port the input state into the input vertex and then 
perform the same measurement protocol. 

When applying this procedure to an eigenstate 
of one of the Pauli-matrices cTx, Uy or a^, the encoded state is a graph state and can be regarded 
as a code word vector for the respective quantum code (see sec. l2'2.lV Note that the [7, 1, 3]- 
CSS-code is a concatenation of the Steane code depicted in fig. El This means that each of the 7 
output qubits of the first encoding level is encoded again by the same encoding procedure, such 
that the overall encoding procedure maps some input state into an encoded state on the 7 x 7 = 49 
output qubits of the second level. This is graphically reflected by the fact that the right graph in 
fig.^Jis obtained from the left graph in fig.Elby attaching another cube at each 'output' -corner 
in the initial cube and merging the former output and input vertex into one auxiliary vertex. 

4'4. CSS-states and secret sharing. - The class of CSS states corresponds to the class of two- 
colorable graphs, i.e., graphs that allow for a coloring of the vertices with two colors such that 
no two adjacent vertices have the same color (see sec. ll 2t . In general a CSS-code is a stabilizer 
code, whose stabilizer can be generated by stabilizer elements cr that either consists of (Tx- or 



(^'^) Whereas the five-qubit code does not require any auxiliary qubits to be measured, the auxiliary vertices 
of concatenated CSS-code are located at the corners of the inner cube, where the open side indicates the 
input qubit. 




Fig. 11. — Five-Qubit-Code and concatenated CSS-Code- 
The graphs representing the encoding procedure for the 
five-qubit and the concatenated [7, 1, 3]-CSS-code with in- 
put (red), auxiliary (blue) and output (black) vertices. 
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CTz -matrices, i.e., cr 



By performing Hadamard operation to the vertices of one 
coloring partition any two-colorable graph state is easily seen to be of CSS-type(^^). Conversely 
it is shown in ref. 14711 that any CSS-state is also LC-equivalent to some two-colorable graph 
state. 



In (431 CSS states are used to design a 'prepare 
and measure' protocol for quantum cryptography 
that can be used for conference key agreementi^^) 
and quantum sharing of classical secrets{^^). 

We remark that not all Pauli measurements pre- 
serve the two-colorability of the underlying graph. 
Whereas local a^- or cr^ -measurements in two- 
colorable graphs yield graph states according to 
two-colorable graphs, try -measurements of two- 
colorable graphs can lead to graph states which are 
not even locally equivalent to two-colorable graphs IstIi . This might be important in the con- 
text of 1-way QC models, in particular in the context of fault tolerance. There, graph states 
corresponding to the non-clifford part of a computation may be purified via multi -particle entan- 
glement purification, where current protocols require two-colorability of the underlying graph. 




Fig. 12. - The encoding graph for the Steane 
code is two-colorable. 



4'5. Entanglement purification and secure state distribution. - Despite of their central use 
for quantum information processing purposes, in reality, multi -partite entangled pure states, such 
as graph states, will not be available with unit fidelity. The reasons for this are manifold. For 
instance, the operations to create these states are always noisy, the distribution of these states 
among distant parties occurs through noisy quantum channels, and the storage of states is sub- 
jected to decoherence. For the class of two-colorable graph states entanglement purification 
procedures are known to maintain or purify these states. Some of these procedures 13811 provably 
apply to the regime of noisy control operations. The basic idea in entanglement purification is to 
use several copies of low-fidelity entangled states in order to distill a few highly entangled states. 
Thereby, one copy of the state with low fidelity is measured in order to reveal information about 
either one copy (recurrence and pumping schemes) or several other copies (hashing protocols) of 
the imperfect graph state. For the purification step each of the imperfect copies (mixed states) is 
transformed('^^) into a mixed state that is diagonal in a graph state basis associated with the two- 
colorable graph state in question. Given that the initial fidelity is sufficiently high, an ensemble 
of these graph diagonal states can then be purified. In Ij^ genuine multi-party entanglement pu- 
rification protocols based on a recurrence or pumping schemes and based on hashing procedures 



{^^) I.e., a CSS code with stabilizer of full rank rank(5) = iV. 

(^^ ) This protocol allows several parties to share a secure conference key in the presence of an eavesdropper. 
|-30-j -pj^jg protocol allows a party to share a classical secret with n parties, such that at least k other parties 
are needed to reveal this key (for some fixed fc > ^). 

('^^) This purification is done by means of some local twirling operations (see Proposition l23l in sec. llOk so 
with respect to the average with respect to an appropriate local symmetry group. 
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(see also 14711 ') were considered and compared with procedures using bi-partite purification (■^^). 
The different purification schemes have their own advantages. But for local noise, multi-party 
purification protocols are, in most cases, not only more efficient but provide also better achievable 
fidelities for the distilled states in the case of noisy control operations. For practical purposes it is 
useful that the purification regime('^'^) for the recurrence protocols as well as the error threshold 
for the noisy control operations involved seem not to depend on the number of particles in the 
graph state but rather on the maximal degree of the underlying graph. Thus two-colorable graph 
states provide a reservoir of entangled states between a large number of particles, which can be 
created and maintained even in the presence of decoherence. Entanglement purification schemes 
will be discussed in more detail in sec. 

These purification schemes were also modified to blind purification protocols. That is, the 
purification takes place in such a way that even the involved parties (except of a central party) 
do neither know the identity of the state they are attempting to purify nor have a possibility to 
learn it. Such protocols allow for the secure distribution of two-colorable graph states, which 
remain unknown to the different parties and any potential eavesdropper |4^ . Such secure state 
distribution is a quantum primitive in a real-world scenario (i.e., taking imperfections in local 
control operations and channel noise into account) which may be used as a basic building block 
for both quantum and classical security applications. 

5. - Physical implementations 

Let us briefly discuss a few proposals of how to prepare graph states in real physical systems, 
aiming at realizing some of the applications, such as the one-way quantum computer (QCq). In 
general, graph states do not appear as ground states of physical systems, since, in most cases, they 
rely dominantly on two-body interactions Issll . Following the preparation procedure in sec. 12' II 
graph states can nevertheless be obtained in any physical system that allows to implement an 
Ising interaction = tr^cr^. In particular, graph states can be generated, of course, in all 
physical devices for universal quantum computation. Clearly, the complexity of the graph state 
preparation depends on the respective computational primitives of the physical realization js^. 

One may distinguish two classes of physical systems where cluster or graph state prepara- 
tion very naturally reflects the underlying architecture. The first one consists in lattice systems, 
such as cold atoms in optical lattices. Here, a neighborhood relation is inherited by the physical 
neighborhood of constituents in the lattice itself, and by means of appropriate switching of in- 
teractions, one can generate cluster states of some dimension IstIi . The second one comprises of 
physical systems where the neighborhood relation can be freely chosen, such as in purely opti- 
cal systems or hybrid setups making use of 'flying' optical qubits for the entangling process and 
'static' matter qubits for storage of the actual graph state. We will briefly sketch both approaches, 
starting with the former. 



(^^) Here one e.g. establishes some highly entangled pairs shared between different parties, in order to 
teleport the desired -particle state prepared at one party to the remaining parties. 

(^^) I.e., the purification regime is given by the necessary fidelity of the initial low-entangled copies, such 
that some higher-entangled state can be distilled. 
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The preparation of a two- or three-dimensional cluster state (of 'arbitrary' size) can be re- 
garded as a fixed operational step, if the Ising interaction can be implemented homogeneously 
(for all directions) throughout the whole lattice Isil IstIi . This globally tunable Ising interac- 
tion can be realized in optical lattices with ultra-cold atoms via state-selective displacement of 
the atoms and controlled cold collisions t88lil or via tunneling looll . For instance, employing a 
superfluid-Mott insulation quantum phase transition a Bose-Einstein condensate can be loaded 
into an optical lattice achieving almost unit occupation per lattice site Is^ I91I1 . Interference 
experiments i92l indicate that graph states can be obtained within such setups. Graph states 
prepared in optical lattices therefore are a promising resource for scalable quantum computation 
in the framework of the one-way quantum computer, since cluster states can be created with an 
operational effort that is independent of the system size. Existing implementations are currently 
still facing at least two major challenges; On one hand, the number of atoms at each lattice site 
has to be exactly controllable. In particular, the presence of defects, i.e., empty lattice sites, might 
spoil any quantum computation. On the other hand, the one-way computer requires the atoms 
at different sites to be individually addressable by local measurements. A number of proposals 
have already been made to overcome these obstacles (see, e.g., refs. and one can 

say that optical lattices remain one of the prime candidates to study multi-particle entanglement 
and scalable quantum computation. Alternative implementations have also been suggested for 
physical systerns, in which the underlying Hamiltonian is given by an (isotropic) Heisenberg 
interaction 196119711 . These proposals can be realized in XY-spin chains formed within an optical 
lattice of neutral atoms as well as in solid-state systems, such as quantum dots. 

The second type of physical systems does not directly exploit immediate adjacency of the 
respective constituents: In this class of physical systems the key point is the fact that one in- 
troduces a separation between the act of creating entanglement and the act of performing the 
actual computation. This is of central importance in particular in setups where the elementary 
gates function on a probabilistic basis, or where erasure errors such as photon losses constitute 
the predominant obstacle that has to be overcome. This very much applies to architectures of 
quantum computing based on linear optics. Linear optical setups are attractive, as photons are 
relatively robust with respect to decoherence, and accurate state control is possible using linear 
optical elements L98.1 . However, in order to achieve universal quantum computing based on dual 
rail encoding - where logical qubits are encoded in state vectors 1 0) 1 1) and 1 1) 1 0) of two modes - 
measurements are necessary, rendering any scheme probabilistic. Notably, the scheme of Knill, 
Laflamme, and Milburn in ref. I98il employs the non-linear sign shift gate, acting as 

(107) xo|0)+a-i|l)+X2|2) K^a;o|0)+.Ti|l)-a-2|2) 

as a primitive probabilistic gate, realized using additional modes and measurement. Here, |0), 
and 1 2) denote the state vectors of states containing 0, 1, 2 photons. It has been shown that, 
although these gates operate on a probabilistic basis (only certain measurement outcomes are 
accepted as being successful), the overall computational scheme can be uplifted to almost unit 
success probability. This is done using teleportation, based on appropriate entangled resources 
which are in turn built up by invoking these probabilistic gates. This is a powerful theoretical 
idea, heavily exploiting the previously described stabilizer formalism. So in this way, near- 
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deterministic quantum computation is possible in principle using linear optical elements, single- 
photon sources, photon-resolving detectors, and feedforward, meaning later action may depend 
on earlier measurement outcomes Ig^l . Practically, the key obstacle is, however, that the required 
resources are tremendous (approximately 10^ beam splitters required for a single CNOT gate 
operating with a failure rate of 10^''), essentially originating from the fact that the elementary 
gates work with such a small probability of success 19911 . 

Graph state methods can indeed significantly reduce the surmounting number of required re- 
sources (although still giving rise to challenging prescriptions). In the remainder of this section, 
we will thus briefly review recent work discussing grap h states and one-way quantum compu- 
tation in the context of linear optical setups. Ref. llOOll is the firs t wor k making use of states 
related to graph states to reduce the overhead in resources. In ref. jlOlll . the presented scheme 
has explicitly been phrased in terms of cluster and graph states, making use of the setting of the 
scheme by Knill-Laflamme-Milburn. 

The number or necessary resources was further reduced in refs. jl02L 103 1. employing two 



types of fusing gates that glue pieces of linear one-dimensional clusters or of two-dimensional 
graph states together Essent ially, the basic resource here are Bell states using non-deterministic 
parity-check measurements 110411 . involving the combination of photons on polarizing beam 
splitters, followed by measurements on the output modes. The computational basis is again 
the one of dual rail encoding, here specifically the basis 

(108) |i^):=|0)|l), |T^):-|1)|0) 

corresponding to a horizontally and vertically polarized photon. A fusion of the first type 
amounts to mixing the input, assumed to contain only one photon per spatial mode, at a polariz- 
ing beam splitter, rotating the output by 45° and measuring it with a polarization-discriminating 
photon detector This is the parity-check operation considered in ref. 1,104]. In the case when one 
and only one photon is detected - occurring with a probability of 1 /2 and considered the suc- 
cessful event - the state is transformed according to a non-trace-preserving completely positive 
maps with either of the two Kraus operators, 

(109) K, = mmm - \v){v\{v\)/V2, 

(110) K2 = i\H){H\{H\ + \V){V\{V\)/V2. 



Given such a successful event, one glues two pieces of a graph state together 110211 . From max- 
imally entangled two-dimensional res ourc es, one can build up linear cluster states. For a first 
experimental demonstration, see ref. jlOSll . To b uild up h igher -dimensional structures, a vari- 



ant of a destructive CNOT gate can be used llOJ], in ref. I102II referred to as fusion of second 



type. This fusion gate of the second type consists of a polarizing beam splitters, two polar- 
ization rotators of 45° at the inputs and outputs each, followed again by polarization resolving 
photon detection. These gates are still probabilistic, yet, failures can be tolerated by merely re- 
building the affected graph section, assuming that failures are heralded, meaning that one has 
classical knowledge of inappropriate measurem ent o utcomes. First steps towards fault tolerance 
with respect to photon losses have been taken In this way, graph states can be prepared 
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from essentially probabilistic ingredients, giving rise to deterministic quantum computation with 
smaller overhead in resources. 

Quite recently, arbitrary graph states on four qubits were experimentally generated in an 
entirely optic al system, with the four qubits being represented by the polarization state of four 
photons llOTII . Although the feasibility of 1-QC was demonstrated through a set of one- and two- 
qubit operations, at the present stage such experiments rather constitute a proof-of-principle than 
a clear-cut route to scalable quantum computation. This is mainly because the scalability of the 
preparation by means of parametric down conversion is currently bound to only a few photons. 
Using such technology as source of entangled photons, hence, the overall success probability is 
exponentially decreasing. 

An alternative route in the second framework is provided by hybrid solutions, where the ad- 
vantages of photons providing coupling mechanisms and of matter systems as long-lived storage 
devices are simultaneously exploited. Based on schemes that allow for the prepa ration of e ntan- 
gled states using matter systems in leaking cavities using flying optical q ubits lllOSi llOOll . one 
can immediately construct schemes that allow for graph state preparation 1 1 lol iTll 1 12il . even 



in setups where in intermediate steps no further local rotations of the state of the matter qubits is 
required. The matter qubits may be spatially separated, and the adjacency pattern of the graph 
states to be prepared is in then, in principle, completely arbitrary. Such schemes can even be 
made essentially deterministic, in that if an entangling operation fails, it can o ften be rep eated 



until a successful operation occurs, without damage to the nascent graph state 111 1 IL ll 1211 . Us 



ing multi-port linear optical devices, the need of doing intermediate local rotations can also be 
largely eliminated in the preparation of the graph states. Also, whole parts of graph states can 
be fused together, in an essentially deterministic manner. In architectures making use of opti- 
cally active electron-nuclear systems, such as in N-V-centers in diamond, one has effectively a 
(D^ (g) (D^-system locally available, allowing for a wa y to e ntangle matter qubits in a way that is 
more robust with respect to unavoidable photon loss 11 1311 . All these schemes have in common 
that the neighborhood relation between constituents is fairly arbitrary, and that arbitrary graph 
states can - in principle - directly be prepared. 

One general lesson to learn from graph state approaches from the perspective of implemen- 
tations is that there are probably no generally valid uncompromising requirements for physical 
systems to allow for scalable quantum computation. Instead, it very much depends on the un- 
derlying physical architecture what computational model is advantageous; hence rendering oth- 
erwise inappropriate architectures useful and promising. Schemes for quantum computing based 
on graph states promise to significantly lessen the challenges posed by the obvious requirement 
in any scheme for quantum computation, that one needs extraordinarily good isolation against 
decoherence effects and precise state manipulation at the same time. 



6. - Reduced states of graph states 

As we will discuss in sec. |5| for pure multi-partite states the reduced state obtained after 
tracing out (forgetting the information in) some part of the system captures many interesting 
entanglement properties of the state. More precisely, consider a pure state \'4))^'^ of a joint 
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system of parties A and B. On the one hand the reduced state 
(111) p^:=trB{p) 

represents the information available for one subsystem A, if it is not provided with any infor- 
mation that corresponds to some measurement statistics on party B. The study of entanglement 
accessible to subsystem A in this way is an interesting issue in itself. On the other hand, for pure 
states p = IV') (V"! the state is entangled with respect to the partitioning {A, B) iff the reduced 
state p^ = tr^(p) is mixed and thus cannot be written as a product state = . 

Moreover, the mixedness of this reduced state |V')'^^- for example in terms of some entropic 
measure, allows one to quantify the amount of entanglement contained in this state between the 
parties A and B. 

The following proposition shows that for graph states the reduced density matrices can be 
represented efficiently in terms of their stabilizer elements or their adjacency matrix 

mini 



Proposition 8 (Reduced state). Let A Q V be subset of vertices for a graph G = (V, E) and 
B ^ V \ A the corresponding complement in V. The reduced state p^ :~ trs (|G) (Gj) is given 
by 

where 

(113) Sa ■■= W e S\supp{a) C A} 

denotes the subgroup of stabilizer elements a € S for \G) with support('^^) on the set of vertices 
within A. Pq is up to some factor a projection, i.e., 



(114) {Pt) =^Pg- 

// projects onto the subspace in H'^ spanned by the vectors 



(115) |r'S')G[^, = ap' |G[^]) (B'CB), 

where G[A\ = G\B is the subgraph ofG induced by A and T' := denotes the \ A\ x 
off-diagonal sub-matrix of the adjacency matrix T for G that represents the edges between A 
and B: 



(116) 



Tab 



('^*) The support of a Pauli operator a = a]^ ® . . . ® crf^ is the set of all indices a G for which a acts 
non-trivially, i.e., ia 7^ 0. 
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In this basis, can be written as 

(117) '°G = i E rB')G[A]{rB'\. 

B'CB 

Proof: Eq. (II 12> immediately follows from eq. ( I30> and the fact that the partial trace of 
(7 = cr^ J <X) • ■ • (8) cr/^ € S can be taken successively over the different vertices in e B and gives 
trb((Tj^J = 2(5i^o- Pg is proportional to a projection, i.e., 

follows, because Sa is a subgroup of S. To show eq. ( II 17> . the partial trace over B can be taken 
in the basis of B given by 

(119) \B')f^^\B',)l:^a^'\0)^. 

bGB 

This basis decomposition corresponds to successive local ct^ -measurements of all vertices in B. 
The set B' C B or the corresponding binary vector determines the measurement outcomes, 
i.e., if b ^ B' or likewise the corresponding component of the binary vector = then the 
measurement outcome at this vertex is +1 and —1 otherwise. According to Proposition^ after 
measurement of cr^ the state of the remaining vertices is the graph state vector jG \ 6) in the case 
of the outcome +1 , and cr^''\G \ h) if the outcome is —1. This can be summarized to 

(120) {a^^f'\G\h), 

since G {0,1} represents the measurement result { + 1,-1}. Because the subsequent mea- 
surements commute with the previous local unitaries, the final state vector according to the result 

B' = {B'^)beB e Ff is 
(121) 

\{{a^^f^ \G\B)Mb')'^ = al^-«^'^^)|G[A])^®|B')f = af^' |G[A])^®|i?')f , 

because the sum is performed modulo 2, and for all e i? the binary vector corresponding to 
Nb is the 6-th column of the adjacency matrix T. Distinguishing the two parts of the adjacency 
matrix T b and Tab that correspond to edges within B, and edges between A and B, one arrives 
at the state vector associated with the measurement result given by B' 

(122) {-l)^B'\TsB') ^T_^sB' ^ ^ (_^yB'|r,i3') \Y abB%^j,^ ® |B')f ■ 

Because the possible measurement results are attained with probabihty 1/2, this proves eq. ( I117l i. 
□ 
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7. - Equivalence classes under local operations 

For a characterization of the entanglement in graph states we now examine the equivalence 
classes under local operations. Concerning locality we restrict to the finest partitioning, i.e., each 
vertex a G V represents a single party and the quantum operations £ in question are a certain 
subclass of the class of all completely positive maps (CPM) that are separable with respect to 
this finest partitioning. Since we are interested in the equivalence classes of graph states under 
these local operations, we can consider the situation in which a pure state jV'i) is mapped onto 
another pure state {^2) by the CPM £ with non-zero probability. It is generally quite a subtle 
problem to characterize the class of all transformations £ that can be implemented by means 
of local operations and classical communication (LOCC), and in the following we restrict to a 
subclass of LOCC -protocols where £ factors out as the tensor product of a local operator Ei for 
each party: 

(123) £{p) ^El®...®E%p {E\)^ ®...® {E^y . 

This means that the pure state is converted (in general: stochastically) into the state \^'2) = 
(8) ... (8) E^ In the following we will consider three different classes of local operations, 
namely 

• SLOCC: invertible stochastic local operations and classical communication, i.e., the op- 
eration Ea E SL(2, C) at each qubit is an arbitrary invertible matrix; SLOCC-equivalence 
occurs typically with nonunit probability; 

• LU: arbitrary local unitaries, i.e. the operation Ea G SU(2) at each qubit is some unitary 
operation; LU-equivalence occurs with unit probability; 

• LC: local Clifford unitaries, i.e., the operation Ea G Ci at each qubit is one of the Clifford 
unitaries introduced in sec. l2'2.2l LC-equivalence occurs with unit probability; 

For X S {LC, LU, SLOCC}, two states and \tp2) are called X— equivalent if there exists an 
X— operator X such that X|-02) ^ \tpi)- We will frequently write 

(124) IV'i) IV'2) 

if both states are X-equivalent. 

Let us first give a brief overview of entanglement classes for arbitrary states, i.e., not neces- 
sarily graph or stabilizer states. We start with the case of two d-level systems, i.e., we consider 
bi-partite entanglement. Here the Schmidt decomposition serves as a well-known standard form, 
from which the conditions for different types of interconvertibility can be read off. A bi-partite 
state ji/'i) can be transformed into the state |-02) 



(125) 



iV'i) — *x IV'2) for X G {LU, LOCC, SLOCC} 
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by means of (i) SLOCC, (ii) LOCC or (iii) LU operations iff for the corresponding Schmidt 
decomposition of 

(126) \^,)^^ =Y.Xl\^)^\^)^ and |^,)^^ = f] A^b^ |,)^ 

1=1 3=1 

(i) Ri > i?2, (ii) the coefficient fist iX})^i is a majorization of (A^)j^j^ or (iii) the coefficient 
list {Xl)^i coincides with {X'j)f^i up to permutations, respectively. It follows that in the bi- 
partite case the SLOCC-equivalence leads to d distinct equivalence classes corresponding to the 
different values for the Schmidt rank R. On the other hand, there is an infinite number of bi- 
partite LU-equivalence classes, which are parameterized by the d Schmidt coefficients A;. 

Moving to multi-partite entanglement, in the case of three qubits it was shown in ref. jl isil 
that there are 6 orbits under SLOCC-operations. But true three-qubit entanglement is only con- 
tained in the two distinct classes represented by the GHZ-state \GHZ) (see sec. 14' U and the 
W-state 

(127) \W) = (|0,0,1) + |0,1,0) + |1,0,0)) . 

Note that all two-qubit states are SLOCC-equivalent to a two-qubit graph state corresponding 
to either the empty graph (product state) or connected graph (Bell state). But the W-state is an 
example for a pure state on three qubits that is not SLOCC-equivalent to some graph state. This 
is because the GHZ-state \GHZ) is the only 'connected' graph state with three vertices, as we 
will see below. 

For A'^-qubit systems with > 4 the number of orbits under SLOCC-operations is infinite 
and is in fact specified by an exponentially increasing number of parameters. The latter scaling 
behavior for the number of required parameters is due to the fact that a generic orbit is speci- 
fied (up to some irrelevant complex constant) by 3iV complex parameters describing the group 
SL(2,C)'^^ in question whereas a pure state on (C^)®^ is specified by 2^ — 1 co mple x pa- 
rameters (neglecting an overall complex phase). For 4 qubits it was shown in ref. 11 1611 that 
there exists a standard form for generic states under SLOCC-operations that is determined by 3 
complex parameters and 8 further standard forms corresponding to classes of 'degenerate' pure 
states. Needless to say that the number of LU- or LC-equivalence classes necessarily has to 
be specified by even more parameters since the corresponding matrix groups are contained in 
SL(2,C). 

In 11 1711 it is shown that any multi-partite state can be transformed by SLOCC-operations into 
some standard form that is unique up to LU, while maximizing all entanglement monotones('^^). 
In this way the problem of deciding, whether two states are SLOCC-equivalent, can be reduced 
to the problem of deciding, whether the standard forms of these states are LU-equivalent. In this 



('^^) More precisely the entanglement monotones in question are linearly homogeneous positive functions 
of the state that remain invariant under determinant-! -SLOCC operations. 



50 



M. HEIN, W. DUR, J. ElSERT, R. RAUSSENDORF, M. VAN DEN NEST and H.-J. BRIEOEL 



formalism all pure states with maximally mixed single-qubit reduced density matrices 

(128) p-^try\^(\^b){ij\)^^la yaeV 

are standard forms of this local filtering sequence and they are maximally entangled in that these 
states maximize all entanglement monotones. 

Let us return to the question of equivalence under local operations for graph states. We start 
with SLOCC-equivalence. According to Proposition |S] any single-qubit reduced state is 
maximally mixed for all graph states corresponding to connected graphs. Therefore, 'connected' 
graph states are already in standard form under SLOCC in the above sense. Thus, when re- 
stricting to connected graph states all SLOCC-equivalence classes coincide with LU-equivalence 
classes; this property can easily be extended to the case of general (not necessarily 'connected') 
graph states by considering each connected component separately. We arrive at the following 
result: 

Proposition 9 fSLOCC- equals LU-e quiv alence). Two graph states \ Gi) and \Gi) are SLOCC- 
equivalent iff they are LU- equivalent Ul4l : 

(129) |Gi) ^sLocc IG2) ^ |Gi) IG2) . 



This result was first obtained in ref. II 14ll . We remark that the question whether for graph 
states also 

(130) |Gi)^Lu|G2) =^ |Gi)^Lc|G2) 

holds is still an open question. Note that the backward implication is trivial since the group of 
LC-unitaries i s a p roper subgroup of all LU. For a large subset('^^) of graph states, however, it 
was shown in 11 1811 that both notions of equivalence coincide. The hypothesis of a general coinci- 
dence for all graph states is further supported by results H 1911 about the corresponding invariants 
under these operations, which will be briefly reviewed later in this section, and the classification 
of graph states with up to seven vertices that we will discuss below. Note that a general coinci- 
dence of SLOCC-, LU- and LC-equivalence for graph states would be particularly advantageous 
for the following two reasons: since in this case all 3 local equivalences would correspond to 
LC-equivalence, firstly, checking whether two given graph states are locally equivalent could 
then be done efficiently (see sec. l2'2.4l : secondly, all 3 locally equivalences would entirely be 
described by the local complementation rule, yielding a description of local equivalence of graph 
states in simple, purely graph theoretic terms. 



{'^^) A stabilizer element a G 5 with minimal support supp((T), i.e., no other stabilizer element cr' G 5 has 
a support supp(cr') that is a proper subset of supp(cr), is called minimal element of S. Let Ai denote the 
subgroup generated by all minimal elements in S. Now the notion of LU- and LC-equivalence coincide for 
all stabilizer states, for which a^, ay and cr^ occurs at each qubit in Ai. 
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To distinguish the different equivalence classes under LU (or, equivalently, SLOCC) we will 
first derive a simple (though not complete) set of invariants that can be efficiently computed for 
graph states from the underlying adjacency matrix of the graph. As mentioned above, for any 
pure state in the joint system of two parties {A, B) with arbitrary dimension- 

ahty dA = dimcH^i and ds ~ dimcH^ the Schmidt rank('^'') SR"^('0) := rank (^P^^ is an 
entanglement monotone with respect to (^4, i3)-local operations. In the case where {A, B) is a 
bi-partition of a many-qubit system we will for simplicity consider not the rank SR'"^(7/') of the 
reduced states but rather the logarithm of it with respect to the basis 2, i.e.. 



(131) 



SR^(V) :-log2 [rank (4)] 



With this notational simplification the Schmidt rank SRa(^) of pure state ItA) has the transparent 
interpretation('^'*) in terms of the maximal number of Bell pairs that 

• are required to prepare | iji) and 

• can be extracted from 

with finite probability of success (i.e., under SLOCC-operations). 

We will now study the entanglement 
in graph states \G) by considering the 
Schmidt rank SRyi(G) with respect to differ- 
ent bi -partitions (A, B). For a bi-partition 
of a graph G = {V,E) we can 
again use the decomposition for the adjacency 
matrix T of eq. il 16> into Ta, and Fab 
according to edges within A, edges within 
B and those edges between A and B. In 
this notation the corresponding Schmidt rank 
SRa(G) is simply given by the binary rank 
(i.e., the rank over GF(2)) of the \A\ x \B\- 
off-diagonal sub-matrix T' ~ Tab- 



Proposition 10 (Schmidt rank). Let (A, B) 
be a bi-partition for some graph state \G). 
Then the Schmidt rank of the graph state with 
respect to this bi-partition is given by 




(132) 



SRa{G) = rankw,_ {T') 



Fig. 13. - An example of an equivalence class which 
is a proper subset of the class No. 4 in List A but 
which is not LC-equivalent to any of the graphs de- 
picted in fig. m With the LC-rule in sec. I2'2.2l it is 
straightforward to check that all graphs within this 
class are LC-equivalent. 



(^^) Note that from the Schmidt decomposition it follows that SR'^{tP) = SR'^{%p). 

This follows straightforwardly from statement (i) in this section about the interconvertibility of pure 
bi-partite states under SLOCC-operations. 
(^^) I.e., A U B = 1/ and A n B = 0. 
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A 


{1} 


{2} 


{3} 


{4} 


{1,2} 


{1,3} 


{1,4} 


A 


{1,2,3,4} {2,3,4} 


{1,3,4} 


{1,2,4} 


{1,2,3} 


{3,4} 


{2,4} 


{2,3} 



SRa(Gi) 





1 


1 


1 


1 


1 


2 


2 


SRa(Gi) 





1 


1 


1 


1 


2 


1 


2 



Table III. - The list of Schmidt ranks Si?A(Gi) and SRa{G2) for the graphs Gi in fig.^and for the 
graphs G2 in fie.U3\ 

Alternatively, the Schmidt rank is determined in terms of the ratik{'^^) of the subgroup Sa of 
stabilizer elements with support in A by the formula: 

(133) SRa{G) = \ A\ - rank{SA) ■ 

Proof: Because of Proposition|3j the linear independence of the vectors \T' B')QyA] {^^ ) is in 
one-to-one correspondence to the linear independence of the corresponding vectors T' B' over 
F^. Hence we find for the rank of the reduced state 

rank (p^) = dime span { \T'B')g[a\ \B' ^B] 
= dimp^ span {V B' | B' C 5} 

(134) = ranky^ (F') . 

Taking the log2 we obtain the eq. J132> . Eq. jl33t follows from the fact that according to 
eq. M 14t tq-^Pg is a projection and thus we can compute the rank alternatively as 



(135) 




Taking again the log2 we obtain eq. since the minimal number of generators for Sa is 

givenbylogad^Al). □ 

Thus, for any partition (A, B) the Schmidt rank SRa(G) is an invariant under arbitrary local 
unitaries that can be formulated in purely graph theoretic terms. We now consider the list of 
Schmidt ranks with respect to all possible bi-partitionsC'^). This yields a set of invariants which 
has already been considered in graph theory under the name connectivity function js^. For 
example, a comparison of the invariants for the graphs depicted in fig.0]and fig.ll3lshows that 



("O) See eq. 01} in sec. lITll 

(*^) In the graph state basis for the graph G[A] of Proposition[8| 

(*^) Since SRa(G) = SR^(G) the different bi-partitions are fixed by choosing the smaller partition, say 
A, of the bi-partition [A, B). This gives 2^~^ possible bi-partitions. 
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6 




10 



Fig. 14. - The Petersen Graph. The depicted labeled graph is not LC-equivalent to the graph which is 
obtained from it by exchanging the labels at each end of the five "spokes", i.e., the graph isomorphism which 
permutes the vertices 1, 2, 3, 4 and 5 with 6, 7, 8, 9 and 10, respectively. However, the lists of Schmidt ranks 
(or, equivalently, the connectivity functions) of these graphs coincide. 



the corresponding lists of Schmidt ranks within each of these figures coincide but differ between 
the two figures (see Tab. Illl> . This impfies that the corresponding sets of graph states are not 
equivalent neither under LC-operations nor under general local unitaries. 

We note that the Schmidt rank list does not provide a complete set of invariants that would 
characterize all equivalence classes under LC-operations. For the Petersen graph{'^^) shown in 
fig.^Jand the isomorphic graph, which is obtained from it by exchanging the labels at each end 
of the five "spokes", no local complementations exists that transforms one graph into the other, 
although the Schmidt rank lists for both graphs coincide. 

In illi a complete set of polynomial LU invariants for graph states is obtained (in terms 
of binary trees). As stated in eq. ( I133t the set of invariants {SRA(G)}yicy, corresponding to 
polynomial invariants of degree k — 2, can for stabilizer st ates b e formulated in terms of the 
dimension of a subspace of the stabilizer S. It was shown in jl 19ll that, similarly, the invariants 
of degree k > 3 cor respond to dimensions of certain subspaces of the k — 1-fold direct product 
^x(fc-i) j.gj^ jl 19I1 also provides a finite complete set of polynomial invariants for the smaller 
group G = Ci of local Clifford unitaries (see sec. l2'2.2l . For graph states these LC-invariants 
of degree k are again given by the dimension of a subspaces in 5^'-'''^^-'. Note that the algebra 
of polyn omial LC-invariants is in general larger than the algebra of polynomial LU-invariants. 
In 11 1911 it is shown that the set of polynomial LU-invariants for the degree k = 2 and k ~ 3 
are equivalent to the corresponding sets of LC-invariants. These results support the conjecture 
(see sec. l2'2.2l and sec.O that for graph states the notion of LU-equivalence and LC-equivalence 
coincide. 



(*^) This counter-example has first been discovered in ref. ll2Cll . 
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In the remainder of this section we will now discuss LU-equivalence for graphs with N < 7 
vertices C^). We have examined the graph states of all non-isomorphic, connected graphs 
with up to seven vertices. More precisely, from the set of all possible graphs with 7 vertices 

(2(2) « 2 X 10^ possibilities), we consider the subset of 1252 graphs on up to 7 vertices that 
are non-isomorphic with respect to graph isomorphisms. Note that a graph isomorphism phys- 
ically comes up to an exchange of particles. We furthermore restrict to those 996 states that 
correspond to connected graphs. This is because a state \G) corresponding to a disconnected 
graph G is simply the tensor product \G) = \Gi)^^ (g) . . . (8) IGa/)"^" of the graph states \Gi) 
corresponding to the connected components Gi of G, where {Ai, . . . , Am) {M < N) is some 
partitioning of the vertex set V. Thus all entanglement properties of the composite state | G) are 
essentially determined by the entanglement properties of its components | G, ) . In particular, such 
a {Ai , . . . , ^A/)-product state is LU- or LC-equivalent to some other graph state |G') iff |G') al- 
lows for a decomposition |G') = \G'i)"^^ ig) . . . (g) \G'm)^^' with respect to the same partitioning 
{Ai, . . . , Am) with components \G^) that are LU- or LC-equivalent to the respective states \Gi) 
for |G). Of the 996 isomorphism-classes of corresponding graph states, 46 classes have turned 
out to be not invariant under local unitary operations. 

Within each of these classes all graph states are equivalent modulo local unitaries and ad- 
ditional graph isomorphisms. Thus, if we exclude the graph isomorphisms, as e.g. in quantum 
communication scenarios, the number of inequivalent classes of graph states is even larger (see 
Tab.lTvl. 

In List A and B of Table |V] we give a list of simple representatives of each equivalence 
class together with a table summarizing some interesting properties of these states. For this we 
have generated the of 996 non-isomorphic, connected graphs with the MATHEMATICA package 
described in secQJand tested for local equivalence considering only LC-unitaries (see sec. l2'2.4l . 
By considering the Schmidt rank with respect to all possible bi-partitions, the corresponding lists 
of Schmidt ranks for each representative turned out to be different even if we allow arbitrary 
permutations of the vertices. This shows that the found sets of locally invariant graph states are 
minimal even with respect to the larger group of all LU. 

We have also listed the sizes of the corresponding equivalence classes under LU and graph 
isomorphisms, as well as whether two-colorable representatives exist. By the rank index given 
in List A and B of Table |VJ we simply compressed the information contained in the Schmidt 
rank list with respect to all bi-partite splittings, counting how many times a certain rank occurs 
in splittings with either two or three vertices in the smaller partition. For example, the rank index 
RI3 = (20, 12, 3) of graph number 29 means that the rank 3 occurs 20 times in all possible 
3-4-splits, the rank 2 twelve times and the rank 1 only three times. Similarly, because of RI2 = 



^44 -J -pj^g classification of non-equivalent graph states naturally generalizes to the case of stabilizer states or 
codes, for which the graph represents a particular standard form (see sec. irO and sec. l4'3L For the quest 
of good error correcting (self-dual) codes such a cla ssification has recently attracted some attention. A 
similar classificatio n can also be found in ref. I6clll2lll and has been extended to graphs with up to = 12 
vertices in refs. 
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non-isomorphic 




non-isomorphic 


N 


non-isomorphic 


and non-LC-equivalent 


non-isomorphic 


and non-LC-equivalent 




graphs 


graphs 


connected graphs 


connected graphs 



1 


1 


1 


1 


1 


2 


2 


2 


1 


1 


3 


4 


3 


2 


1 


4 


11 


6 


6 


2 


5 


34 


11 


21 


4 


6 


156 


26 


112 


11 


7 


1,044 


59 


853 


26 


8 


12,346 


182 


11,117 


101 


9 


274,668 


675 


261,080 


440 


10 


12,005,168 


3,990 


11,716,571 


3,132 


11 


1,018,997,864 


45,144 


1,006,700,565 


40,457 


12 


165,091,172,592 


1,323,363 


164,059,830,476 


1,274,068 



Table IV. - In the first column the number of graph states with N = 1, ... ,12 vertices is listed that are 
non-isomorphic under graph isomorphisms. The third column instead contains the corresponding number 
of no n-isom orphic graph states that correspond to connected graphs. The values for both columns are taken 
from where they can he found under the sequence number A000088 and A001349. Similarly column 

No. 2 and No. 4 contain the corresponding numbers of graph states that are not equivalent und er graph 
isomorphism and LC-operations. The values in the second column were computed in ref. \7(il together 
with a database of representatives for each equivalence class. By checking the list of Schmidt ranks we 
have shown that the values in the second and fourth column for N < 7 vertices coincide with the numbers 
of non-isomorphic graph states when considering the larger group of LU- or SLOCC-operations. The 
values for N — 10, 11, 12 in both columns were again taken from ref. It^/ . where also a database of 
representatives for each equivalence class can be found. 



(18, 3) the rank 2 (1) occurs 18 (3) times in all 2-5-splits of the graph number 29. As it can be 
seen from Tab.lIIII although the classes of graph states in fig.|4]and fig.[0]have different Schmidt 
rank lists and thus are non-LU-equivalent, both classes have the same rank index RI2 = (2, 1), 
since the rank indices are invariant under arbitrary permutations of the vertices. Thus no graph in 
fig- Sis locally equivalent to any graph in the equivalence class represented in fig. El But both 
belong to the same equivalence class represented by graph No. 4 in List A when considering both 
local unitary transformations and graph isomorphisms. In fact a permutation of vertices 2 and 3 
maps graph No. 4 in fig.[0]onto graph No. 1 in fig.E] 

In Tab. IIVI we have summarized the number of equivalence classes for connected and all 
(i.e., possibly disconnected) graphs and compared them with the corresponding values when 
disregarding LC-equivalence. From a quantum information point of view, which only consid- 
ers states up to LU-equivalence, the class of non-equivalent graph states provides a significant 
reduction of the set of all graph states. The table nevertheless shows that the obtained set of 
non-LU-equivalent graph states is still sufficiently rich to form an interesting subclass of states 
to serve as a starting point for the study of multi-party entanglement. Besides the question of 
local equivalence we will see in the remainder of this thesis that many other interesting entangle- 
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No. 3 No. 4 




No. 18 No. 19 




No. 21 No. 22 




1 No. 


LUclass 


\y\ 






PP 


Rh 


Rh 


2 -col 1 


1 


1 


2 


1 


1 


1 






yes 


2 


2 


3 


2 


1 


1 






yes 


3 


2 


4 


3 


1 


1 




(0,3) 


yes 


4 


4 


4 


3 


2 


2 




(2,1) 


yes 


5 


2 


5 


4 


1 


1 




(0,10) 


yes 


6 


6 


5 


4 


2 


2 




(6,4) 


yes 


7 


10 


5 


4 


2 


2 




(8,2) 


yes 


8 


3 


5 


5 


2 


3 




(10,0) 


no 


9 


2 


6 


5 


1 


1 


(0,0,10) 


(0,15) 


yes 


10 


6 


6 


5 


2 


2 


(0,6,4) 


(8,7) 


yes 


11 


4 


6 


5 


2 


2 


(0,9,1) 


(8,7) 


yes 


12 


16 


6 


5 


2 


2 


(0,9,1) 


(11,4) 


yes 


13 


10 


6 


5 


3 


3 


(4,4,2) 


(12,3) 


yes 


14 


25 


6 


5 


3 


3 


(4,5,1) 


(13,2) 


yes 


15 


5 


6 


6 


2 


2 


(0,10,0) 


(12,3) 


yes 


16 


5 


6 


6 


3 


3 


(4,6,0) 


(12,3) 


yes 


17 


21 


6 


6 


3 


3 


(4,6,0) 


(14,1) 


yes 


18 


16 


6 


6 


3 


3 


(6,4,0) 


(15,0) 


yes 


19 


2 


6 


9 


3 


4 


(10,0,0) 


(15,0) 


no 




No. 


1 LUclass] 






SRlTlilX 


PP 


Rh 


Rh 


2 - col 


20 


2 


7 


6 


1 


1 


(0,0,35) 


(0,21) 


yes 


21 


6 


7 


6 


2 


2 


(0,20,15) 


(10,11) 


yes 


22 


6 


7 


6 


2 


2 


(0,30,5) 


(12,9) 


yes 


23 


16 


7 


6 


2 


2 


(0,30,5) 


(14,7) 


yes 


24 


10 


7 


6 


2 


2 


(0,33,2) 


(15,6) 


yes 


25 


10 


7 


6 


3 


3 


(12,16,7) 


(16,5) 


yes 


26 


16 


7 


6 


3 


3 


(12,20,3) 


(16,5) 


yes 


27 


44 


7 


6 


3 


3 


(12,21,2) 


(17,4) 


yes 


28 


44 


7 


6 


3 


3 


(16,16,3) 


(18,3) 


yes 


29 


14 


7 


6 


3 


3 


(20,12,3) 


(18,3) 


yes 


30 


66 


7 


6 


3 


3 


(20,13,2) 


(19,2) 


yes 


31 


10 


7 


7 


2 


2 


(0,34,1) 


(16,5) 


yes 


32 


10 


7 


7 


3 


3 


(12,22,1) 


(16,5) 


no 


33 


21 


7 


7 


3 


3 


(12,22,1) 


(18,3) 


no 


34 


26 


7 


7 


3 


3 


(16,18,1) 


(18,3) 


yes 


35 


36 


7 


7 


3 


3 


(16,19,0) 


(19,2) 


no 


36 


28 


7 


7 


3 


3 


(20,14,1) 


(18,3) 


no 


37 


72 


7 


7 


3 


3 


(20,15,0) 


(19,2) 


no 


38 


114 


7 


7 


3 


3 


(22,13,0) 


(20,1) 


yes 


39 


56 


7 


7 


3 


4 


(24,10,1) 


(20,1) 


no 


40 


92 


7 


7 


3 


4 


(28,7,0) 


(21,0) 


no 


41 


57 


7 


8 


3 


4 


(26,9,0) 


(20,1) 


no 


42 


33 


7 


8 


3 


4 


(28,7,0) 


(21,0) 


no 


43 


9 


7 


9 


3 


3 


(28,7,0) 


(21,0) 


yes 


44 


46 


7 


9 


3 


4 


(32,3,0) 


(21,0) 


no 


45 


9 


7 


10 


3 


4 


(30,5,0) 


(20,1) 


no 



Table V. - List K: List of connected graphs N = 2, 3, 4, 5, 6 vertices that are not equivalent under LU transformations and 
graph isomorphisms. List B.' List of connected graphs with seven vertices that are not equivalent under LU transformations 
and graph isomorphisms. The corresponding tables list for each equivalence class the number of vertices \V\ and edges 
\E\, the maximal Schmidt rank SRmax, the Pauli persistency PP {see sec. I8'3t , the rank index RI3 and RI2 (for splits 
with 2 or 3 vertices in the smaller partition), the number of non-isomorphic but LU equivalent graphs \LUclass\ and the 
two-colorable property 2 — col. 
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ment properties of graph states have a concise and efficient translation in terms of the underlying 
graph. This allows for an exemplary study of multi-party entanglement in the regime of many 
parties. 



8. - Entanglement in graph states 

As discussed in the previous sections, graph states provide an interesting class of multi- 
partite states that are relatively easy to survey even in the regime of many parties. Since the 
graph essentially encodes a preparation procedure of the state, we will now examine the question 
how the entanglement in a graph state is related to the topology of its underlying graph. More 
precisely we address the issue of quantifying and characterizing the entanglement of graph states. 

We start in sec. 18' II with a review of results from 11501 Il24il about the 'non-classicality' of 
graph states and how the entanglement present in these states can be experimentally verified 
considering Bell inequalities or entanglement witnesses. Then, classical correlations and entan- 
glement between pairs of particles are discussed in sec. I8'2I The main part of this section is 
finally devoted to the quantification of entanglement in graph states in terms of the Schmidt mea- 
sure, which will be introduced in sec. 18 31 We present bounds and rules that render the evaluation 
of this measure feasible for some interesting classes of graph states and discuss some examples 
of practical interest. 

8' 1 . Bell inequalities and entanglement witnesses. - The notion of entanglement in quantum 
mechanics as it was posed in ref. |l|, yj] by Einstein, Podolsky, Rosen and Schrodinger in the 
year 1935 has - since Bell's reformulation 0] in 1964 - frequentl y bee n used synonymous with 
'non-classical correlations', although today it is well-known l23[ll25ll that one has to consider 
a finer distinction. Up to some so-called 'detection loophole', first experiments jsl were able 
to verify that some quantum states can indeed reveal correlations which cannot be predicted by 
some local hidden variable (LHV) models. In these models any observable has a predetermined 
value (realism), regardless of whether it is measured or not. Moreover, the choice of which 
observable is measured does not 'affect the other parties' (locality). The two constraints can be 
phrased in terms of Bell inequalities, which bound the possible correlations that can be explained 
within these LHV models. For a precise formulation of the concepts of LHV description and the 
derivati on of the corre sponding Bell inequalities in the multi-party setting, we refer the reader 
to refs. jl26l 127 , 128 1 and for a brief revie w to ref. t l29til3Qil . Different Bell inequalities can 



also be regarded as entanglement witnesses \ 131 132k for different types of entanglement i 



m a 



multi-party entangled state. These witnesses can be quite useful to detect entanglement in the 
vicinity of graph states. In the following we will shortly review the results of L5(Xi 12411 on Bell 
inequalities and entanglement witnesses for graph states. 

But let us first consider an extension of the GHZ-argument that rules out a LHV de- 
scription of the spin statistics for GHZ-states, to the case of general graph states iiEii: The 
non-trivial graph state with two qubits is LC-equivalent to the singlet state and thus violates the 
original Bell inequality proposed by Bell in ref. 0]. For any connected graph state on more than 
three vertices any connected subgraph on three vertices a, b, c gives rise to a contradiction within 
any possible explanation of the observed correlations between spin measurements at different 
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particles. Consider, for example, the case where all three vertices a, 6, c are pairwise adjacent, 
i.e., {a, b} , {b, c} , {a, c} G E. Due to the non-commutativity of the spin observable algebra 
one easily computes that the product of the corresponding correlation operators Ka — cr^cr^". 



Ka — CT^O'^'' and Ka = cr^cr'^" yields 



(136) 



KnKhKr 



As discussed in sec. l2'2l these stabilizer elements provide con- 
straints to the four different measurement settings 



(/) m>f» = 1 (///) m>. 



(//) m^mf " = 1 {IV) - mlml^mlm 



1 



= 1 , 




where m° = ±1 and m% = ±1 denote the measurement out- ^'8- ~ connected sub- 

■r^u » • J- ■ graph on three vertices gives 

comes if the qubit at vertex a is measured in spin-x- or z- _ " 

rise to a violation of local real- 

direction. More precisely any LHV model, which assigns prede- 
termined('^^) values nix, ni, to measurement outcomes to z-and 

z-measurements at the different vertices in Na U Nf,U Nc with some probability, must be such 
that these assignments obey all the four equations (I)-(IV). Due to the commutativity for the 
multiplication of the measurement results the product of eq. (I), (II) and (III) gives 



(137) 



and thus a contradiction with eq. (IV). A similar argument holds for th e cas e where only two 
pairs of the three vertices a, 6, c are adjacent. Thus we have obtained 

Proposition 11 fNon-classicality of graph states). Any graph state corresponding to a connected 
graph violates local realism. More precisely, for a connected graph state with more than two 
vertices any connected subgraph on three vertices a, 6, c yields a contradiction when trying to 
explain the correlations between the different Pauli-spin-observables present in the reduced state 
Pq (on the subset A = Na U A^h U Nc) by means of some LHV model. 

More quantitatively one can compare (convex) functions of the correlations in a state ob- 
tained for different measurement settings with the results which can be found for these settings 
when considering all LHV models. Upper bounds to these functions give rise to so-called Bell 
inequalities, which the correlations arising from different measurement settings have to obey, if 
they can be explained by a LHV model. In fact for N parties and two dichotomies^) observables 
at each par ty the set o f states that allow for a LHV description, can be ch aracte rized by only one 
inequality jl27l 128 1. which is a generalization of the CHSH inequality jl34ll for the two-party 
setting. For the multi-party setting the different Bell inequalities can in general capture different 
types of entanglement according to different partitionings. For any non-trivial graph state Bell 



(*^) That means independently of the chosen measurement setting in (I)-(IV). 
(*®) I.e., the respective observables have two measurement outcomes ±1. 
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inequalities with three dichotomic measurements per site have been derived in ref. jsoll that are 
maximally violated by this graph state: 

Proposition 12 fBell inequalities for graph states). A (connected) graph state \G) with stabi- 
lizer S maximally violates {'^'^) the Bell inequality 



(138) 



where C is the maximum of the absolute value of the mean value (X^o-es f'^ifeen over all deter- 
ministic{^^) LHV models. In the case that C > 2^^^ the Bell inequality detects only states that 
are NPT with respect to every partitioning. 

ref. S also provides general rules for the computation of the value C, which a priori de- 
mands to check an exponentially (with N) increasing number of LHV models. Numerical results 
for graphs with up to 10 vertices show that Bell inequalities for rings and chains give rise to a 
large relative violation ^ for these states while the Bell inequalities for the corresponding GHZ 
states yield a small relative violation. 

This Bell inequality can be expressed in terms of an entanglement witness f'^) of the type 



(139) 

i.e., a state p violates the eq. ( I138> iff 



W ^ly-\G){G\, 



(140) 



(W) = tr(Wp) < 



To detect entanglement itself and not 'non-classicality' with respect to LHV descriptions, the 
constant C can be chosen according to a maximization over the smaller set of separable states 
only (''"). A concrete procedure to measure this witness in an experiment can be found by de- 
composing W into a sum of locally measurable operators. But, for an arbitrary graph state, the 
witness in eq. ( I139> seems to be deco mpos able only into an exponentially (with N) increasing 
number of local measurement settings 112411 . Toth and Giihne therefore proposed more practical 
entanglement witnesses and Bell inequalities for graph st ates that in ma ny cases, such as CSS 
states, can be evaluated in only two measurement settings 1 124t 135 , 136ll : 



r)I-e.,|E.g5(Gk|G)| =2^^. 

Due to convexity it suffices to consider deterministic LHV models, which assign definite values ±1 to 
all observables (with unit probability). 

(*^) An entanglement witness W is an observable with a positive or zero expectation value tr(VV/9) > 
for all separable states and a negative expectation value tr(VVp) < for some entangled states. In a 
measurement of this witness the entanglement present in the latter states can thus be detected ('witnessed'). 
(^'') Note that there exists non-separable states that still allow for a LHV description ISIl . 
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Proposition 13 fEntanglement witnesses for graph states). Let \ G) be a graph state correspond- 
ing to a connected graph. Then 

(141) := Iv-Ka- Kb 

is an entanglement witness for the \ G) that detects entanglement in the reduced state Pq (A — 
Na U NbU {a, b}) with only two measurement settings and thus can rule out full separability of 
the total graph state. The entanglement witness 

(142) W2 := (iV - l)ly - ^i^a 

detects genuine multi-party entanglement. IfG is M -colorable, then the evaluation of the witness 
eq. M42\ requires at most M local measurement settings. 

For further entanglement witnesses that are particularly robust against global white noise 
and that can be derived for s ome specia l cases such as GHZ states and linear cluster states, we 
refer the reader again to ref. 1 124 , ll35ll . In practice it is quite useful that the witness eq. (I142> 



for genuine entanglement gives also lower bounds to the fidelity {G\p\G) with the ideal graph 
state, when p denotes the outcome of a actual preparation procedure for the graph state. This 
provides an efficient method to verify that, in a given experiment, the entanglement is really 
present in a form that is sufficiently close to a desired graph state. An alternative approach 
to detect multi-partite entanglement, which is particularly suited for implementations in optical 
lattices and magnetic micro-traps, can be found in ref. 1 1 37|. 

8'2. Two-particle correlations and localizable entanglement. - In this section we consider the 
entanglement properties of the reduced state p\^'^^ of two qubits a and h that is obtained after 
tracing out or disregarding the information about the remaining particles c ^ V \ {a, 6} in a 
graph state \G). As discussed in sec.|6|the reduced state 

(143) P{a,b} = J2 

is essentially given by those stabilizer elements a E S^a.b} that act non-trivially only on the 
qubits a and b. Let us first examine the classical correlations between two non-isolated vertices 
a and & in a graph state | G) 

(144) Q^;f (G|< ® a^^lG) - (G|<|G) (G|a^^|G) t,j - 1,2,3 . 

Note that these correlation functions Q°^j only depend on the reduced state Pq^'''^ of the particles 
a and b and are given by 

(145) Qi,-^[PG ^^''0~[o if ' 
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since the expectation values e.g. {G\(t^\G) = \i{pQ(Tf) = tr(ilaO'f) vanish for both (non- 
isolated !) vertices a and h. For the maximal classical correlation between two vertices a and b 
in a graph state | G) 

(146) Q;J,L = ^ max^ IQ^"! 

we find that it vanishes whenever the neighborhoods Na \ b and A*";, \ a of the two vertices with 
respect to the remaining graph are non-empty and distinct. 

Pro/joi/f/on 74 (Two-party classical correlation). For two non-isolated{^^) vertices a,b G V 
in some graph G = (V, E) we have 

n47^ n-" ^ ^ ^ {Na\b),{N,\a)^{band{N,\b)^{N,\a) 

^ ' ^'""■^ 1 1 otherwise 



Proof: According to eq. (I145> we have to show that all stabilizer elements ct e 5 have support 
supp((T) on more than the two vertices a and h iff [Na \b) ^ {Nt, \ a) and {Na \ b), {Nb \a) 7^ 0. 
Since all generator elements are generated by different combinations of the correlation operators 
Kc (C C V), the sufficiency can be derived as follows. In order to generate a stabilizer 
element cr = Hcgc ^^^^ supp((T) C {a, 6} at most the two correlation operators on vertices 
a and b can be considered, i.e., C C {a, b}, because any other correlation operator Kc for 
c E V \ {a, b} leads to a non-vanishing Pauli operator ctJ or ajj on the central vertex c outside of 
{a, b}. Moreover any support supp(A'a) \ {a, 6} 7^ or supp{Ki,) \ {a, 6} 7^ of the correlation 
operators Ka and A'b outside of {a, 6} can only be compensated if both vertices have the same 
neighbors outside of {a, b}, i.e., (Na \ b) ~ {Ni, \ a). In this case we are left with the two 
possibilities 



(148) KaKk 



if {a,b}eE 

,6 



a-al if {aM^E 



□ 

Although a graph state might contain some non-vanishing classical two-party correlations it 
does generally not include any entanglement between any two qubits, unless -as discussed below- 
the remaining parties are allowed to assist the revealing of such entanglement. 

Proposition 15 (Two-party 'quantum correlation'). For any vertices a,b € V the reduced 
state P(^'''^ of some graph state \G) is separable unless the graph contains the isolated edgei^^) 
{aM- 



(^^ ) If one vertex is isolated all correlation functions vanish, since the corresponding state is a product state 

e.g. |G) = |+)1G\a)^'\^ 

{^^) If G contains an isolated edge {a, 6} the state \G) decomposes into a pure Bell state on the vertices 
{a, 6} and some other graph state on the remaining vertices. 
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Proof: If G does not contain the edge {a, 6} as an isolated edge then according to Propo- 
sition|8]the reduced state p^'^^ is either the rank- 2 -projector Pq'^^ = \ {Pi2 + '^'^ Pi2<y^) 
(Pi2 ■■= |G[{1, 2}]) (G[{1, 2}]|) for some set C = {1}, {2}, {1, 2} or it is a rank-4-projector and 
thus the maximally mixed state p^" ''^ = j^ab- In the latter case the above statement is trivial. 
That also the rank-2-projectors correspond to separable states can be derived from the fact that 
they are PPT (which means that their partial transpose is positive) according to eq. ( I198> . which 
is a sufficient condition for separability in 2 x 2-systems. |--| 

Although there is per se no entanglement between arbitrary two particles in a connected 
graph state, such entanglement can be revealed between any two parties a and h, if the remaining 
parties are allowed to perform local measureme nts >36lll. More generally, the notion of localizable 
entanglement 'LE°'^ {p) was introduced in refs. Il38lll39lll4(ill for multi-spin states p, defined as 
the maximal amount(^'^) of entanglement that can be created (or localized), on average, between 
two spins at position a and h by performing local measurements on the other spins. For a general 
state p it has been shown in ref. Iil38al that the localizable entanglement LE"'' is related to the 
maximal classic al correlati on Qmax and the entanglement of assistancei^^) AE"'' as measured by 
the concurrence jl41 142 1 : 



(149) QtM < ^^"'(P) < AE"^(p) 



Despite of the separability of the reduced states, measurements on the remaining particles can 
nevertheless create maximal entanglement between any two vertices in a 'connected' graph state, 
which corresponds to maximal localizable entanglement in this case. 

Proposition 16 (Localizable entanglement). Consider any two vertices a,b G V in a graph 
state correspondin g to a connected graph G = (V, E). In all measurement branches of the 
following protocol h3al a maximal entangled state is created between the vertices a and b: 

1. Choose any path (oq = a, ai, . . . , a„_i, a„ — b) connecting the vertices a and b. 

2. Measure the spin of all vertices except ai in z-direction. 

3. Measure the spin of all vertices at for i ~ 1, . . . ,n — 1 in x-direction. 

Thus the localizable entanglement LE"'' of a 'connected' graph state is maximal. 

Proof: If the graph is connected there exists a path connecting any two vertices a,b £ V. 
According to Proposition the CTz -measurements in the second step simply remove all vertices 
but those on the path (ap = a, ai, . . . , a„_i, a„ = b). Note that the different local unitaries 
according to the different cr^ -measurement outcomes consists only of -operators and thus do 
not alter the measurement direction of the subsequent (Ta; -measurements in the 3. step. The 
same proposition also implies that a sequence of measurements on the inner vertices (i ~ 



(^^) E.g. in terms of its concurrence. 

{^^) The entanglement of assistance extends the concept of locahzable entanglement in that it allows also 
joint measurements to be performed on the other spins. 
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1, . . . , 71 — 1) of this path removes these inner vertices but keeps the connectivity between the 
two neighboring particles Ui-i and a.i^i in this chain. E.g. an CTj: -measurement at vertex ai of 
the initial chain Gq yields the shorter chain Gi = (Go \ ai) U {qq, 02} on the remaining vertices 
(tto = a, 02, ... , fln-ij CLn = b). The LC-unitaries corresponding to the different measurement 
outcomes again do not rotate the subsequent measurement directions. Thus we 'inductively' 
arrive at the final graph G„_i corresponding to the maximally entangled state \Gn-i) of eq. ( I23t . 

□ 

Concluding, the results of this section indicate that the entanglement present in graph states is 
not based on bi-partite 'quantum correlation' but the entanglement is rather delocalized among 
all particles. 

83. Quantifying entanglement. - We have seen that graph states are entangled quantum states 
that exhibit complex structures of genuine multi-particle entanglement. The main aim of this 
section is to apply the quantitative theory of multi-particle entanglement to the study of corre- 
lations in graph states. Needless to say, despite considerable research effort, there is no known 
computable entanglement measure that grasps all aspects of multi-particle entanglement in an 
appropriate manner, if there is any way to fill such a phrase with meaning. Several entangle- 
ment measui^formi^ have yet been suggested and their properties studied 
l l4ill44lll4ill4illHll4ill^ 



In this section the underlying measure of entanglement is taken to be the Schmidt measure 

1I143I1 . which is a proper multi-particle entanglement monotone that is tailored to the character- 
ization of such states. As holds true for any known measure of multi-particle entanglement, its 
computation is exceedingly difficult for general states, yet for graph states this task becomes fea- 
sible to a very high extent. We present various upper and lower bounds for the Schmidt measure 
in graph theoretical terms, which largely draw from stabilizer theory. These bounds allow for an 
evaluation of the Schmidt measure for a large number of graphs of practical importance. 

The Schmidt measure has been employed to quanti fy th e degree of entanglement, as a gener- 
alization of the Schmidt rank in the bi-partite setting jl43ll . This measure is sufficiently coarse 
to be accessible for systems consisting of many constituents and to allow for an appropriate 
discussion of multi-particle entanglement in graph states. 

Any state vector e Hi (g) ... (X) H^r of a composite quantum system with N components 
can be represented as 

R 

(150) IV') =^Cd^->®---®l^f), 

i=l 

where G (D for i = 1 , . . . , i?, and | tA" ) € H„ forn = 1 , . . . , . The Schmidt measure associated 
with a state vector | is then defined as 

(151) EsdV)) = log2(i?mm), 

where Rmm is the minimal number R of terms in the sum of eq. ( I150> over all linear decompo- 
sitions into product states. It can be extended to the entire state space (and not only the extreme 
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points) via a convex roof extension(^^). It should be noted that the Schmidt measure E5 is 
discrete, e.g. in the case of two-level systems H„ = it can only take the values 

(152) Esim G {log2(77i) I m = 1, . . . , 2^} , 



and thus fails to be continuous, which requires some care when extending the measure to general 
mixed states via the convex roof construction. However, since the set of graph states is discrete 
itself, the case does not occur that two graph states become arbitrarily close in terms of some 
distance measure but not with respect to their Schmidt measure. In any case the Schmidt mea- 
sure is a general entanglement monotone with respect to general local operations and classical 
communication (LOCC), which typically leave the set of graph states. If the Schmidt measure of 
a state vector \ip) is evaluated with respect to a partitioning {Ai, Am), it will be appended, 

(153) e(,^--^->(|V)), 

in order to avoid confusion. For a graph G = {V, E), the partitioning with M = N and A„ = 
{n} will be referred to as finest partitioning. If no upper index is appended to the Schmidt 
measure, the finest partitioning will be implicitly assumed. 

Among the properties that are important for the rest of this section are the following jl43l 

[l50il : 

Proposition 17 (Schmidt measure). 

(i) E5 vanishes on product states (only), i.e. 

(154) Es(|V'))=0 ^ \4') = \ip^)(^...(^\ii^'') . 



(ii) E5 is non-increasing under SLOCC Il43r . i.e. 

(155) |V> —SLOCC IV'') =^ Es{W))<Es{m 
Similarly for the LU-equivalence we find 

(156) |V)^Lu|V'') =^ Es{W))^Esm) ■ 



(^^) Note that every positive function -Epme defined on the set of pure states that vanishes exactly on product 
states (see property (i)) and is non-increasing under SLOCC (see property (ii)), can be extended to the entire 
state space by 

E{p) :=inf{^A,£;p„re(|V'.» I p = ^ A,|V'«){V'^I with A, > ^ A, = 1} . 

i i i 

Note that infimum of the average entanglement can be regarded as taken with respect to all preparation 
procedures of the mixed state. The entanglement measure E can be proven to (a) vanish exactly on all 
separable states and is (b) convex as well as (c) non-increasing under SLOCC operations. 
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(iii) E5 is non-increasing under a coarse graining of the partitioning, i.e. 
(157)(^i,...,Am) < {B,,...,Bm') =^ Ef'■-^"\\i^))>Ef'■-''''"\\^^)) . 

Thus if two components are merged in order to form a new component, then the Schmidt 
measure can only decrease. 

(iv) E5 is sub-additive, i.e., 

(v) For any bi-partition [A, B) E5 coincides with the Schmidt rank Egd?/;)) ~ SRAdV")) = 
log2(rank(tr^[|V') ("01]))- In particular, E5 is additive within a given bi-partitioning, i.e., if 
A = Ai U and B = Bi U B2, then 

(159) E^^'^)(|^i)®|V'2))=E^^-^^)(|Vi))+E^^-^^)(|V2». 

It should be noted that for general pure states of multi-partite quantum systems the Schmidt 
measure is - as any other measure of multi-partite entanglement - exceedingly difficult to com- 
pute I37I1 . In the following we will provide lower and upper bounds for the Schmidt measure 
of graph states in graph theoretic terms, which will coincide in many cases, and will then apply 
these rules, and calculate the Schmidt measure for some of graphs and graph classes that are of 
interest for applications.. 

Let us first derive a lower bound to the Schmidt measure, namely the maximal Schmidt rank 

max 

(160) SRn,ax(G) max SRa(G) • 

If one maximizes over all bi-partitionings [A, B) of a graph G = (V, E), then according to 
eq. (I157> one obtains a lower bound for the Schmidt measure with respect to the finest partition- 
ing. 

Since the Schmidt ranks SRa(V') for the different bi-partitions are already entanglement 
monotones with respect to {A, i3)-local SLOCC-operations, it is straightforward to see that 
SRmax is a proper though also discrete entanglement measure(^^) that captures the notion of max- 
imal bi-partite entangleme nt co ntained in a multi-party entangled state. In fact this measure was 
also considered in ref . fl6l Il5lll in the context of an efficient simulation of a quantum algorithm 
on a classical computer. There it was shown that, if through a pure-state quantum computation 
all underlying pure states of the A^-qubit quantum register have polynomial bounded 'Schmidt 
rank', i.e., 2^'*™"' < poly(iV), then this quantum computation can be efficiently simulated by a 
classical algorithm requiring only polynomial increasing memory space and computational time. 



(^®) SRmax has properties (i) - (iv). 
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The classical algorithm allowing for a simulation of quantum computations with only slightly 
entangled quantum register uses an efficient decomposition of the pure state as in eq. dlSOI l that 
allows to directly read off and manipulate the Schmidt coefficients for different bi-partitionings 
through the quantum computation. This algorithm proved to be very useful also for an efficient 
simulation of the dynamics in one-dimensional quantum many -body systems due to its tight con- 
nection to existing methods such as density matrix renormalization group and 



opens a way to even simula te sy stems in higher dimensions, where these methods so far have not 
seemed to be very suitable 1741 17511 . 

For graph states the maximal Schmidt rank SR,nax actually coincides with continuous en- 
tanglement measures such as the entropy of entanglement or the purity of the reduced density 
matrices. From eq. ( II 17> one can compute that the entropy or the purity of the reduced density 
matrices for \G) according to a bi-partition [A, B) gives 

(161) SR^(G) = -tr[p^log2(p^)] - log2(tr[(p^)2]) . 

This again expresses the fact that, for a non-empty graph, |G) is a 'maximally' {A, i3)-entangled 
State vector with 2 s Schmidt coefficients. 

Finally, the Schmidt rank of a graph state is closely related to error correcting properties of 
a corresponding graph code (see sec. I4'3> . Let ^ be a partition, according to which \G) has 
maximal Schmidt rank. Then, according to ref. , choosing a subset X C A, the graph code, 
which encodes an input on vertices X in output on vertices Y = V \ X according to G, detects 
the error configuration E = A \ X , i.e., any errors occurring on only one half of the vertex set E 
can be corrected. In particular, all strongly error correcting graph codes in ref. must have 
Schmidt measure N/2. The following proposition gives at least a sufficient condition 
when a partition has a maximal Schmidt rank with respect to the corresponding bi-partite split. 

Proposition 18 (Maximal Schmidt rank). A sufficient criterion /or a bi-partite split {A, B) to 
have maximal Schmidt rank is that the graph Gab corresponding to the edges between the 
partitions contains no cycles, and that the smaller partition contains at most one leaf with respect 
to the subgraph Gab- If Gab is not connected, then it is sufficient that the above criterion holds 
for every connected component of Gab- 

Note that a leaf is a vertex of degree 1, i.e., a vertex to which exactly one edge is incident 
a. 

For the upper bound we consider a sequence of local projective measurements that finally 
completely disentangles the state vector \ip) in each of the measurement results. Let m denotes 
the number of measurement results with non-zero probability. Clearly, any of the states resulting 
from the different measurement outcomes is a product state and thus the whole measurement 
procedure gives rise to a decomposition of the initial state 10) as in eq. ( fBOl l. Thus we obtain 
the upper bound 

(162) Esim < log2(m) . 

In particular, for any sequence of measurements in the Pauli basis ax, (Jy or cr^ that yields an 
empty graph, the number of local measurements in this sequence gives an upper bound on the 
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Schmidt measure of the corresponding graph state. This is because -apart from the trivial case 
of a CTj; -measurement at an isolated vertex-, both measurement results ±1 of a local Pauli mea- 
surement are attained with probability 1/2 and yield locally equivalent graph state vectors \G') 
and \G"). More generally we find 

(163) E5(|G")) < Es(|G)) < Es{\G')) + 1 . 



In the following we will call the minimal number of local Pauli measurements to disentangle 
a graph state its Pauli persistency PP(G). The notion of persistency was introduce in ref. in 
the context of general projective measurements in order to study the stability of entanglement in 
cluster states with respect to local measurements. Since each measurement simply deletes all 
edges incident to a vertex, any subset V' QV of vertices in a graph G, to which any edge of G 
is incident, allows for an disentangling sequence of local measurements. In graph theory those 
vertex subsets are called vertex covers VC(G). Thus we have found the upper bounds 



(164) 



Es(|G)) < PP(G) < VC(G) 



For graphs with many edges a combina- 
tion of and ay will give better bounds than 
restricting to cr^ measurements only. For ex- 
ample, according to the measurement rules in 
Proposition Q any complete graph (in which 
all vertices are adjacent) can be disentangled 
by just one ctj, -measurement at any vertex (see 
fig. I16> . As we have seen, this corresponds 
to the fact that these graph states are LC- 
equivalent to the GHZ-type graph states (see 
also sec. 14' 1> . in which every vertex is adja- 
cent to the same central vertex. 

Let us briefly summarize the relevant 
bounds for our further considerations in a 
proposition. 




Fig. 16. - A single ctj, -measurement at an arbitrary 
vertex in the complete graph No. 7 suffices to disen- 
tangle the corresponding state. Similarly, a single (Jz- 
measurement at the central vertex in the graphs No. 
1-6 or a single -measurement at the non-central 
vertices is a disentangUng measurement. This is due 
to the fact that all graphs (No. 1-7) are locally equiv- 
alent by LC-unitaries, which transform the measure- 



ment basis correspondingly. 
Proposition 19 f Bounds to the Schmidt measure). 

For any graph state \G) the Schmidt measure Ea is bounded from below by the maximal Schmidt 
rank SR„,ax and from above by the Pauli persistency PP or the minimal vertex cover VC, i.e. 



(165) 



SRmax(G) < Es{\G)) < PP(G) < VC(G) . 



An application of the LC-rule (see Proposition |3, of course, does not change the Schmidt 
measure. But also other local changes to the graph, such as the deletion of edges or vertices, 
have only a bounded effect on the Schmidt measure isTil : 



Proposition 20 fEdge-A^ertex rule). 
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• By deleting or adding edges between two vertices of a graph G the Schmidt measure of the 
resulting graph G' can at most decrease or increase by 1; 

• If a vertex is deleted, the Schmidt measure of the resulting graph G' decreases, but at most 
by 1. 




Fig. 17. - An example for the (4, 5, 3)-cluster state and the graph corresponding to the adjacency matrix 
(see eq. <116H for a bipartitioning [A, B) with maximal Schmidt rank SR^ — SRmax- Here the 
vertices in A are depicted by small boxes ■ . 

Let us now apply these findings to evaluate the Schmidt measure for some important classes 
of graph states: 

Proposition 21 fExamples). 

• The Schmidt measure for any multi-partite GHZ states is 1. 

• The Schmidt measure of a ID-, 2D-, 3D-cluster state is [^J- 

• The Schmidt measure of an entangled ring with an even number N of vertices is given by 
N/2. 

• The Schmidt measure of a tree is the size of its minimal vertex cover VC. 

In all these cases, the upper and lower bounds, i.e., the maximal Schmidt rank and the Pauli 
persistency, coincide (^^) Although we refer to ref. 113 711 for a detailed proof, for a even ring or 
a cluster state with at least one side of 'even length', we have indicated the partition A with 
maximal Schmidt rank SRa = SRmax (see Propositionll8> in fig.[n]or fig.ll9lrespectivelv. 



{^'') Since the maximal Schmidt rank for any state can be at most [-^J , the Schmidt measure of those cases, 
where upper and lower bound coincide, is also bounded by this number. Thus all graph states which allow 
for a determination of the Schmidt measure along these lines have Schmidt measure of at most [-^J . 
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# Smallest Bipartition =5 # Minimal Vertex Cover = 4 

Fig. 18. - The graph No. 1 represents a tree. Its bi-partitioning (^4, B), for which in graph No. 2 the vertices 
in A are depicted by disks, neither is a minimal vertex cover nor yields maximal partial rank. Instead the 
set of vertices A, represented by large disc's in graph No. 3, is a minimal vertex cover with maximal partial 
rank. 



Fig. gives an example for a tree for 
which the Schmidt measure does not coincide 
with the size of the smaller bi-partition, the up- 
per bound according to Propositionl21l 

We have also computed the lower and up- 
per bounds to the Schmidt measure, the Pauli 
persistency and the maximal partial rank, for 
the non equivalent graphs in List A and B in 
Table IVl They are listed in the corresponding 
tables in Sec|7| 

For connected graphs the Schmidt rank 
cannot occur for any bi-partite splitting 
[A, B), since this would correspond to an empty graph Gab between the partitions. Because 
the rank index is invariant under permutations of the partitions according to graph isomorphisms 
it provides information about whether two graph states are equivalent under local unitaries plus 
graph isomorphisms as treated in sec.|7| But note that graph number 40, 42 and 44 are examples 
for non-equivalent graphs with the same rank index. Nevertheless, comparing the list of Schmidt 
ranks with respect to all bi-partitions in detail shows that no permutation of the vertex set exists 
(especially none which is induced by a proper graph isomorphism on both sides), which would 
cause a permutation of the corresponding rank list, such that two of the graphs could be locally 
equivalent. 

For 295 of 995 non-isomorphic graphs the lower and upper bound differs and that in these 
cases the Schmidt measure also non-integer values in log2{l, 2^} are possible. Moreover 
note that only graph number 8 and 19 have maximal Schmidt rank with respect to all bi-partite 
splits. Entanglement here is distributed symmetrically between all parties, which makes it "dif- 
ficult" to disentangle the state by few measurements. From this one can understand why the gap 
between the lower and upper bound occurs in such cases. As discussed above, from all graph 



No. 1 No. 2 




Fig. 19. - Graph No. 1 is an entangled ring on 18 ver- 
tices. Graph No. 2 represents the graph correspond- 
ing to the adjacency matrix (see eq. <1 16t ) for 
a bipartitioning {A, B) with maximal Schmidt rank 
SRa = SRmax. The vertices of partition A are de- 
picted by boxes. 
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22 23 24 25 26 27 2S 29 30 



codes with less than 7 vertices only these two are candidates for strongly error detecting graph 
codes introduced in ref. 

In the remainder we will briefly discuss two exam- 
ples which were already introduced in sec.|3|in the con- 
text of quantum error correction and the one-way model 
for quantum computation, namely the graph that is used 
to realize the QFT on three qubits in the one-way quan- 
tum computer The vision behind this is to flesh out the 
notion of entanglement as an algorithmic resource, as it 
has been put forward in ref. lIsaBsilillil. 

Example 1: Concatenated [7, 1, iyCSS-code. 

As discussed in sec. 14 31 the graph G depicted in 
Fig. 1201 represents an encoding procedure for the con- 
catenated [7, 1, 3]-CSS-code. The corresponding graph Fig- 20. - Resource graph state for the con- 
state has Schmidt measure 28. For encoding, seven a^- catenated [7, 1, 3]-CSS code, 
measurements at all vertices of the inner square except o have to be performed. In this resulting 
graph obtained without measuring the vertex o represents the resource for the alternative encod- 
ing procedure. It has maximal Schmidt measure 25, whereas the corresponding and 1 code 
words have Schmidt measure 24. They can be obtained with probability 1/2 from |G") by a 
tTz -measurement at the vertex o. 

Example 2: Quantum Fourier Transform ( QFT) on 3 qubits. 

The graph No. 1 in fig. |^ is a simple example of n° ' 
an entangled graph state as it occurs in the one-way 
computer of refs. I32I I33I I34I Esll . This specific ex- 
ample represents the initial resource (part of a cluster) 
required for the quantum Fourier transform QFT on 3 
qubits lillil 13 El. It has Schmidt measure 15, 
where the partition 

A = {2, 4, 7, 9, 11, 13, 15, 18, 20, 22, 24, 26, 28, 30, 32} 

is a minimal vertex cover with maximal Schmidt rank. 
In the process of performing the QFT, all vertices except 
the output vertices 5, 16, 33 are measured locally. Dur- 
ing this process, the entanglement of the resource state 
(with respect to every partitioning) can only decrease. 
Similar as with the graph state vector |G") obtained from 
fig. |20| graph No. 2 represents the input-independent 
resource needed for the essential (non-Clifford) part of 
the QFT protocol lillilliHil. It has Schmidt mea- 
sure 5, where the partition A = {2, 9, 10, 11, 15} now 
provides a minimal vertex cover with maximal Schmidt 
rank. 

Fig. 21. - Graphs associated with the QFT 
on 3 qubits in the one-way quantum com- 
puter. The boxes denote the input (left) and 
output (right) vertices. 
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9. - Weighted graph states 

In this sectio n we extend the concept of a graph describing some state to the class of weighted 
graph states We return to the definition of graph states in terms of the underlying inter- 

action pattern. In sec. 12' H it was shown that any interaction pattern, for which the temporal order 
of the interactions between the particles is irrelevant and thus can be described by a graph, can 
only contain interactions Hab that are up to local z-rotations Ising-type interactions or likewise 
are given by the phase gate Hamiltonian Hat = ® In order to describe the state 

by a simple graph we have fixed the interaction time for the particles to a phase tt in the previous 
sections. We will now allow the particles to interact according to the same Hamiltonian Hab but 
for different interaction times (pat- This corresponds to the situation of a disordered system as it 
occurs e.g. in a spin glass or semi-quantal Boltzmann gas as described below. The interaction 
pattern is now summarized by a weighted graph, in which every edge is specified by a phase ipab 
corresponding to the time the particles a and b have interacted. The weighted graph state \G) is 
thus given by 

(166) \G)= II (7„fc|+)^ 

{a,b}eE 

where the operations Uab depend on the interaction phases ipab- 

(167) Uab e-'^-^- = e-'^i^'-'^'^M^'-"') . 

The corresponding adjacency matrix T for this weighted graph collects the weights Tab ~ <fab = 
(fiia and gives rise to a concise decomposition of the weighted graph states with respect to the 
standard basis \W)^ = \WiY ■■■\Wn)^ {W = {Wi,...,Wn) € F^) jsil: 

(168) \G) =2'f n Uab\W). = 2-f e^'^'^-^-'^m. . 

WCV{a,b}eE WCV 

This easily follows for each basis state | W) ^ by induction over the involved interaction unitaries 
IJab using the fact that 

(169) UabWaTW^)^ = e'^"^^-^''|W^a>lM^6>'' = e'^ ('^"^"''^''+^''^''"^" ^ | H^a>1 W'b) ^ . 
We remark that for any weighted graph states |G), the set of states 

(170) \W)=a^\G) 

form a basis for (C^)^. This can most easily be seen by realizing that aY commutes with [/„;,, 
and cr^ 1+)^ gives rise to a orthonormal basis consisting of pure product states. The application 
°f W.{a b}eE ^o-b transforms this orthonormal basis into a basis of weighted graph states, where 
each of the basis states has equivalent entanglement properties and the states are connected by 
local unitary operations. 
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In contrast to a straightforward extension of the interaction picture to the case of weighted 
graph states, no such generaHzation of the stabilizer formaUsm in terms of generators within 
the Pauli group is possible. This implies that many results obtained in the previous sections 
do no longer hold for weighted graph states. Nevertheless one is still able to show liHl that 
entanglement in a weighted graph state is closely related to the connectivity properties of the 
underlying graph. More precisely, the weighted graph state is entangled (non-separable) with 
respect to the bi-partition {A, B) iff there are non-vanishing interaction phases Lpab > for some 
a £ A and b E B between these partitions. This is because the rank of the reduced density 
matrix as computed in Proposition 1221 in sec.|5]is larger than one(^'*) in this case. Moreover, 
between arbitrary two parties A and B entanglement can be created (localized) by performing 
local operations on the remaining vertices (i.e., inV \ {AU B)) iff the parties are connected by 
a path with non-vanishing interaction phases. The localizable entanglement is not zero in these 
cases, since the protocol in Propositionll6lof sec. l8'2l can also be applied to weighted graph states 
in order to reveal entanglement between some vertices a E A and b € B liUl . 

Moreover for the following analysis of entanglement present in weighted graph states it is 
crucial that the reduced density matrices of these states can still be determined efficiently: 

Proposition 22 (Reduced state for weighted graph states). Let A (- V be a subset of vertices 
for a weighted graph G and B = V \ A the corresponding complement in V. Then the reduced 
state Pq ~ tr b{\G) {G\) is given by 



(171) 



n Uab-p^ui, 

{a,b}£EA 



where Ea = E H {A x A) denotes the set of edges within A and 
(172) Pg= E Ca,aMi)UA2\- 



Ai,A2CA 



The matrix elements are 

CAiA2 

(173) 



9Af E 



i(Ai-A2)-T'-B' 



B'CB 

beB 



n 



-{A,-A-2)-T', 



Here T'^^ denotes the b-th column of the matrix T' = T ab representing the edges {a, 6} with 
weights ipab between the partitions A and B (see eq. M16\ ). 

Proof: From elementary facts('''^) about the partial trace Pq = 1-s:b{\G){G\) it follows that, 
for the computation, we might as well apply the interaction unitaries Uaa' in eq. Mbbl acting 



(^*) More precisely, for any (fab 7^ the rows according to index and a = {a} differ: We find for the 
matrix elements ^ / g^, since = Caa = 1 and = C„0 = e'^'-e"- '^""/^ TlbeN^ cos[^]. 
(5^) I.e., trslC^ ® (7^p75^ (g. (U'^)'^] = C\.ib(p) D for arbitrary unitaries U and matrices C, D. 
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solely on vertices in A after the partial trace and neglect all unitaries acting only on B. Thus 
Pq in eq. ( I171> is the reduced density matrix corresponding to the graph Gab with adjacency 
matrix T' — Tab- In order to determine Pq we split up the basis vectors \ W) in eq. ( I168> . e.g. 
IVFi)^ = \ Ai)'^\Bi)^ , according to the partitioning {A,B) and compute 



WuW^CA 



(174) 



2N 



A2)-r'-B' 



\Ai)t{A,\ 



The matrix element CaiA2 can be further simplified as stated in eq. ( I173> : 



(175) CA.A2 = E n e^'^^-^^y^'^ = R (l + ^^^^^^^^'^^ 

(176) = e'2 '^'•es^"" -^/-b I I cos 



B'CB beB' 

beB 



beB 



using i(l + e'"^) = cos(|). □ 

Since the entanglement properties between the partitions {A, B) are invariant under LOCC, 
we can disregard the unitaries remaining on A and thus directly examine pQ. From eq. ( I173> 
one sees that the total effect of the interactions with particles in B on the matrix elements (co- 
herences) CaiA2 consists in multiplying the effects for each individual particle b G B. Thus 
Pq can alternatively be obtained as the Hadamard product(^") of the reduced states Pq^ due to 
sole effect of particle b E B, i.e., according to the induced graph Gb '■= G[A U b] on the vertex 
6 and all vertices within A. Note that the computational effort to calculate the reduced density 
matrix scales exponentially in \A\ but only linearly in the number |-B| of particles in the remain- 
ing system, while for general pure states an exponential increase (with N) in computation time 
and memory cost is required for both parties. Hence all quantities that depend on the reduced 
density operator can be determined efficiently. For instance, from p^ of one and two qubits, we 
can calculate all two-point (and also higher order) correlation functions Q"'', 

(177) Q»f =«a^)-«)(a^) ^, J = 1,2,3, 

the entanglement of formation(^^) betw een p airs of particles, as well as lower and upper bounds 
on the localizable entanglement LE*^^ 113811 that was already discussed in sec. l8'2l for simple 



(^'') I.e., componentwise matrix multiplication of matrices written in the standard basis. 
(®^) For a two-qubit mixed state p"*" the entanglement of formation Ef(p) = f(C{p)) is related to the 
concurrence C(p) := max{0, Ai — A2 — A3 — A4}, where \i are the eigenvalues of the Herm itian matrix 
(y^((7y(TyP*(Ty(Ty),/p)^''^ in decreasing order, by some monotonically increasing function / Il54ll . 
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graph states. Note that the maximal classical c orrel ation Q^j^^x between two particles is given 
by the largest singular value of the matrix Q"j' Il38ll . Recall that the localizable entanglement 
LE"'' is the maximum amount of entanglement that can be established between a pair of particles 
a, b, on average, by performi ng lo cal measurements on all other particles. More over t he relation 



Q°t^ < LE"'' < AE"'' holds {iM, where AE"'' is the concurrence of assistance \U, 

With an efficient calculation of the reduced state we can also determine the entwpyi^^) of 
bi-partite entanglement 

(178) Sa{\G)) S{pi) EE -tr[p^ log2(p^)] 

between a small number of verti ces A and the rest, as well as the multi-partite entanglement 
measure Emw pro posed in ref. jl47l 155 1. This measure Emw is given by Emw = 2[1 — 



tiMW pro poseg m i 
tr(pg;')] Il47lll5j 



This method for calculating reduced density matrices can be readily extended to the case, 
in which the particles are initially prepared in an arbitrary product state (^^gy (aa|0) + /3a|l)) 
instead of |+)^ or even to slightly entangled (^■^) initial states. Similarly the VBS-picture can 
be generalized to all states produced by the interaction Hamiltonian Hab acting on arbitrary 
product input states • • • |(/)Ar)^. However, further modifications of the involved states are 
required. The (unnormahzed) VBS-like state is of the form \^) — ^ \Xa^bi) with \xa'b^) = 
Uaiy I ^V^y I V^)'' , where |0) = a|0) + and| ^) ^|0) + v^|l). But, for the 
sake of simplicity, in the remainder of this article we restrict to weighted graph states, i.e., states 
arising from input states |+)^. 

Since the class of weighted states comprises particles interacting for different interaction 
times, it provides an interesting model for the study of the entanglement dynamics in many- 
particle systems. In the remainder we will briefly review a few results about the static and 
dyna mic entanglement properties of spin lattices (^*). For studies on spin gases, see refs. Isil 



153|]; these results are based on the efficient method to calculate the reduced density operators of 
a small number L < 10 of arbitrary spins as described above. 

In harmonic systems Jsii Issll . it may be remarked, also higher-dimensional systems and 



C'^) Note that for pure bi-partite states the entanglement of formation and distillation both are given by the 
entropy of entanglement Sa ■ 

(^^) I.e., pure states with small Schmidt rank with respect to the respective bi-partition. 
^64 -J ji^g investigation of entanglement properties of strongly interacting many body systems has proven to 
be a fruitful approach. Clearly, the ground states of interacting many-body systems at zero temperature are 
correlated. These correlations are not only reflected by scaling laws for two point correlation functions: In 
fact, it turns out that characteristic scaling laws concerning ground state e ntanglement can be estab lished 
reminding of the behavior of such two point correlation functions I53.l53lfl38lll56l [l5l[l5i[l5l. This 
observation refers on the one hand to entanglement properties of two distinguished constituents of a many- 
body system. On the other hand, it holds for block entanglement of a number L of consecutive constituents 
and the rest of an infinite chain js^.'lSS*!. Notably, the specifics of the scaling of entanglement were shown 
to indicate quantum phase transitions 1 138, 156, 157, 158, 159]. Long-range correlations can even be found 
in systems wi th gapped Hamiltonians in the sense of a divergent characteristic length of the localizable 
entanglement ll39lll6Cll . 
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systems of non-integer dimension can be studied, leading in particular to 'area-theorems', i.e., 
statements on the relationship between the degree of entanglement of a distinguished region of 
the fuU lattice system and its boundary area. These statements that gapped harmonic systems 
indeed imply the validity of an area theorem hold true on harmonic systems even defined on 
general graphs 

El. 

Let us start our considerations with the case where the graph has some lattice structure. Thus 
we consider N spin 1/2 systems (qubits) with pairwise interactions, described by an Ising-type 
Hamiltonian (see eq. ( I19» 

(179) H = Y, b)^ii- ® (1 - <^l)- 

a<b 

We assume that the spins are arranged on a d-dimensional lattice with fixed geometry and are 
initially completely polarized in .T-direction, i.e., j^Po) = |+)^- As indicated in sec.|5]the meth- 
ods we develop can also describe disordered systems with random coefficients f{a,b) and can 
take arbitrary (product) input states into account. We are interested in (entanglement) properties 
of the state 

(180) :=e-**^|*o)- 

We consider the situation where the coupling between spins obeys a certain distance law, in the 
sense that the coefficients /(a, b), describing the strength of the coupling, only depend on the 
distance Vab ■= \\a — b|| betwe en particles a and b, i.e., /(a, b) = f{rnjX In the example of ions 
stored in microtraps jl61 , 162 1 one finds for instance f{rab) = \ 162 1. 



' ab 

In the following the (bi-partite) entanglement between blocks of a small number L < 10 of 
neighboring spins as measured by the entropy of entanglement will be denoted by Sl = S{pl)- 
Clearly, < S'l < L, where Sl — L indicates maximal entanglement between the blocks. 
For different distance laws, we have investigated the scaling of block-wise entanglement Sl in 



ref. II153II and observed the following: 



The maximal two point correlations in a spin chain decay slower than exponential for all 
power laws f{rab) oc rl ^". Therefore, the correlation length ^ and also the entanglement 
length £,E diverge lll38ll . This indicates long-range quantum correlations for all power 
laws, as we find that only exponential fall-off functions /(a, b) = g"'*' "'' lead to a finite 
correlation length. 

In the limit N ^ 00 and L ^ 00, entropy Sl saturates as a function of L for power 
laws f{rab) ~ with a > 1. This result generalizes to d-dimensional lattices. When 
considering blocks of L particles contained in a d-dimensional ball, Sl can at most grow 
like the volume of that ball, whereas we find that for a > d the upper bound on Sl grows 
at most like the surface of the ball. 



• With the results of the previous section about the Schmidt rank one can also examine the 
special case of simple graph states analytically. Here, the interaction Hamiltonian has a 
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fixed interaction length A and constant interaction strength f{rab) = 1 if Tab < A and 
zero otherwise. One finds for the resulting states I^'tt), i-C-, for t ^ tt, that the entropy Sl 
equals L if the radius of the hypersphere is smaller than A. Otherwise, Sl scales essentially 
Uke the volume of a surface shell with thickness A that is Sl « AL^''^^^/''. 

• One can also consider the dynamics of entanglement that means the change of entangle- 
ment and correlations of the state |5't) with time. The scaling of the entropy with the 
block size L is essentially still governed by the specific form of the distance dependence 
for any finite t, because infinitely remote regions still influence a subsystem in a simi- 
lar way as discussed before. For large times t, more and more of the interaction phases 
'Pab = /(a, b)t start to oscillate (as they are effectively taken modulo tt) and approach 
in the limit of large t a (quasi)-random distribution. In the limit of an infinite chain and 
t ^ oo, the entropy of the reduced density operator of any finite group A of particles is 
maximal, S{pa) = \ This can be seen by considering the off diagonal elements of re- 
duced density operators, which all contain infinite products of cosines of (sums of) random 
angles. All these products tend to zero for N oo, leading to a maximally mixed state. 

Finally we will slightly change the physical 
setup and consider a ip/n ^fli that is a system of in- / ^, / \ \ f 

teracting spins with coupling strengths that are now 
stochastic functions in time Isill . One example is 
given by a semi-quantal Boltzmann gas of parti- 
cles where each particle follows a classical trajec- 
tory but carries a quantum degree of freedom that 

is affected whenever two particles collide. During Fig. 22. - Schematical drawing of a semi- 
a collision of two particles, the internal degrees of quantal Boltzmann gas |5l||. 
freedom interact and can become entangled. This model was introduced in ref. IjtIi and closely 
analyzed in ref. liHl . There an ideal gas was considered in thermal equiUbrium with elastic colli- 
sions, whose mean free path is comparable to the size of the enclosing volume. For a fixed inter- 
action pattern of the particles the time dependence of the weighted graph state |^'t), representing 
the internal degrees of freedom, is given by the adjacency matrix T{t) as a function of time. From 
this the corresponding statistical state can be obtained through an averaging process by assign- 
ing a probability to every collision history and assuming that the collision particles acquire an 
interaction phase ipab inversely proportional to their relative velocities. Thus the weighted graph 
state, corresponding to the resulting statistical szenario, is specified by the density, the volume 
and the temperature of the underlying Boltzmann gas as well as by the initial spatial distribution 
of the particles, which was chosen to be homogeneous. Whereas for low temperatures T the 
initial entanglement generation rate is proportional to Vt, in the case of high temperatures the 
entanglement generation is governed by the slow collision events and proportional to \JT~^ Isill . 
Moreover in the long time limit and for sufficiently many particles the state l^'oo) (equilibrium 
state) is maximally entangled with respect to all bi-partitions (A, i?), i.e., Sa — ^ a- 

In I51I1 another model is also discussed, where the gas particles hop randomly between dif- 
ferent sites of a lattice. This semi-quantal lattice gas allows the (numerical) study of highly cor- 
related classical trajectories. This framework provides microscopic models to compare Marko- 
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vian and non-Markovian as well as correlated and uncorrelated decoherence processes. Non- 
Markovian decoherence (^^) of a set of particles is dominated by (repeated) collisions with the 
same set of environmental particles (high density), whereas Markovian noise is mainly driven 
by independent collisions with different particles (dilute gas). Similarly, correlated decoherence 
processes can be modeled by collisions of some particles with the same particle from the envi- 
ronment. 

We also remark that the fact that for weighted graph states reduced density operators of 
(few) particles can be calculated efficiently can be used to de velop a novel method for ground 
state approximation of strongly correlated quantum systems 116311 . In particular, (superposi- 
tions) of locally transformed weighted graph states can be used to approximate the ground state 
of a strongly correlated system in any spatial geometry or dimension. As expectation values 
of local observables, including energy, can in this case be efficiently calculated, a (numerical) 
variation over weighted graph states can be performed and good approximation to the ground 
state may be obtained. The fact that weighted graph states exhibit rich entanglement features, 
including infinite correlation and entanglement length, maximal localizable entanglement and 
maximal block-wise entanglement, may even allow for a successful treatment of critical systems 
or systems in dimension d> 2. 

10. - Graph states in the presence of decoherence 

Let us now return to the case of simple graphs. Real implementations for graph states are 
subjected to decoherence. Thus the corresponding real state p is in general some mixed state 
that is -depending on the quality of the preparation procedure- more or less close (^^) to the 
ideal pure state |G). First, we study the stability of entanglement in graph states when they are 
exposed to noise and discuss the lifetime of entanglement with respect to the number of particles 
in the system. Finally, a method to overcome the effect of decoherence is introduced, namely 
entanglement purification, which was already mentioned in sec. l4'5l where we also indicate pos- 
sible applications to fault tolerant quantum computation (see also sec. l4 2t and to secret sharing. 
In this section we will focus on one of these applications and show how the purification pro- 
tocols can be modified such that the distribution of the (two-colorable) graph states can remain 
unknown to the different parties {blind purification). 

Let us start with a few remarks underlying noise models that we consider in the following. 
An arbitrary noise process acting on an iV-qubit system that ideally is prepared in the pure graph 
state pg — \G) (G|, is frequently represented as a completely positive map V (CPM) that can be 
decomposed in terms of Pauli matrices acting from left and right: 

3 

(181) p = V{pg) := V E., 

^ — ' ii 

(i = l, -..,«) 



(^^) Here the decoherence of the state for a set of particles is due to collisions with the remaining particles, 
after tracing out the environmental degrees of freedom. 
C^®) E.g. in terms of the fidelity {G\p\G). 



^1 „N „ ^ N 

■"N CTfci ■■■^kM PG^^h "'^In 
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Under certain conditions on the noise process the completely positive maps T) ^ Vt are 

the solutions to a master equation describing the dynamics of the noise process in time. But, in 
the following, we keep our considerations on the level of CPMs, whose Kraus coefficient matrix 
may be time dependent, i.e., Ei-i, = Eti, . [t], if some dynamical description(^^) is 

impos ed. 

In jl65ll it is shown that by randomly choosing some Pauli matrix for each individual particle 
and applying them before and after the actual decoherence process occurs, an arbitrary noise 
process can be 'depolarized' such that the overall noise process is described by a tensor E , , 
that is diagonal, i.e. 

3 

(182) p' = V'ipo) = 

(i = l N) 

We refer to such channels V' that are diagonal when decomposed with respect to the Pauli matri- 
ces, as (iV-party) Pauli channels. For graph states we can make use of the following relations (^*) 



(183) = aMG) al\G) = af"+1G) 

in order to rewrite all a^.- and -operators in the decomposition eq. jl82t in terms of a^- 
operators. In this way one verifies that the resulting state p' in eq. ( I182> is diagonal in the 
graph state basis \U)g = <^^\G) (see Proposition|3jl. In the following we call such states graph 
diagonal states. In order to achieve the simplified standard form for the map describing the de- 
coherence process, i.e., for arbitrary input states, one has to perform some operation before the 
actual noise process has affected the state. Considering noisy preparation procedures for graph 
states, a 'before' does not make much sense and thus the resulting state 'D{pg) can hardly be 
regarded as some standard form for an imperfectly prepared graph state p. But it was shown 
in ref. (3 that the resulting standard form can nevertheless be obtained by performing another 
twirling operation after the (noisy) preparation procedure that gives rise to the same depolariza- 
tion for the corresponding mixed state p js^-tlSO]: 

Proposition 23 (Graph-diagonal states as standard forms for mixed states). Any mixed state 
p of N qubits can be depolarized into some graph diagonal state 

(184) P = >'u\U)g{U\ 

(7cy 

for some graph G with N vertices by uniformly choosing the 2^ stabilizer elements a E S and 



(^^)For details we refer e.g. to ref. 

(^*) These relation directly follow from the stabilizing properties of the correlation operators Ka corre- 
sponding to the graph state \G). 
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applying them to the state, i.e., the corresponding twirling protocol is 
(185) P' ^ ^Y.^'P''- 



For the rest of this section, we restrict to decoherence that arises if all qubits individually (or 
independently) are affected by noise described by the same Pauli channel 

3 3 

(186) V{p) = ^p,a,pa, with (^p, = 1). 



These decoherence models are of particular interest in quantum information theory, especially in 
the study of fault-tolerant quantum computation, and contain for example: 

1 . for pq = ^^j^ and pi = P2 = P3 = ^-j^ the depolarizing channel 

(187) V{p) = pp+{l-p)^l; 

2. for po = Pi = P2 = and = the dephasing channel; 

(188) V{p) ^ pp+^^{p + a,pa,); 

3. for pq = P2 = P3 = and pi = the bitflip channel 

(189) V{p) = pp+^{p + 

The corresponding coefficients Xu are given by the following proposition l49il . 

Proposition 24 (Effect of individual Pauli channels on a graph state). Under decoherence de- 
scribed by the same individual Pauli channel V, the graph state \ G) transforms into a mixed state 
p = riaey T^'^{pg) that is diagonal in the graph state basis \U)g ■ The diagonal elements Xjj 
in eq. il84\l can be computed to be oftheform(^^) 



nam \ - J^l „|(7'\(rc/'+(7)| |c/'n(r(7'+;7)| |(rc/'+(7)\c/'| 

C/'CV 



(®®) In both expressions we have made use of the notational simplifications described in sec. im For 
example TU' denotes both the set and the binary vector that is obtained by the multiplication (modulo 2) 
of the adjacency matrix V with the binary vector corresponding to the set U' . 
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where Qi := ^ for i = 1,2,3. In the case of the depolarizing channel {q := qi = q2 = qs = 
) this simplifies to 

(191) Ac/=pl^l Yl . 

U'CV 

In the next subsection, we will analyze the time dependence of the entanglement properties 
of the decohered state p{t) for different initial graph states p(0) = \G) {G\. Although the results 
of the following subsections can be generalized to more general decoherence models, we will 
restrict our review to the special case of noise that effects each particle in the graph state individ- 
ually. In other words the evolution at each qubit is described by the map given by eq. ( I186> 
with Pauli operators aj acting on qubit a. We will be interested in the evolution of a given pure 
state I ^P) of qubits under this decoherence model. That is, the initial state | ^P) suffers from 
decoherence and evolves in time to a mixed state p{t) given by 

(192) pit) ^V^V^ ...V'^\G){G\ . 

The depolarizing channel with p{t) ~ e^"* is of particular interest, since the decohered state 
due to an arbitrary noise channel can be further depolarized to some state, which might also be 
obtained directly by some depolarizing channel. Moreover, among the stated noise models the 
depolarizing channel is the only channel that is basis independent, i.e., invariant under unitary 
transformations. We will frequently use the Pauli channel and will describe the entanglement 
properties of p{t) in terms of the parameters pi. Nevertheless one has to keep in mind that the 
time dependence itself is already included in the parameters pi = Pi{t). 

10' 1. Stability of entanglement. - Let us now examine the stability of entanglement in graph 
states under the influence of decoherence. As we have seen, multi-particle entanglement is a cen- 
tral property for many practical applications in quantum information. For all these applications 
it is therefore of great interest, to determine the lifetime of entanglement when it is exposed to 
noise. A second motivation aims at a more fundamental direction. With nowadays available tech- 
nologies it is possible to prepare and observe entanglement on microscopic scales; but it is also 
often argued that this task might become exceedingly difficult when considering a macroscopic 
number of particles. So this subsection addresses the question, whether multi-particle entan- 
glement can be stable on a macroscopic level. For the lifetime of entanglement it is not only 
necessary to specify the underlying decoherence model, but also the very notion of multi particle 
entanglement (^") itself. This is mainly due to the fact that multi party entanglement is a subtle 
issue in quantum information theory (see e.g. 167ll ). Apart from some special cases, the ex- 



istence of an entanglement measure that is satisfying for information theoretic purposes as well 



{™) The lifetime discussed in this chapter differs conceptually from the often used Ti - or r2 -decoherence 
(or -relaxation) rates in that the latter are rather related to the stability of quantum coherences and classical 
correlations rathe r than to the stability of entanglement. For some some interrelations between these notions 
see e.g. ref. Il6dl . 
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as applicable and calculable for mixed states, is still an open problem(^^). In the following we 
will therefore concentrate on the discussion of two qualitative entanglement criteria. Through- 
out the chapter we will consider N two-level systems (qubits) with corresponding Hilbert space 
H = (C^)**^. The N particles are distributed among N parties 1, . . . , iV. Starting with a pure 
GHZ or graph state we will consider the A^-party separability and distillability properties of the 
decohered state p{t) (see eq. ( I192h : 

On the one end of the scale the state p{t) can still be N -party distillable entangled, as it is the 
case for the corresponding pure states in question. Hereby we call p{t) A^-party distillable, if any 
other true A^-party entangled state | $) can be obtained (distilled) asymptoti cally from multiple 
copies of p under local operations and classical communication (LOCC) 142114511 : 

(193) ^Locc |*)($| . 

We remark that in the multi-copy case all true TV-party entangled states are equivalent since 
they can be transformed into each other by LOCC. That is, the condition that any true N party 
entangled state can be created can be replaced by the condition that some A^-party entangled 
state, e.g. the initial pure state, can be created. Disregarding the practicabiUty of the underlying 
distillation protocol, the state p{t) is then as useful as any other entangled state and therefore 
can in principle be regarded as a universal resource for quantum information processing such as 
quantum communication. 

On the other end of the scale, p{t) might have also become completely separable or classical 
in the sense that it can be described by a classical mixture of product states, i.e., p is N-party 
separable, if 

(194) p{t)=Y,Pkpl®pl®-.-®Pk- 

k 

If a state is completely separable, it is no longer entangled with respect to any partitioning. In 
between these two extremal cases, p{t) can contain different types of blockwise entanglement. 
We can consider different partitionings of particles into M groups {M < N), where each group 
forms a subsystem with a higher dimensional state space and consists of several particles. M- 
party distillability [separability] can then be defined with respect to a given partitioning in a 
similar way, where the notion of local operation has to be adapted accordingly. We will call p{t) 
M -party distillable, if there exists at least one partitioning, with respect to which p{t) is 7\/-party 
distillable. 

Based on the notion of M-party separability and distillability, one can define the lifetime 
of entanglement. A A^-party state l^*) (^| which is subjected to decoherence for time t evolves 
into a mixed state p{t). The lifetime of A^-party distillable entanglement is given by the time 
after which the state p{t) looses the property of A^-party distillability. This implies that lower 



C^) We note that for systems with only a few number of qubits (A'' < 7), quite recently Il68ll the ef- 
fect of decoherence on GHZ- and W-states was studied in terms of an entanglement measure which is a 
generalization of the concurrence. 



82 



M. HEIN, W. DUR, J. ElSERT, R. RAUSSENDORF, M. VAN DEN NEST and H.-J. BRIEGEL 



bounds on the lifetime of distillable entanglement can be obtained by showing that the state p{t) 
is distillable, while an upper bound can be obtained by proving non-distillability of p{t). When 
considering partitions of the system into AI groups, the Ufetime of M-party entanglement with 
respect to a given partition is defined accordingly. We refer to the lifetime of A/-party entangle- 
ment as the time after which p{t) is non-distillable with respect to all A/-party partitions. In a 
similar way, one can define a lifetime with respect to the separability properties of p{t). 

In order to determine entanglement properties of t he mixed states in question, we will contin- 
uously make use of the partial transposition criterion ll69Hl70ll . an entanglement criterion which 
provides necessary conditions for distillability and separability. The partial transposition is de- 
fined for bi-partite systems only, while a system can in general consist of several parties. Making 
use of the concept of partitionings of the system, in particular considering all bi-partitionings, 
one can use the partial transposition criteria also for multi-partite states. Let A denote a subset 
of TO parties ai, . . . , In general, given an operator X acting on (E) C^^, we define the 
partial transpose of X with respect to the first subsystem in the basis {|1), |2), . . . , I^a)}, X'^^, 
as follows: 

(195) X^- {hk\X\j,l) \j,k){i,l\. 

A Hermitian operator X has a non-positive [positive] partial transpose (NPT) [(PPT)] if X^-* 
is not positive [positive] respectively. That is, X^-* is NPT if there exist some such that 

{^\X^^\^) < 0. 

The positivity of the operator p^-* gives a necessary criterion for separability, whereas the 
non-positivity of p'^-^ is necessary for the distillability of the density oper ator p. In particular, if a 
bi-partite density operator is PPT, then it is certainly not distillable jl70ll . This implies l43l that 
if a multi-particle density operator p is PPT with respect to at least one bi-partite partition, then 
p is certainly not A^-party distillable. On the other hand, positivity of all bi-partite partitions is a 
necessary condition for iV-party separability. In the case of two dimensional systems Cf® 
the PPT [NPT] criterion is necessary and sufficient for separability [distillability] jl69l I7lll . A 



detailed discussion of the application of the partial transposition criteria to multi-partite systems 
can be found in ref. I42I1 . 

Although the computation of the spectrum for the partial transposition of a general den- 
sity matrix requires a numerically demanding diagonalization procedure (^^), for graph diagonal 
states this spectrum can be directly determined from the spectrum {Xu)ucv of the graph diago- 
nal state itself li^ : 

Proposition 25 (Partial transposition for graph diagonal states). For any graph diagonal state 
p il84l the partial transposition p"^-* with respect to some partition A is again diagonal in 
the (same) graph state basis \U)g- In order to compute the corresponding eigenvalues, let 
r' — T/iyie denote the adjacency matrix of the graph between the partition A and its com- 

C^) Note that the size of the spectrum grows exponentially with the number of particles. 
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plement A'^ (see eq. U161 }. Then: 

(196) = ^'u\U)g{U\ with 

ucv 

\kerT'\ ^ ,Ux,Ay) x 

% = ^^^^ \u+x+Y) , 

(kerr')-i-x(Mr') 

where Ay G A is any set with T' Ay = Y, and the kernel ker or the orthocomplement _L are 
taken with respect to the subspace 'P{A) spanned by the sets in A. 

Let us give two examples for formula ( I196> : If V is invertible, then kerF' = {0} and 
(kerr')-^ = V{A) holds. Moreover (IT961 can be simplified by parameterizing Im F' with Y = 
T'A2, where A2 C A: 



(197) ^u = ^ E (-1)^^"^^^ A(c,+^,+r'A.) . 

Ai,A2CA 



If A = {a} for a non-isolated vertex a G V the eigenvalues of the partial transposition with 
respect to A are 

(198) A'j/ = - {Xu + Xu+Na + Xu+a — Xu+N^+a) ■ 

Similarly for the partial transposition with respect to the split A ~ {a, b} versus rest, where 
a,b E y are two non-adjacent vertices with linearly independent neighbor sets Na and TVf,, one 
obtains; 



(199) ^'^ = I f E E ^ 



where 



M+ = {0, a, 6, a + b, Na, Nb, Na + Nb, a + Nb,b + Na, a + b + Na + Nb} and 
M- ^ {a + Na,b + Nb,a + Na + Nb,b + Na + Nb,a + b + Na,a + b + Nb} . 

If a and 5 are adjacent the same formula holds but with neighbor sets N'^ = Na \ b and iV^ = 
Nb\a restricted to A'^. 

This procedure to compute the eigenvalues of the partial transposition described in (iii) does 
not require the diagonalization of a 2^ x 2^-matrix and therefore allows the evaluation of the 
PPT criteria with respect to different partitions, as long as the vector consisting of the initial 
eigenvalues Xjj (which is already of length 2^) is small enough to be stored and -in the case that 
it occurs as a result of Pauli channel- as long as this vector can also be computed fast enough. 
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Consider, for illustration, rings up to size 
= 10 suffering from decoherence due to 
the depolarizing channel, which are examined 
with the help of the partial transpose with re- 
spect to all possible partitions. Fig. de- 
picts the critical value for p, after which the 
state p first becomes PPT with respect to some 
partition, which implies that at this point the 
state p is certainly no longer -party distill- 
able. For fig.|23lthe critical value pcnt has also 
been computed, after which the state p has be- 
come PPT with respect to all partitions, i.e., 
after which p contains at most bound entangle- 
ment with respect to any partition. In contrast 
to the case of -party GHZ states, for which 



Peril 

0.8r 



first PPT 
all PPT 

uppei^ 

uppei^ 
N 



4 5 6 7 8 9 10 

Fig. 23. - For the case of particles in rings of size 
A'^ < 10, which individually decohere according to 
the same depolarizing channel with parameter p: the 
critical value pcrit, after which the first [last] partition 
becomes PPT A [□], the lower and upper bounds ac- 
cording to ref. (431. 



the one-versus-rest partition is the first to become PPT, the numerical results for small N indi- 
cate that in rings this split seems to be most stable against decoherence due to noise described 
by individual depolarizing channels, and that the smallest eigenvalue of the partial transposition 
with respect to these one-versus-rest splits {a} is given by Ajv„+a- 

If the initial state is a GHZ state, e.g. in the star graph representation with center at qubit 1, a 
direct decomposition with respect to the logical basis (|0), |1)) is sometimes more advantageous 
for the computation of the partial transposition and its eigenvalues. In particular, a certain sub- 
class of such GHZ diagonal states, namely the one where Xjj = Xu+i for all [/ C V except 
and {1}, allows for a very precise determination(^'^) of its entanglement properties l42il . 

This class of GHZ-diagonal states is, for example, naturally obtained, if the qubits of the ini- 
tial GHZ state are effected individually by white noise. In this case, most of the abov e entangle - 
ment properties of the decohered mixed state p{t) can be determined analytically i4i [l7i[l7l : 

Proposition 26 (GHZ states in the presence of individual white noise). 

Consider a GHZ state that is exposed to decoherence described by individual depolarizing 
channelsDt with parameter p{t) = e^"*. Then the resulting mixed state p{t) = DY \GHZ){GHZ\ 



C^) More precisely, starting with a pure star graph state \ G) this class of states p is given by 

p = Ao1G)(G1+Ai|1)g(1|+ J2 ^u{\U)g{U\+\U + 1)g{U + 1\) 

ucv\i 

with A = A() — Ai > and thus determined by 2^^^ parameters. In (0| it is shown that (i) any mixed 
state p' can be depolarized to this class of states p by means of random local operations (leaving Ao = A0, 
Ai = A'lj^j and all \u = l/2(Ai7 + Ay+i) invariant), (ii) p is PPT with respect to the bi-partite split 
{A, B) iff A < 2\a\i, (iii) p is separable with respect to the partition {A\, . . . , Am) iff all bi-partite 
splits {A, B) > {Ai, . . . , Am) are PPT, (iv) a maximally entangled state is distillable between the pair a 
and b iff all bi-partite splits {A, B) such that a £ A and b £ B are NPT and (v) a maximally entangled 
Af -party GHZ state is distillable between the vertices ai, . . . , um iff all bi-partite splits (A, B), such that 
neither all ai , . . . , um £ A nor all ai , . . . , om G B, are NPT. 
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is N -party distillable [separable] iff the partial transpose {po) ^ w.r.t. all possible partitions A 
are non-positive [positive ]. More precisely, 

1. p(t) remains iV-party distillable entangled as long as the most fragile splits, 
l-versus-(N — 1) particles, remain NPT. 

2. p[t) becomes A^-party separable as soon as the most stable splits, 

-versus- ^^p- particles, become PPT. 

3. The lifetime of A^-party distillable entanglement decreases as the number of particles 
N ^ oo. 

The Af -party distillability [separability] of the resulting mixed state p{t) with respect to a fixed 
partitioning is determined by the subsystem that contains the smallest number of particles, since 
the corresponding partial transposition is the first one to become positive. In particular: 

4. The maximum number M of subsystems that remain entangled, i.e., the effective size of 
M-party distillable entanglement decreases in time. 

5. As N oo any partitioning into groups of size m leads to vanishing [finite] lifetime of 
the corresponding M ~ ^-party entanglement, if the size m of each group is finite [tends 
to 00 ]. 

Encoding the qubits of a GHZ state by some quantum error correcting code and performing error 
correction at the end of the noise process, the lifetime of [blockwise] entanglement between the 
logical qubits can be increased effectively. Thus the lifetime M-party entanglement in encoded 
GHZ states can also remain finite on a macroscopic level as long as the level of encoding (e.g. 
concatenations) is sufficiently large. 

Thus, we find that without error-correction the lifetime of iV-party entanglement in GHZ- 
states vanishes on a macroscopic level. This result can also be extended to more general deco- 
herence models 1.49.. ,.172:1 . In the remainder of this subsection we now discuss the lifetime of 
A'^-party entanglement in graph states and show(^^) that for a significant subclass, such as the 
cluster states, the lifetime of distillable entanglement is, in contrast to GHZ-states, essentially 
independent of TV. 

To this aim, we establish a lower bound on the lifetime of graph states by considering an 
explicit entanglement distillation protocol. The distillability of a iV-party entangled state can be 
shown by providing a procedure that allows one to generate maximally entangled pairs shared 
between any pair of neighboring particles. This is sufficient, as any A^-party entangled states 
can be then produced by local operations from these entangled pairs. The distillability of such 
pairs serves only as a tool to show A^-party distillability, and no conclusions about the nature 
of entanglement contained in the state can be drawn. In particular, one should not conclude that 
entanglement contained in a cluster state were in some sense only bi-partite. 



C^) An anal ysis o f long-range entanglement in the thermal state of a three-dimensional cluster state can be 
found inref. fTTl . 
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For decoherence of individual particles described by Pauli channels, one can make use of the 
following facts |0]: (i) measuring all but two particles a, 6 in the eigenbasis of a, results into 
the creation of another graph state with only a single edge {a, b} (see Proposition0in sec.|3}; (ii) 
the action on a specific graph state of any operator O which is a tensor product of Pauli operators 
can be equivalently described by an operator O' consisting of only operators acting on the 
same graph state. This implies that a PauU-diagonal map V'^ acting on qubit fc of a graph state 
I G) can be described by a map M whose Kraus operators contain only ct^ operators acting on 
qubit k and its neighbors. This follows from (T^|[/)g = {~l)^''a'^Ka\U)G, where := a%Ka 
is an operator which contains only products of operators at neighboring particles of particle 
a, and similarly for cjy. 

Hence measurements of on all but particles a, h commute with the action of maps describ- 
ing the decoherence process. The resulting state pab of particles a, b is only influenced by noise 
acting on particles a, h and their neighbors A^q, iVf,. The measurements effectively decouple par- 
ticles a, h from the remaining particles. Distillability of pab can be determined by employing the 
partial transposition criterion, where for a two qubit system negative partial transposition akeady 
ensures distillability. As pab is only determined by noise operators acting on particles a, 6 and 
their neighbors, this already implies that for all graph states with constant degree (e.g. cluster 
states), the distillability of pab will not depend on the size of the system N . In fact, one finds a 
threshold value for parameters describing the decoherence process that only depends on the local 
degree of the graph. To be precise, the influence of independent neighbors and joint neighbors 
is slightly different. We remark that in certain cases, e.g. for GHZ states, the local degree itself 
may depend on N . If this is however not the case, the lower bound on the lifetime of distillable 
entanglement obtained in this way is constant and shows no scaling with the size of the system 
N. 

For local depolarizing channel, one finds for in- -in(p„i,) 
stance that pab is distillable entangled if + 
p\N^+Nb\ j^p\Ni\+i > I Solving this polynomial o.e \ 
inequality yields: for GHZ states: Kt > ln(2)/N, in \ 

agreement with the analytic results for GHZ states; for u, 

ID, 2D, 3D cluster states: nt > 0.3331, 0.1886, 0.1318 

0.2 ^ 

respectively. The results are summarized in fig. 1241 V 3" 

J o I 1 GHZ exact 

The results can be extended to any noise model with 4 — . ghz bound 

20 40 60 80 100 N 

some finite range, i.e., all Kraus operators act non- 

trivially only on a finite, localized number of qubits. Fig- 24- - Under individual coupling due to the same 

u»» depolarizing channel the lower bounds on nt to the 

Also m this case, for graph states correspondmg to . ^ .„,,,. , ^ , 

1 ^ hietime oi distillable A* -party entanglement tor the 

graphs with a constant degree there exists a lower bound ^ 2d., 3D-cluster state remain constant for ar- 
on the Ufetime of distillable entanglement that is finite bitrary system sizes N. For the iV-party GHZ state 
and constant, independent of the total number of parti- the lower bound as well as the exact value for nt, 
cles A'^. until which GHZ state remains distillable entangled, 

Several upper bounds on the distillable entangle- -^t^^tly decreases and goes to zero as iV oo. 
ment in general graph states can also be provided 14911 . which again do not depend on the size 
of the system as long as the maximal degree remains constant. These bounds are obtained by 
employing the partial transposition criterion, or by considering the noisy interactions resulting 
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from interactions between adjacent particles in the graph and phase noise acting on them. In the 
second case, an upper bound on the Hfetime is obtained by ensuring that the resulting CPMs are 
not capable to generate entanglement. Considering blockwise entanglement a detailed analysis 
for general graph states has not been accomplished yet. Nevertheless, the scaling behavior of 
M -party entanglement is restricted to a range between the upper and lower bounds, which in 
the case of cluster and similar graph states are independent of the number N of particles. In 
this sense also the scaling behavior of blockwise entanglement in these states must be essentially 
independent of the size of the system. Finally, most of the above results can also be extended to 
the class of weighted graph states (see sec.|5}. 

To summarize, we found that under quite general assumptions about the underlying noise 
model the lifetime of true iV-party-entanglement in GHZ-states decreases with the size of the 
system, whereas for cluster and similar graph states the lifetime of N-party entanglement is 
essentially independent of the size of the system. These results suggest a remarkable robustness 
of certain classes of macroscopic entangled states -namely all graph states with constant degree- 
under various decoherence processes. 



10'2. Entanglement purification. - Although the results of the previous section suggest a 
robust scaling behavior of the entanglement in graph states with respect to system size, we nev- 
ertheless need methods to protect theses states from noise. Many applications make use of pure 
graph states, and hence a method to generate high-fidelity approximations to such pure graph 
states is of significant importance. The reasons that in a realistic situation one ends up with 
mixed states rather than pure states are manifold. For instance, the qubits constituting the graph 
state may interact with uncontrollable degrees of freedom of the environment, leading to deco- 
herence; in distributed scenarios which have to be considered in the context of certain multi -party 
communication settings, the multi-particle entangled states are distributed through noisy quan- 
tum channels. In this section, we will discuss ways to maintain or recover the entanglement of 
such (noisy) graph states. In particular, we will describe multi-party entanglement purification 
protocols that allow one to create, from many identical copies of noisy entangled graph states, 
few copies with increased fidelity. The fidelity can, under the idealized assumption of perfect 
local control operations, be made arbitrarily close to one. Even in the presence of noisy local 
control operations, a significant enhancement of the fidelity is possible. The entanglement pu- 
rification protocols show in fact a remarkable robustness against noise, where channel errors of 
the order of several tens of percent are correctible, and errors in local control operations of the 
order of percent are tolerable. 

In ref. Issll multi-particle entanglement purification protocols for all two-colorable graph 
states have been developed. These protocols are generalization s of a protocol for GHZ-state 
purification introduced in ref. 1 175 1 and further developed in ref. lll76ll . In the following we will 
briefly discuss the recurrence protocols for the purification of general two-colorable graph states 
. We remark that also hashing and breeding protocols are not discussed here, which operate 
jointly on many copies and which allow for purification with a finite yield were developed in 
refs. 1I38II47I1 . Given a noisy graph state p' corresponding to a two colorable graph G, one can 
always transform p' to a standard form p diagonal in the graph states basis corresponding to \G) 
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without changing the diagonal elements, 

(200) p = Y^X^\w){Wl 

w 

where \W) = erf |G) and \w = tv{\W){W\p') = tv{\W){W\p) (recall that {\W)} form a 
basis). The transformation to this standard form is achieved by applying randomly the local op- 
erations corresponding to the correlation operators Ka of \G) 113 ST . For notational convenience, 
we distinguish the two sets A and B corresponding to the two colors in the two-coloring of the 
graph, with indices Wa and Wb, and identify \ Wa^ Wb) = \W). 

We consider a recurrence purification protocol that consists of two sub-protocols, PI and 
P2, each of which acts on two identical copies pi = p2 = p, P12 ■= pi P2- The basic idea 
consists in transferring (non-local) information about the first pair to the second, and reveal this 
information by measurements. In sub-protocol PI, all parties who belong to the set A apply 
local CNOT operations to their particles, with the particle belonging to p2 as source, and pi as 
target. Similarly, all parties belonging to set B apply local CNOT operations to their particles, 
but with the particle belonging to pi as source, and p2 as target. The action of such a multilateral 
CNOT operation on two graph states can be determined within the stabilizer formalism and is 
given by 

(201) \Wa,Wb)\Va,Vb) ^ \Wa,Wb + Vb)\Va + Wa.Vb), 

where Wb + Vb again denotes bitwise addition modulo 2. A measurement of all particles of p2 
follows, where the particles belonging to set A [B] are measured in the eigenbasis {|0)a:, 
of cTj; [{|0)z, |l)z} of az] respectively. The measurements in sets A [B] yield results (—1)^" 
[(—1)^''], with ^a, £ {0, 1}. Only if the measurement outcomes fulfill {£,a+J2beN Cb)iiiod2 = 
Va -which implies Wa + Va ~ 0- the first state is kept. In this case, one finds that the re- 
maining state is again diagonal in the graph-state basis, with new coefficients 

(202) >^Ua,Ub = X! -^^Ua,Wb^Ua,Vb^ 

{{Wb ,Vb)\Wb+Vb = Ub} 

Here K is a normalization constant such that tr(y5) = 1, indicating the probability of success of 
the protocol. In sub-protocol P2 the role of sets A and B are exchanged, and a similar expression 
for the action of the protocol on initial graph diagonal states can be derived js^l- The total 
purification protocol consists in a sequential application of sub-protocols PI and P2. While 
sub-protocol PI serves to gain information about Wa, sub-protocol P2 reveals information 
about Wb- Typically, sub-protocol PI increase the weight of all coefficients \o,Wb, while P2 
amplifies coefficients Xwa-O- In total, this leads to the desired amplification of Ao.o, where the 
fixed point of the protocol is Ao,o 1, i e-, iterative applications lead to the distillation of states 
with fidelity arbitrarily close to unity. 

A numerical analysis of the purification regime of these recurrence protocols was performed 
in ref. Is^. Results analogous to the robustness of entanglement under decoherence are found 
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for the purification regime of the protocols. That is, for GHZ states the acceptable amount of 
channel noise (per particle), such that purification is still possible, decreases with increasing par- 
ticle numbers. For cluster states, or more generally all two colorable graph states with constant 
degree, this threshold value turns out to be independent of the particle number N. Channel errors 
of the order of several tens of percent are correctable. The applicability and performance of the 
protocol was also analyzed when taking noisy local control operations into account. In this case, 
no maximally entangled states can be produced, but still the fidelity of the states can be increased 
provided the local operations are not too noisy. The reachable fidelity and fixed point is deter- 
mined by the noise in local control operations. Again, for GHZ state the threshold value for local 
control operations such that purification is possible is more stringent for higher particle numbers, 
while errors of the order of several percent are tolerable in the case of cluster-state purification 
of arbitrary size O81I . This remarkable robustness of entanglement purification schemes opens 
the way for applications based on high fidelity graph states in real world scenarios. 

103. Multipartite secure state distribution. - Based on the multi-party entanglement purifi- 
cation protocol described in the previous section, a modification of the scheme was found that 
allows one to purify two-colorable graph states in such a way that all (but one) of the involved 
parties do not know — and have no means to learn — which state they are trying to purify. Such 
a process is the basic tool to achieve secure state distribution (see 14811 ') and constitutes a novel 
quantum primitive with possible applications in distributed secure applications. The goal is to 
generate a high fidelity approximation to a secretly chosen graph state among spatially separated 
agents. This is done by distributing multi-party entangled states from a central station via noisy 
quantum channels to local agents, and increasing the fidelity of the states via entanglement pu- 
rification. It is important that at no stage of the protocol, the involved local agents should be able 
to learn about the identity of the state they are processing, and at the end they should be able 
to hand over a high fidelity approximation of this state to end-users (who have placed the order 
for this state). The secrecy in this context is of particular importance, as different graph states 
represent different resources for security (and other) applications, and the end-user may want to 
keep the ressources unknown to the other parties or to potential eavesdroppers. 

The basic idea of the modification is to keep the identity of the state at any stage of the 
protocol secret, and generate effectively states that are -from the point of view of the local 
agents- described by completely mixed states. To this aim, random local basis changes using 
Pauli operations are independently applied to the individual copies before distributing the states 
and hence before purification starts. This ensures not only that initial states are randomized 
(recall that already application of operations at different locations produces all possible basis 
states), but also that measurement outcomes and measurement statistics are randomized. The 
purification protocol can be adopted to account for these additional basis changes. There are a 
number of technical details regarding the exact protocol to ensure unconditional security, and 
we refer the interested reader to ref. jiill for details. We remark that two alternative ways to 
achieve the secure distribution of two-colorable graph states were also proposed in ref. l48il . 
one based on the purification of Bell pairs with subsequent teleportation, the other based on the 
purification of enlarged graph states where randomization (and secrecy) is achieved by additional 
entanglement with some central station. The three different protocols are applicable in different 



90 



M. HEIN, W. DUR, J. ElSERT, R. RAUSSENDORF, M. VAN DEN NEST and H.-J. BRIEGEL 



regimes, where generally direct multi-particle entanglement purification turns out to allow for 
higher target fidelities than schemes based on bi-partite purification. 

11. - Summary 

Graph states form a rich class of entangled states that exhibit key aspects of multi-partite 
entanglement. At the same time, they can be described by a number of parameters that grows 
only moderately with the system size. They have a variety of applications in quantum informa- 
tion theory, most prominently as algorithmic resources in the context of the one-way quantum 
computer, but also in other fields such as quantum error correction and multi-partite quantum 
communication, as well as in the study of foundational issues such as non-locality and decoher- 
ence. 

In this review, we have given a tutorial introduction into the theory of graph states. We have 
introduced various equivalent ways how to define graph states, and discussed the basic notions 
and properties of these states. The focus of this review has been on their entanglement properties. 
These include aspects of non-locality, bi-partite and multi-partite entanglement and its classifica- 
tion in terms of the Schmidt measure, the distillability properties of mixed entangled states close 
to a pure graph state, as well as the robustness of their entanglement under decoherence. We 
have also reviewed some of the known applications of graph states, as well as proposals for their 
experimental implementation. Some of the latter material, specifically about implementations, 
should thus be taken as preliminary and reflecting only the current state of research. 
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